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1. IKTRODUCTIOH 

Very little conclusive Insijitr^ is currently availsble for predicting 
the behavior of reduced-order adaptive controllers (ROAC) , i.e. an 
adaptive controller based on a plant »odel of lower dimension than that 
of the actu.il plant. In [1] the ROAC problem is seen to be an un- 
avoidable consequence of the application of any finite-dimensional, 
lumped-parameter system (LPS) adaptive control strategy to an infinite 
dimensional, dist r'buted-parameter system (DPS). In [2] .»nd [ j] the relative 
orders of the plaiu, model, and controller are noted to be critical. In 
[4] the problem of ad.iptlve parameterization of a C order controller for 
a M order model of a N order system Is addressed when N ■ M > C. However, 
as nv'tcd In 12) and [ll, the practical, poorly understood case is when 
N M ^ C. Tills condition, i.e. the system order exceeds that of the plant 
model used for adaptive controller par.ameterizatlon, will be considered to 
constitute the ROAC problem in this report. 

Tliough numerous strategies have been espoused [5] for extracting a 
low order model from a too complex system description, none currently seem 
fully applicable to the real-time, recursive requirements of on-line 
adaptive control algorithms. These reduced-order modeling techniques seem 
to fall Into two broad categories: (i) Extract those "modes" (or component- 
subsystems) from the lull system description that are most influentlil in 
the performance of the model in Its subsequent use. This strategy is 
followed, for example. In (6) and [7]. (11) Parameterize the reduced- 

order structure to provide the best prediction of the desired output. 

This latter strategy will produce a model the "modes" of which need not 
correspond’ to any of those of the full system as noted in [3]. This latter 
strategy includes the model reference approach prominent In adaptive 
systems [8]. Two methods of interpreting the misbehavior of ROAC, one 
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based on system Input-output description [3] [9] and one on state- 
variable description [1] llO], hove enterged recently from attempts to 
develop adaptive controllers for flexible spacecraft. 

Consider the Implementation of the proper, single-input, single- 
output (SISO), autoregressive, moving-average (ARMA) system 

N 

y(k) ■ I [a y(k-i) b. u(k-i)] (N even), (1-1) 

i-1 


where u is the system input and y the output, in parallel, i.e. partial 
fraction exp.nnded, form as 


Vj(k) 


y,(k-i) + U(k-I)] 


( 1 - 2 ) 


N/2 

y(k) - I y,(k), 
i-1 


(1-3) 


where y^ Is the output of the t-th mode. As Interpreted in I3j the first 
strategy in the preceding paragraph when used for Identification would lead 
to approximation of (1-3) as 

M/2 of [1,N/2] 

“ I y,(k) (1-4) 

with (1-2) describing the dynamics of each yjj^(k). The M/2 modes chosen in 
(1-4) may bo selected as those M/2 from the M/2 In (1-3) that provide the 
"best” fit of y to y. Note that the order M of (1-4) refers to the order 
of the underlying reduced-order ARMA model. If the dimension of (1-4) (and 
(1-1)) were based on the number of quadratic "modes” then M (and N) would 
be replaced by M - y (and ii - ^). This quadratic "modal" designation of 
order is common in the flexible spacecraft literature [3]. As noted in 
111* [31, and [111 the extraction of the yj. chosen in (1-4) from the y in 


(1-3) renuilna an open quuation. Howt^ver, these would be required 
for identif ication of the appropriate (1-2) parameters and 
Alternatively, the second stratcfty would approximate tl~3) with 

M/2 of 11,-M-] 

^(k) - I y,(k) (1-5) 

t 

in, e.g., a least-squares sense by selection of the y^(k). Tlicae estimated 
"moual" outputs need not match the corresponding values in (1-2) and, 
therefore, need not lead to identification of the corresponding and 

in (1-2) as noted in [3]. The successive y ^chosen to fit (1-5) to (1-3) 
may not even obey a time- invariant difference equation of the form of (1-2) 
with replaced by y^. The problem with closing the adaptive control loop 
via a simultaneous identification and control strategy could be intensified 
by feeding back the v. Instead ot the unavailable y. to meet a modal 
control objective. 

Tlie alternate interpretation of the ill-effects of ROAC, derived from 
the flexible spacecraft control problem [1] [10], begins with the separation 
of a SISO, state-space model for, e.g., (1-1) into the reduced-order and 
unmodeled (or residual) segments as 


Xj,(k+1) 


r 

1 


1 

•f 

*’n 

Xj^(k+1) 



1 

Xj^(k) 




( 1 - 6 ) 


y(k) 





XrCR) 


(1-7) 


As outlined in [1] the derivation of (1-6) - (1-7) can be viewed as a 


projection operation on the full systen. rrow (1-7) the ”*pi Hover" 
of the residual modes Into the observation of y(k) via c„ x_(k) and from 
(1-6) the "rtpillover'' of the control designed for the reduced model Into 
the residual modes Xp(k+1) via bj^ u(k) are clearly displayed. Also the 
possible coup I ins of reduced model modes and residual modes via A^^^ and A^^ is 
immediately apparent. Shown in [1] and [10] is the predictable fact that 
if Aj^j^ •• 0 and * 0 then, assuming the residual modes remain stable, 
any full-order (N) , stable adaptive controller identifying A^^, bj^, and c^j 
explicitly or implicitly from y and u alone would retain its stability. 

If Aj^j^ » 0 but Cj^ and bj^ are nonzero then the degradation of ROAC has two 
sources. Unmodeled components in y vi.n Cj^ will generate an error that 
is Indistinguishable trora parameter error thereby causing adaption. Tl»e 
application oi ii>ntrol u to the residual states Xj^ via bj^ will contribute 
further to this unniodeled component of y. Not only will the parameter 
eslimati's bo Incorrect from use of y and u to identify only Aj^, bj^, and 
Cj^ but .ilso the state estimates provided by an adaptive observer will be 
Incorrect, which If fed back could lead to an unpredictable "controlled" 


response. 

These two problems of Inappropriate parameter and state estimation are 
the same ones noted in the first interpretation with y approximation of 
y in (1-5). The ability to extract the y, and obtain y corresponds to 
an effective zeroing of Cj^. This report will. In part, attempt to imple- 
ment, compare, and contrast these two strategies embodied in (1-4) and (1-5). 

The next section details the specific objectives of this study. 

Section III presents the example autoregressive, moving-averiige plants 
that are to be used in the simulations. Section IV presents the adaptive 
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control algorithiM to b«i us«4 and their sources In the literature. 
Section V outlines the forsMts for the slimlated tests Including the 
description of numerical figures of laerlt to be tabulated in Section VI. 
Section VII offers interpretations of the test results and Section VIII 
draws conclusions relevant to the ROAC problem. The laat sections of 
this report include the referenced literature and the appendices 
including computer program listings. 
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II. OBJECTIVES 

The principal objective of this study is to test the usefulness of 
the folklore of reduced-order modelling with respect to adaptive control. 
In particular four "facts" will be tested: 

(i) Heavily damped modes may be neglected relative to more lightly 
danq}ed modes In reduced-order-nodel derivation. 

(ii) Finite bandwidth actuator? limit the mmber of modes necessary to be 
modeled. 

(ill) An "oiitlnuil" reduced-order controller neglects the modes contribut- 
ing the least degradation in the control system performance measure. 

(Iv) Indirect and direct adaptive control are essentially equivalent and 
interchangeable . 

The Implication to the ROAC problem of each of these statements will be 
developed in the following paragraphs. The simulations of the following 
sections will be chosen to test the veracity oi* these implications. The 
conclusion of this report will summarize the useful "facts" that either 
escape unscathed or emerge from these tests. 

Since each of these facts has been accepted into the reduced-order 
and/or adaptive control folklore it is difficult to pinpoint particular 
references succinctly stating these points. However, several classical 
control texts contain the source of point (i) in the concept of dominant 
roots or poles. For example: "Ihe complex conjugate roots neat the 

origin of the s-plane relative to the other roots of the closed-loop 
system are labeled the dominant roots of the system since they represent 
or dominate the transient response. The relative dominance of the roots 
is determined by the ratio of the real parts of the complex roots and 
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will rttsuU in reaaunabls do«inance for rctloa «xc««dins five. ... 

Dominance .. also depends upon the relative magnitudes ... of the 
residues evaluated at the complex roots, [which] depend upon the loca- 
tion of the xeros in the s-piane" [12] . Or: "The relative dominance 
of closed-loop poles is determined by the ratio of the real parts of 
the closed-loop poles, as well as by the relative magnitudes of the 
residues evaluated at the closed-loop poles. The magnitudes of the 
residues depend upon both the closed-loop poles and xeros. If the ratios 
of tlie real parts exceed five, .and there are no zeros nearby, then the 
closed-loop poles nearest the Jui axis will dominate in the transient- 
response behavior bec.auso these poles correspond to transient-response terms 
which doi .tv sl«’w!\" (I I, ;>. .’'^>11. This scp.ir.at Ion i tuicopf h.ts itcen fi>rnjii- 
llzed via singular pert urh.it ion theory [14]. Despite the concomitant 
warning in both (12] .md 111] for caution in the use of this rule and the 
explication of its application to the closed- loop system, point (1) is 
commonly (though admittedly Inappropriately) used for model reduction 
prior to control design. TTie (mis) implication for discrete systems is 
tliat if the open-loop slngulailtles are separable Into a group in the 

5 

z-plane outside a radius of r (<1) and the other group inside the radius r 

Pa T a 1 T' 

(since |z|»e ^ ^ , where T is the sasqile period) then the first group 

alone provides a highly accurate input-output model of the entire system 
and therefore provides a useful dimension and parameterization for reduced- 
order controller design. One objective of this study is to test the useful- 
ness of this guideline for ROAC. Predictably, such a rule will be valid 
only when the closed-loop system retains its open-loop singularity 
separability. 
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A frequently voiced (througl. uiwritten) crlticie* of the pre- 
occupation with the effects of spillover Is based on the assertion that 
the frequency content of the input to a plant has a strong effect on 
the most suitable reduced-order isodel. Since actuators do not have 

infinite bandwidth and are commonly isodel ed as essentially loy pass 

« 

filters [15] the higher frequency modes of the system will receive such 
iuRlgnlf leant excitation that they are ignorable as noted in point (il). 
This again extrapolates a reasonable open^loop response mechanlsia to a 
closed-loop situation. Again, if the requirement that plant open-loop 
singularities were shifted only slightly by the feedback were included 
then this proposition would be strengthened. As such it narrows the 
more classical concept by utilizing the shape of the input spectrum in 
addition to the plant singularity constellation. (>ne problem with the 
use of this idea in the ROAC problem is the nonlinear, time-varying 
character of the feedback, which does not fit the linear system 
Character of this guideline. Tliat is, during adaption a large, 

“high f requcncv** ci>ntrul effort ct'tild imbalance the rull**off provided 
by the low pass actuators* Another problem Is the unmodeled phase shift 
induced by the actuators unless the actuator outputs are available as 
the identifier inputs. Due to the prevalence of this seemingly untested 
proposition, construction of a meaningful example for its examination 
represents another report objective. The influence of input frequency 
content on reduced-order controller selection will also be tested for 
various reference signal frequency content distributions. 
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Point* (i) and (ii) r*ly on the near-*qulval«nce of the two reduced- 
order modeling etrategie* interpreted Ir introduction J modal «electlon 
and full-behavior approx Ima tiers. The euggcstlon in point (iii) recognize* 
the diet inct ion between theee two etrategie* and repreeent* the 
*en*lblc result of pursuit of the first strategy embodied in (1-4). Also, 
clearly the modes of the plant may not be ; a separated as required in the 
discussion of point (i)> either in the open- loop or closed-loop system, 
which requires a selection mechanism. As shown in [7] for lightly damped 
systems forced by infinite bandwidth inputs the modal costs arr> pro- 
portional to the product of the modal time constant, observability norm, 
and dlsturbability (or controllability) norm. Since for a partial fraction 
expansion the obse v 'blJ ity and diaturbability norms are related to the 
njod.ll residual, this augget^tion c.in be viewed as a more sophisticated 
version of point (i). However, as demonstrated by example In [7] the low 
frequency inodes need not always provide the best open-loop reduced-order 
model. Again iUIh mod.il selection procedure is principally Intended for 
an open-loop fit. For reduced-order control .ippllcatiou the truncated modes can 
be viewed as the desired degrees of freedom omitted from the controller. However, 
a misinterpretation of this procedure would suggest that a reduced-order 
controller composed of a prespei ifled number of well-selected modal controllers 
would be optimal. This would suggest (falsely) that the second ROAC strategy of 
full approximation in the Introduction, as represented by (1-5), could 
never prove better than the first of selective extraction, as represented 
by (1-4). The third objective of this study will be to test this impli- 
cation by comparing single-mode controllers chosen to control single-mode 
reduced-models of a two mode system and a single-nvjde controller chosen 
to control the full system to "match" a single-mode objective. Such a 
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coa^arison will clearly rely on the cxa«|>la choaan and la axp«ctcd»as 
with the flrat two objectlvca, to yield initially aabigitoua results. 

Consider for example a second-order systcsi controlled by a constant 
output feedback gain to meet a first order response Matching the dominant 
pole of the root locus. The appropriate gain need not be provided by 
either separate nodal controller but a close fit to the objective does 
exist fur some reduced-order controller paranetcrization. 

In the special case of full-order model following, indirect and direct 
adaptive control have recently been shown to be equivalent (16) [17], 
Indirect adaptive control uses the current plant parameter (and state) 
estimates to solve for the controller parameters (and feedback signal). 
Direct adaptive control updates the controller parameters such that the 
control system response matches that of a prescribed model. Discovery of 
this equivalence has led to highly touted claims of the interchangeability 
of the algorithms resulting from these two approaches. The implication is 
that the two approaches are also equivalent for ROAC. Comparison of the 
indirect self-tuning scheme with the initial plant MA gain known, the direct 

input matching scheme, and the approximate direct output error identification 
interpretation with ■ fj in [17], all of which could be based on an equa- 
tion error parameter estimator, does not dispute this claim for the specific 
model-following problem and equilibrating choice of adaptive gains. If 
the initial plant MA gain is ulso estimated then indirect and direct schemes 
can be expected to behave differently. Only recently has a simple provable 
direct adaptive control scheme, based on output error parameter estimation, 
been developed [18]. Like the approximate strategy in [19] a particular 
choice of designer selected constants reduces this strategy to the common 
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equation error baaed scheme in jl7). The complex j^neratlon of an additional 
reference model input required in [2U] is avoided in [18]. 

As an aaide, note that even though the reduction of a gradient based 
solution [21] and a stability theory based solution [22] to the equation 
error formulation of ARMA process identification prove identical this does 
not suggest the equivalence of gradient based [23] and stability theory based 
[24] solutions to the output error identification problem. At least currently, 
they are seen as special cases ofageneral, possible solution [25]. These 
two forms of adaptive parameter estimation underly the various prov )le 
adaptive control schemes suggesting differences for ROAC. Furthermore 
the distinction of equivalent convergence points versus convergence 
paths must be made. For full-order adaptive identification the asymp- 
totic results may be similar though the transient behaviors are quite 
distinct. This transient disparity Is amplified in reduced-order 
identification where the convergence points also become distinctly 
different [^<>|. Therefore, it would appear that the parameter estimation 
formulation and the explicit or Implicit adaptive control strategy both 
lead to different ROAC behaviors. The fourth ol,jectlve is to test the 
equivalence (or disparity) of Indirect and direct adaptive controllers 
based on gradient and stability theory formulated equation and output 
error parameter estimators. 

Achievement of the preceding objectives may provide suggestions for 
Improving the presently available full— order adaptive controllers for 
ROAC use. Fulfillment of this hope is the ultimate objective of this 


report. 


IIl.l 


III. 

Six test exanplcs will be used to Aset tlie objectives stated in the 
preceding, section. They represent the following six categories: 

(i) Open-loop and closed-loop (in near satisfaction of desired nodel- 
followlng objective) system singularities are separable on a time 
basis due to significant damping differences. 

(11) Open-loop system singularities are time separable but closed-loop 
singularities are only marginally separable. 

(ill) Neither open-loop nor closed-loop singularities are time separable 
but a near fit to desired beliavior exists due to a lightly damped 
nenr-c.mcel lat ion. 

(Iv) Neither open-loop nor closed-loop singularities are time separable 
and no reduced-order controller gain closely satisfies the desired 
control objective. 

(v) Two nearly oscillatory open-loop modes of dlctinct frequency, which 
are only moder.atcly shifted in the closed-loop, are preceded by an 
infinite bandwidth actuator. 

(vi) T«. nearly oscillatory open-loop modes of the same distinct frequencies 
as in (v) are preceded by a low pass actuator providing at least 
-12 db (75%) attenuation at the resonant frequency of the second 
mode relative to the resonant frequency of the first. 

Categories (i)-(lv) will test points (i) and (ili) of the preceding section. 
Categories (v) and (vl) will test point (li). All of the categories will 
teat point (iv). 

Consider, for the first four categories, controlling a stable second 
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order plant with transfer function 

« b ^ e ^ ( b-fe)z - (ba^ +ea^) 

Ut*) z-a^ z -82 ''T*-ajy(s-a2y 
and therefore difference equation description 

y(k) - (aj+a2)y(k-l) - (8ja2)y(k-2) 


(3-1) 


+ (b+e)u(k-l) - (b 82 +eaj^)u(k- 2 ) 


(3-2) 


and D.C. gain 

YU)' 

U(z) 


. _b__. + 

z -1 ^“®2 


(3-3) 


under the assumption of a first order model for (3-l)wlth the objective of 
following a first order model 

s(k) =• cr(k-l) + d.s(k-l) . ( 3 - 4 ) 

Using till* projection technique of (l-h)-(l-7) , (3-1) can be rewritten as 

u(k) 

(3-5) 


(k+li 

m 

‘"1 

0 “ 


Xj (k) 

+ 

V 

X2(k+1) 


0 

“ 2 . 


X2(k) 

1 

r 

L- 


y(k) 


U 1] 


Xj(k) 

X2(k) , 


(3-6) 


where, from (1-6)- (1-7), - a^, - 82 * A^^j^ « Aj^^^ « 0, b^^c^ « b, and 

bj^Cj^ “ e. Note that any proportion of c could be distributed to bj^ and 
Cj^. Therefore Increased c corresponds to Increased spillover. Clearly 
when e » bj^Cj^ ■« 0, a first order model of (3-1) would be exact. In such 
a case, as noted earlier, any stable adaptive control scheme would be 
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successful. Also, equally apparent is the degradation of the reduced- 
order control systsa, adaptive or nut, as c becoaes nonsero and the 
residual mode X 2 contributes significantly to y. 

If e were zero, then control of (3-1) via 

u(k-l) ■ gr(k-l) + fy(k-l) (3-7) 


where 


g 

f 



would convert the plant output to 


(3-8) 

(3-9) 


y(k) ■ ajY(k-l) + bu(k-l) cr(k-l) + dy(k-l) . (3-10) 

Therefore the model-lollowlng error 


s(k) - y(k) - dla(k-l) - y(k-l)] (3-11) 

would decay to zero if (3-41 were stable. If t were not zero, use of 
(3-7) would convert (3-1) to 

yjA) . 1 

lT(z) " ^|{z-ap (z-a 2 > - f((b+f)z - (ba 2 +taj)) j 


(•’•Tj) (”-'P2^ ’ 


(3-12) 


the roots of the characteristic equation of which, i.e. p^ and P 2 » could 
be determined via a root locus. Note that f has a limited effect on the 
poles of (3-12) while g can be chosen to select the U.C. gain. 

Alternatively the control of (3-1) (or (3-2)) could be chosen for 


(3-5)- (3-6) as 
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u(k) ■ Rr(k) + IXj(k) , 


0-13) 


If the reducecl-o> der nodel state la aBSuned available. This asHuauption la 
equivalent to asBuming the poBsibillty of »ea«urlng yj,(k) in (1-4) in the 
fimt reduced-order model strategy in the introduction. The control of 
(J-13) can be vltwed as p .t 1 state feedback (Impractlcally requiring 
state availability) 


u(k) * gr(k) + (f 0) 


(k) 

X2(k) 


(3-14) 


From (3-3), ( 1-b) , and (3-14), (3-13) converts (3-1) to 


-.11 M 
R(7.) ^ ^ 


zl 


ti j 0 


•*2 


[f 01 


-1 


f (b+i )z - (,ba .,-f I a j ) 1 




(z-q^) (z-q 2 ) 


(3-15) 


^ |(z-a^-bf ) (z-a 2 ) 

Note that relative to (3-12) only one pole is arbitrarily shifted, since 
rather than both being shifted under the root locus constraint. Note that ii. 
cither case, (3-12) or (3-15), c can be chosen in constructing the test 
examples for a particular open-loop (and closed- loop, due to pole-shifting 
only by the controllers) zero location h 


ba»+ca, 

h - ~ =P' c 

b+E 


b(a2-h) 

h-a. 


(3-16) 


_Exarap2e_i: a^ - 0.95, b - 0.065, a., - 0.2, r - 0.01, h - 0.3, c - 0.4, 

d - 0.8 


Recall that from (3-1) a^ and a^ are the open-loop pl.int poles, from (3-4) 
d is the desired closed-loop pole, from (3-12) and (3-15) p^, P 2 and q^. 
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are the closed-loop aystetB poles due to output and partial atate feed- 
back, reapect Ively, and from (i-16) h la the open- and closed-loop system 
zero. The closed-loop migration due to various f in (3-7) and (3-lJ) ap- 
pear in Fig. 3-1 and 3-2, respectively. Note that, alnce a^ * 0.77 > h > 82 
and using (3-7) with f “ -2. A ^p^ * ( 0 . 8 )^ • 0.33 ■> h > P 2 " 0.17 and using 
(3-13) with f » -2.3 *>q^ - ( 0 . 8 )^ » 0.33 > h > q 2 ■ 0.2, this example falls 
in category (i). Note the near equivalence of the f’a from (3-7) or (j-13) 
for the same objective. (Refer to Appendix A for the program listings and 
tabular output of these and the following root locus plots.) 

■ 0.9, b ■ 0 . 1 , 82 “ 0 . 1 , e -0.01, h « 0.011, c * 0.7, 
d - 0.65 

The closed-loop root migration due to various f in (3-7) and (3-13) appear 

in Fig. 3-3 and 3-4, respectively. Note that =0.59 > a,, > h and 
5 

d ‘ 0.12 > 32 ^ h, i.e. the dominant open- loop pole and its desired 
location both dominate the other open-loop singularities, but for f - -2.4 
in (3-7) the second closed-loop pole P 2 = 0.13 is crossing the dominance 
threshold with respect to p^^ ■ 0.65 as shown in Fig. 3-3, thereby putting 
this example in category (ii). To retain the open-loop dominance the 
partial state feedback of (3-14) in Fig. 3-4 must be used. 

Exa mple 3 ; a^^ » 0.9, b ■ 0.08, 82 ■ 0.8, c « 0.02, h » 0.82, c * 1, 

d - 0.5 

The closed-loop root migration due to various output feedback f in (3-7) 
is plotted in Fig. 3-5. Figure 3-6 shows the closed-loop roots for various 
partial state feedback gains in (3-13). Reversing the designation of the 
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Example 3. 
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isportant »ode for partial atata faedback by althar swapping Che alaiaancs 
of the state feadback gain vactor in (3- 14} or swapping the designation of 
Sj, b and S 2 . t results In the cloe«d-loop roots ohown in Pig. 3-7. Mote 
In Fig. 3-3 for f * -3.8 the near-cancel 1st Ion of the slower node while 
the faster node natchus that desired. The degree of closeness of this near 
cancellation is not asitthed In either Fig. 3-6 or 3-7. Note that In Fig. 
3-6 if f Is chosen as -S to cause qj • 0.3, uae of the sane f in Fig. 3-5 
would result In «»vercompenMst ion. Neither in the open-loop plant nor for 
any f placing one pole ne.ir the desired location of d • 0.5 are the 
singularities sep.irable as in category (111). 

t **1 " 0.‘^t h ■ 0.5, ■ 0.7, t -0.3, h ■ 0.4, c ■ 0.8, 

d - 0.6 

Hie pole migration due to the output feedback f from (3-7) displays in 
Fig. 3-8 the oscillatory character of the second-order poles in the 
region of desired pole radius. Even partial state feedback proves un- 
successful us shown in Fig. 3-9 due to the retention of one of the poles 
at its nearer unity value than the desired pole. Therefore, this example 
falls into category (iv). 

Consider, for the last two example categories, digitally controlling 
a stable fourth-order discrete equivalent of a continuous plant composed 
of two quadratic, all-pole, modes 

Y(z) ®11^ ^12 . ^'21* ^22 ,, ,,, 

U(V) " "2 * J * 

^ “"ll^"“l2 ^ "“21*" “22 

where the a and 8 are derived from the zero-order-hold equivalent 
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Fig, 3-7: Closed-loop Root Migration due to Partial Stati Feedback 

for Example 3, 
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^11^ ^12 


« (1 -k’’^)Z 


_ 

• +2CjW^8+w^ j 


(3-18) 


with (27] 


f / i] 

tt ^2 ■ 2 e cos L jT/ 1 - c“ 


12 


-e 

X, 


-2CiU,2T 


(3-19) 

(3-20) 


^11 - "2 f ^ ^ po8(iUjT/l^;j) + 

,/~2; 


sin (UjTvi-Cj^) 


(3-21) 


A r -CjU) Tp-; ^ T ^ /-— 

®i2 “ “2 I® ■ co8(u,jT/l-cp + (Cj/Zl-Cj) 

“i ^ 


8in(u)jTZi-t^) 


(3-22) 


The reduced-order, model-following objective is to track the output of 

. (3-23) 


l(z) 

R(z) 


+ 6, 


^ -V''2 


where 7 ^^ and represent the discrete equivalent of desired s-plane 

pole locations In terms of 4 find m and translated via conversions similar 
to (3-19) and (3-20). For ^21 “ ^22 ” ^ objective can only be achieved if 

the numerator of (3-17) is cancelled and replaced by that of (3-23). There- 
fore assume that the numerator zero of (3-23) matches that of (3-17) with 
^21 * ^22 “ ° and only the poles require shifting. This converts the objective 
model in (3-23) to 

S(z) '^^11^ + ^2^ 


R(z) 


(3-24) 


z -Y 2 Z-Y 2 


The control effort 


IH.17 


u(k; ■ 6r(k> + rij^u(k-l) + ri2u(k-2) + Vjy(k-l) + V2y(k-2) , 

(3-25) 


where [28] 




-1 

0 

0 

0 

-1 

“irn 

r\2 


“11 

-1 

-®11 

0 


“l2"^2 

^1 

“ 

“12 

“11 

'hi 

-^11 


0 

_^2 


0 

“12 

0 

"^12 


0 


(3-26) 


converts (3-17), with ^21 “ ®22 * ® (3“24). An algebraic, sequentially cal- 

culated solution of (3-26) avoiding the matrix Inversion Is [3] [9] 


“ ^1 - “ll 


& 1°12 <^i”ii - °12 

pll“ll ®12 " “12^11 '^^i3 


H2 » + Y2 - ^ll\ - “12 


V2 = n2CXj2/Bj2 • 

If P21 *^29 ^ 0 «'ind (3-25) Is used on (3-17) then 

rY(z)/U(z) 


(3-27) 

(3-28) 

(3-29) 

(3-30) 


Y( 
R(z) 


•’ F .//1 -I - 2 x 
. j 1/(1 - n,z - n.,,z ) 




|l - + V2z“^)/(1 - Hj^z"^ - n^z ^) Y(z)/U(z) 


~“ 21 ^~“ 22 ^ ■<• ^( 3222 + 622 ^2 

« 


2 2 2 
(z -nj^z-ri2)(z -a^^z-aj2)(z ~“21^~“22^ 


(V 1 Z+V 2 ) [(SjiZ+3j 2)(2 ■“21^~“22^ *'■ ® "“ll^' 


+ b^z^ + ^ 2 ^^ ^1^ ^0 

_ __ __ _ _ 

a,z + a^z + a,z + a_z + a«z + a, z + a^ 

6 5 4 3 2 1 0 


'-“ 12 II} 


(3-31) 
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where 


b5 - 6(3ji 4- ?2i> 


"■ ■*■ ^12 “ ®2l“ll ®22^ 


“ ^^~®Il“22 “ ^12“21 " ^2l'^12 " ^22“ll^ 


- '^^~^12“22 ~ ^22“l2^ 


bj - bg . 0 


‘*6 " ^ 


(3-32) 

(3-33) 

(3-34) 

(3-35) 

(3-36) 

(3-37) 


^5 “ - “11 " "l 


(3-38) 


-“22 “ll“21 - “12 ^“21 ■*• ’ll"!! " ’"2-''l^l ‘ ^ 1^21 

(3-39) 

“ll“22 ■*■ “l2“21 '^l“22 - '’l“ll“21 ■*■ '^l“l2 ’^2“21 '^2“ll 


^ ''l^ll“21 -''l^l2'*’''l^2l“ll- ''1^22 - '*2®!! - '^2^21 (3-40) 


a„ = “i 2“22 - ’^l“ll“22 ’’l“l2“21 '’2“22 “ '^2“ll“21 ^2“l2 


+ Vj^8^ja22 + '’if^l2“21 ''l^2l“]2 ''l^22“ll 


''2^1l“21 - ''2^12 '’2®2l“ll - ''2®22 


(3-41) 


<l-| “ ^1^12“''2 — '3')®il^oo ~ ^9^19^91 ^1^19^99 3" ^’1 ^00*^ 


11 


2‘*12“21 ''1^12^22 ‘1'22"‘12 


+ V22i 1<^22‘^ ^2^12“2l'*‘ '2®2l“l2 ''2®22“ll (3-42) 


^0 “ ~'’2“l2“22 '*■ ^2^12“22 '^2^22“l2 


(3-43) 
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Note that haa the aaiae numerator aa (See Appendix B for the 

R(z) U(*) '* 

Bupportlng algebra .and alraulnted check for (3-31)-(3-43) for the fol- 
lowing example.) 

E_xamj>J_e 3: T - 0.5, - 1, - 1, “ 0.2, u.j - 0.5, - 0.02, 

0)2 ■ 0.5, Yj “ 1.687, ^2 *' ”0.741. 

The plant In (3-17) resulting from this parameterization has the pole- 
zero pattern shown In Fig. 3-10. Note * 0.05, l.e. both modes 

have the same settling time and are therefore not separable on a time 
basis. However, they are clearly separable on a frequency basis. The 
Cj, and T were chosen to avoid aliasing, retain frequency separability, 

and mtidel the lightlv damped situation of flexible spacecraft [3](10]. 

Further, note tli.il the low frequent y mode has a DC gain which Is 100 times 
that of the high frequency mode. Attempting to increase the damping ratio by a 
factor of three to 0,6 for the low frequency mode (.and therefore third the sett- 
ling time with Uj uneh.mged) leads to the stated obleetlve, which from (3-26) 
or ( J-27)-( 5- 30) yields the tol lowing controller parameterization for 
(3-25): vS - 0.5, - -0.159, n., - -0.0715, » -0.A75, and V 2 * 0.556. 

Factorization of (3-31), where b^ « 0.952, b^ * 1.188, b,^ •• 1.128, 
b., “ 0.83, bj * ‘^5 “ ”0.164, ■ -0.891, a^ “ -0.466, 

a,, “ 0 . 676 , . 1 ^ « 0.0625, and a^ = -0.0338 from (3-32)-( 3-43) , yields the 
pole-zero constellation in Fig. 3-11. (See Appendix C for the factoriza- 
tion and plotting routines used.) Note the closed-loop retention of 
frequency separability. Simply changing X 2 to 10, thereby reducing the 
DC gain ratio from 100:1 to 10:1, leaves the plant poles of 0.923 + j0.23l 
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«nd “0.762 + JO. 570 unnltered but chanRCii the seroa fro« -O.'lbO and -0.1A4 
+ J0.9A2 to -0.959 and 0.675 + JO. 671. Since ^2 does not affect the cal- 
culation of Ojj or Sjj in (3-19)- 0-22), the and selected In (3-27)- 
(3-30) will be the sane. However, this control law results in an unstable 
system. (See Appendix C.) Therefore observation spillover can be seen to 
destabillte an otherwise stable uystem If neglected in reduced-order con- 
troller design. 

The addition of a low-pass actuator in controlling (1-17) would occur 
in the feedback loop after the dynamic output feedback element of (1-23) 
(with * 1), if the reference signal were considered an unmeasurable 
disturbance. However the scaling of r in (3-25) suggests a command signal 
designation of r. Therefore the actuator would precede the pl.mt in the 
forward path. A remaining question centers on the availability of the 
.ictuator output for *'ontrol and identification. Tlte three iw.nnlngful 
possibilities are diagr.immed in Fig. 3-12a-c, where the identifier feeds 
back the parameter irat Ion of the plant for adaptive controller parameteriza- 
tion. Retention of the controller dynamics from (3-25) suggests usage 
of a or c. Tlte frequency limited Identification concept of "fact" (ll) 
in the preceding section implicit in example category (vl) prompts usage 
of b or c. Therefore c, if physically feasible, appears to be the prime 
c.mdldate. However a is the most reasonable if the actuator models a 
physical limitation rather than a control logic device. 

For the simple low-pass actuator model in (15] the zcro-order-hold 
equivalent Is [29] 


H(z) 


(1-z 


-I 





(3-A4) 


Therefore 



(a) Actu.ntor Incorporated in Plant 



(,b) Actuator Incoi poratcd In Controller 



1 


Figure 3-12; Actuator Inclusion Possibilities 
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u(k+l) • u u(k; 4 ,1 - u^(k) , 


when* from (J-25) for Hr. J-I2a and c 


(3-45J 


u.(k) • Ar(k) + n,u (k-l) + n,u,(k-2) + v yU-1) + \ ,y(k-2) 

i i 4* t I 

(3-46) 

and for Hr. 3- 1 2b 


u ?k) 6r(k) + n,u(k-l) + niU(k-2) + v,y(k-l) + v,v(k“2). 

(3-47) 

Applying {3-4'))- (3-46) to (3-1/) with tho luntrollor p.ir.imo tori rat ion 
ariHing from o.g, ( 3-2 7)- ( 3- 30) , i.o. the reduced-order controller Ignoring 
the actuator incluHion, yielda the following over.ill trannfer function 
(.IS shown In Appendix D) 


where 


Y(z) 

R('z ) ' 7 


') , 4 3 2 

‘ItjZ +' 14 ?' + ‘l,z + qjZ 'll* + '»o 


') 4 


z +P^.: +p^z fp^/ +PjZ +p^z +Pj'' + Pf, 


^3 

^12 


HI 

<^(1 

HI 

<5(1 


e“''’’^)tPjj + 


-cT 


- e )(-PijU^l + + e,,) 

e ■ ^12“21 “ ^2l“l2 ” ^'’>’* 11 ^ 


e ^)(-Sj 2<^22 " ® 22 ® 12 ^ 


11 


qi - qo * 0 


-cT 


(5-48) 


11 11 '1 


(3-49; 

(3-50) 

(3-51) 

(3-52) 

(3-53) 

(3-54) 
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- c + «ijj + Hj) - ^ “n'*21 ■ ®12 ''l®21 


+ njOjj - nj 


(3-55) 


• -« (”>122 ‘*ll‘*21 " **12 ■*■ ■*■ " *^2^ * ‘*ll“22 


*'^12^21* ^1^22" ^l‘*n“2l'' Vl2-" V2I* Vll 




(3-56) 


P 3 ' -•’"'^<“ll"22'^°12“21-*-'^l'‘22" Vll“21^''l“l2^V21^Vll^ 


■*■ ‘*12*22 " '^1*1 1*22 “ ’’1*12*21 ■*■ V22 “ ”2*11*21 ■’■ ”2*12 


^1 ^“^11*21 ^ *^22 * ^\2 ** ^21*11^^^ “ ** ^ 


” j 1 ■*■ ^ > j ' ( ‘ " *‘ * 


n-57) 


->'T 


f ('*ij'*22 “ ''1*11*22 " ” 1*12*21 * ” 2*22 ' ”2*11*21 ”2*12 


” *'1*12 *22 ~ ”2*11**22 ” ’’2*12*21 *1 1 1 ** > * 12*21 


— o T 

*1* O — P ) — V^(— *4* ^yyry B 


11 


11 21 ^22 L 


. , 1 -'’Tn 


(3-58) 


Pi - f ("i^' 12*22 ■*■ ”2*11*22 ■*■ ”2*12*21^ “ ”2*12*22 
+ V2(8jia22 ‘"'12*21 ^21*12 ®22*11^^^ ~ ® 


-*oT 

+ v^(t»j,c\,T + 8.,,u,,)U - e ) 


22 12" 


(3-59) 


•"OT *“OT 

Pf^ * p ”2*12*22 ^2^^ ” '■* ^^^12*22 ^22*12' (3-60) 


SimfliJrly tin* ovomll transfer function resulting from application of 
(3-A5) and t3-A7) to (3-17) is 
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Y(.) * ’ 2 *^ * ’l' * ’O 


4 


Y ( 3 - 61 ) 


z + + P 4 *’ + P 3 *' + Pi*’ ■*■ Pi* ■*■ PQ 


where the are aa in ( 3 - 49 >- ( J- 53 ) and 

-oT 


’'*21 “ ^11 ” 


(3-62) 


^ o TT * ^ X 

p^ - -.122 + “ii‘Vi - «i2 ^ ® ^‘*21 + «ii^ ” '’l^’ “ ) 


(3-61) 




|'« ■ ’n'’22 + '’u"2I • •■ ‘•'22 * ■■u''2i ■ ''12’ 


+ •Ijll - +.. I,) - .ij<l VjU -."'’’'’ hSjj tBjj) (3-64) 

*** o *** 

Pj • ^'’11*22 ‘*12 '2p - Hjd - e H-‘»22 '*ll'*21 


Xj ,) - n,U - e j ) - VjU - e ' 


21 n 1 


--T, 


" ■*■ ‘'12 “ ■ *■ ^ 


ii-bb) 


•^T ,1 -‘’Tw ^ 4 

P> * “'•■ *12'*22 " '*1^ ” *’ ^^"‘ 11**22 * '*12“21^ 


1 


ti,U - t ' K-^.,2 + ~ + Vjd - e )(6,,a 


1P22 


(tX 

+ Pj2‘*21 * ‘^2l“l2 '*■ ‘'>->'*11'* “ ■ ** 


ir 


1121 


+ 1^22 + f^i2 ■ '*2l“ll’ 


(3-66) 


?1 • -n^d - ^"°^)'’‘12°22 " ^ '^'^^^”ll®22 * “l2”21^ 

+ v^,d - e '(fcj^a^o + Sj2«21 ^ *^2l'*12 ’’’ ‘'22‘*11- 


+ v,d - t*"'’^dej2‘»22 * 


‘12^ 


(3-67) 


0 


*" X o X 

-n,d - e ' )a^2‘^, , + v.,(l - e )(f, , 5) 


(3-68) 
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# 


T ■ 0.5, Vj " X, I, • 1* Cj 
t«2 “ ■ 1. A • 0.5, -jj 


0.2, h*j ■ 0.5, ^2 
1.687, Yj • -0.761 


0 . 02 , 


Note that 


H(8 


> i „ , o.m'V • V/iW.v-' ■ 


and 


n(H) 


,, ■ -.-.vikiv ■ . 

I* S I I - 


0-70; 


TIii-rctt»re, hUwv I iO, j>* 5871 ^ 


n 


2f,^ ' ui lor 
n 


0.2, the 


[0.20] 

[o.H4j 

frequency rertponue at the resoname of the high frequency mt>de versua 


•xctuator provldea 20 loglv*'v„'; ' -13 db magnitude attenuation of the 

^ U * 07 j 


the resonance of the low frequency mode In Example 5. Therefore Example fa 
falls Into c.itegory (vl). Note th.it r ■ e * 0.61 • e " 

^-5r,2wjT ^ ■ 0.7B. However, the pole-shifting desired is to a 

radius of 0.861. Therefore the closed-loop transient effect of the 
.nctu.itor on the low frequency mode Is expected to be appreciable but 
possibly tolerable. Despite the unity DC gain of H(«), its el feet on 
the pole-shifting will also alter the step response tracking accuracy. 

Ttiese effects are quantified in Table 1-1. Note that, for each actuator 
configuration, a.s imreases the controlled system exhibits Instability 
as one or mc>re of the closed-loop poles moves outside the unit circle; but the 
range of retaining stability Is Increased by each actuator configuration 
relative to (3-fal). (The factorization scheme used is the same as that in 
Appendix C. .‘lee Appendix E for det.iils.) 

Table 3-2 summarizes the six ex.impleH that will be used to meet the 
objectives of the preceding section. The next section presents the adaptive 
control algorithms to be applied. 
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Ex# 

plant 

pole** 

plant 

zero(a)* 

DC gain 

actuator 

pole* 

actuator 
DC gain 

— 

desired 

pole(a)* 

desired 
DC gain 

1 

0.95,0.2 

0.3 

1.31 

- 

- 

0.8 

2 

2 

0.9,0. 1 

0.011 

0.99 

- 

- 

0.65 

2 

3 

0.9, 0.8 

0,32 

0.9 

- 

- 

0.5 

2 

4 

0.9, 0.7 

0.4 

4.0 

- 

- 

0.6 

2 

5 

0.923*j0.231 
-0. 762 0.570 

-0.960 

-0.144±j0.942 

4.04 

- 

- 

0.844-^j0.171 

2.21 

6 

0.923-*j0.231 
-G. 762 ij 0.570 

-0.960 

-0.l44±j0.942 

4.04 

0.607 

1 

0.844±j0.171 

2.21 


* All ainguiarl t iea are z-plane equivalents 


TABLE 3-2: Examples Summary. 
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IV. ABAPTIVE COtrrROL ALGORITW^ 

Presently, adaptive control KcheneB can be categorlxed by the 
following characteristics: (i) basis for parameter update scheme via 

gradient (G) search procedures or stability (S) theory analysis, 

(ii) formulation of parameter estimate error as equation (E) or output 
(0) error, and (iii) designation of parameters to be Identified as plant 
parameters for indirect (I) adaptive control or as controller parameters 
for direct (D) adaptive control. Of the eight combinations only seven 
are presently available with GOD not yet developed. Of the remaining 
seven, two pairs are pairwise identical, GEI/SEI and GED/SED, duo to the 
resulting equivalence of gradient and stability based equation error 
parameter estimation. Tlie remaining five distinct classes providing 
adaptive solution to the model following problems of the pre -ding 
section are detailed below with appropriate source reference. The 
designer selected variables of each are listed following each algorithm 
statement . 

A reduced-order adaptive model-following problem requires specifi- 
cation of a time-varying feedback control law causing the output of the 
plant 

m _ n _ 

y(k) “ b, u(k-j) + [ a. y(k-l) (4-1) 

j-1 J 1*1 

to asymptotically "follow" that of the stable reference mudei 

n 

s(k) * >; r(k-l) + ^ i. s(k-l) (4-2) 

i-1 

without a priori knowledge of the plant parameters and or plant 
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dimensions m .ind n under the assumption that the plant is described 
by 

m n 

y(k) « I b. u(k-j) + I a. y(k-l), (4-3) 

j-l i-1 

where m >ra and n > n. Only if ra-a, n-n, b^-bj.and arc the fol- 

lowing aolutlons expected to cause y->s. If, as In examples 5 and 6 of the 
preceding section, the model of (4-2) Includes zeros, which are restricted 
to match those of (4-3), this converts the problem to one of adaptive pole- 
shifting. Since these zeros are identified in Indirect adaptive control 
the pole-shifting objective presents no apparent difficulty, though no 
globally convergent, even nonreduced-order , indirect adaptive pole- 
shifting scheme has been proven due to sufficient excitation requirements 
130) or matrix polynomial solution singularity problems [31]. Direct 
adaptive controllers must assume the availability of the plant numerator 
coefficients, i.e. the plant zeros, [32] or the stability of the plant 
[33] to achieve globally stable pole-shifting (or replacement). Therefore 
such direct adaptive control schemes (or those requiring numerator 
cancellation [17]) are not presented for examples 5 and 6. These restric- 
tions on direct adaptive control support the statement that presently 
Indirect adaptive control appears more suited to the ROAC problem 
requiring pole-shifting. 

The algorithms, specialized for the examples of the preceding 
section and written in implementation sequence, are: 

For Examples 1-4 : 
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GE.l/SEI 117), 1211, U>1*J. 13M - (36J: 


e(k-n - y(k-n - a(k-I)y(k-2) - b(k-l)u(k-2) 
r 


.i(k) - a(k-l) + 


h+uy* (k-2)+i'u‘ (k-2) 


y(k-2)e(k-l) 


b(k) - b(k-l) + 




[_ h+i.y"(k-2)Hvu"(k-2) J 


u(k-2)p(k-i) 


j.(k-n - 

b(k) 

f(k~l) - 

b(k) 

IVh innor-Ht?lei'ti'd variable requlreroentB: 

1 . ' 0 

' 0 

0 V ii . 2 

h Is ni'mlnally fhosen as one unless auoh a value would oause 
Then h Is chosen within the ran^e in (4-11) such that |b(k)| 

(JOI [231. [2b], [151 ; 

y(k-l) » a(k-l)y(k-2) + b(k-l)u(k-2) 

V(k-l) « y(k-2) + a(k-l)\(k-2) 

'i (k-l) - u(k-2) + a(k-l)Y(k“2) 
e(k-l) - y(k-l) - y(k-l) 


a(k) » a(k-l) + 


h+u\^(k-l)+,’Y^(k-l) 


^(k-l)e(k-n 


b(k) - b(k-l) + 


h+u\"(k-l)-h'Y (k-l)J 


Y(k-l)e(k-l) 


(4-4) 

(4-5) 

(s-6) 

(4-7) 

(4-8) 

(4-9) 

(4-10) 

(4-11) 

b(k) * 0 


(4-12) 

(4-13) 

(4-14) 

(4-15) 

(4-15) 

(4-17) 


IV. 4 


g(k-l) ■ 

6(k) 

(4-18) 

f(k-l) - 

b(k) 

(4-19) 


Designer-selected variable restrictions: 

Ij > 0 and small (4-20) 

p • 0 and small (4-21) 

0 «. h < 2 (4-22) 

h Is nominal ly chosen as one unless such a value would cause b(k) = 0. 

Tlien h Is chosen within the range In (4-22) such that |b(k)| > c ' 0^. 


SOI: [8], [24], [35]; 

y(k-l) “ a(k-l)z(k-2) + b(k-l)u(k-2) (4-23) 

V(k-l) - y(k-l) - y(k-l). +_cL[y.(k-2) - z(k-2) J (4-24) 

h + uz^(k-2) + pu^(k-2) 

a(k) “ a(k-l) + wz(k-2)v(k-l) (4-25) 

b(k) » b(k-l) + pu(k-2)v(k-l) (4-26) 

z(k-l) * a(k)z(k-2) + b(k)u(k-2) (4-27) 

g(k-l) « (4-28) 

b(k) 

f(k-l) - — (4-29) 
b(k) 

Designer-selected variable restrictions; 

> 0 (4-30) 

P ^ 0 (4-31) 
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Ra 


-1 


‘■•.•■‘I 


> 0 


V z ■ 1 


(4-32) 


q can be chosen as zero for any stable a^ and (4-32) will be satisfied. 

1 <, h < 2 (4-33) 

h Is nominally chosen as one unless such a value would cause b(k) ~ 0. 


Then h is chosen within the range in (4-33) such that |b(k)| > c * 0^.. 
GEb/SED tl7J[37]; 

v(k-l) » cr(k-2) + d y(k-2) - y(k-l) (4-34) 

g(k-l) - g(k-2) + (4_35) 

h[l+iiy (k-2)+pu (k-2)] 

f(k-i) - f(k-2) + 

h[l+uy^(k-2)-fpu^(k-2)l (4-36) 

Designer-selected variable restrictions 

p > 0 (4-37) 

p > 0 (4-38) 

|hl •> and sgn(h * sgn(b) (4-39) 

SOD [17]-[20]: 

b(k-l) - jr^(k-2) + yy^(k-2) v4-40) 

Y(k-l) - (U+q)hS(k-2)v(k-2) + dy(k-2) (4-41) 

v(k-l) - (14hS(k-l)rHs(k-l)-y(k-l)+qts(k-2)-y(k-2)]-Y(k-l)} 

(4-42) 

g(k-l) « g(k-2) + pr(k-2)v(k-l) (4-43) 

?(k-l) » ?(k-2) + ;iy(k-2)v(k-l) (4-44) 
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D«slgn«r-selected variable restrictions: 

M > 0 
p ^ 0 

I 1) I 

|h| ■» 1 — and sgn(h) * sgn(b) 


Re 


1-fz"^ 


> 0 


V|z( - 1 


(4-45) 

(4-46) 

(4-47) 


(4-48) 


This last condition is easily satisfied by q ■ -d, which Khould equate SOD 
to OED/SED [18]. 

Each of the precedinp, schemts (4-4) - (4-48) provides g(k-l) and f(k-l) 
for 

u(k-l) - g(k-l)r(k-l) + f(k-l)y(k-l) (4-49) 


for application to (3-2) and generation of y(k). Then these recursions 
arc repeated. 

For Examples 5 and 6: 

GEI/SEI: 

e(k-l) - y(k-l) - tj^(k-l )v(k-2) - a^(k-2)y(k-3) 

- 8j(k-l)u(k-2) - e^(k-2)u(k-3) (4-50) 

v(k-l) - e(k-l)/[h+M^y^(k-2) + W2y^(k-3) + p^u^(k-2) + p^u^(k-3)] 

(4-51) 

a^(k) « aj^(k-l) + p^v(k-l)y(k-2) (4-52) 


02(k) - a2(k-l) + P2v(k-l)y(k-3) (4-53) 

63^ (k) » 83^(k-l) + Pj^v(k-l)u(k- 2 ) (4-54) 

B2(k) “ B2(V.-1) + 


(4-55) 
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DftBlgner-selected variable reNtrictlons ; 


0 

0 


Pj^ > 0 

p2 0 
0 < h < 2 


(4-56) 

(4-57) 

(4-58) 

(4-59) 

(4-60) 


h is nominally chosen as one unless such a value would cause B 2 (k) = 0 or 
3j(k)aj^(k) + 62 ^*^) ~ (k)/P 2 (k) ' 0. Then h is chosen within the 

range in (4-60) such that the absolute value of both terms ■» g = 0^. 


GOI ; 

y(k-l) - a^(k-l)y(k-l) + a 2 (k-l)y (k- 2 ) + 8 j^(k-l)u(k- 2 ) + 62 (k-l)u(k- 3 ) 


(4-61) 

X^(k-l) - y(k-2) + a^(k-l)A^(k-2) + a 2 (k-l) A^(k-3) (4-62) 

X 2 (k- 1 ) - y(k-3) + «j^(k-l) V 2 (k- 2 ) + .i 2 (k-l)X 2 (k- 3 ) (4-63) 

jj^(k-l) “ u(k-2) + aj^(k-l)ir j(k-2) + a 2 (k-l)Y^(k- 3 ) (4-64) 

Y 2 (k- 1 ) - u(k-3) + a^(k-l)Y 2 (k- 2 ) + i 2 (k-l)Y 2 (k- 3 ) (4-65) 

v(k-l) - [y(k-D - y(k-l)]/[h + (k- 1 ) + ^ 2 ^ 2 ^ (k-D 

+ (k-1) + ^ 2 ^ 2 ^ (k-1)] (4-66) 

a^(k) * a, (k-1) + u^vCk-l) X^(k-l) (4-67) 

« 2 (k) =* a 2 (k-l) + U2v(k-1) ^ 2 (k-l) (4-68) 

hj(k) - P^(k-l) + p^v(k-l)Yi(k-l) (4-69) 

P 2 (k) - 82 ^*^"^^ + P 2 v(k-l)Y 2 (k-l) (4-70) 
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Oesignar-aftlccted variable restrlctlona: 


Pj > 0 

and 

small 

(4-71) 

1^2 0 

and 

small 

(4-72) 

^ 0 

and 

small 

(4-73) 

p., 0 

and 

small 

>.4-74) 

0 < h < 

2 


(4-75) 


h is nominally choaen as one unless such a value would cause - 0 or 

hj(k)aj^(k) + 1^2(1^) “ chosen within the 

^ 4. 

range of (4-75) such that the absolute value of both terras ^ t ' 0 . 

SOI: 

y(k-l) - u,(k-l)z(k-2) + «,(k-l)z(k-3) + S, (k-l)u(k-2) + 8,(k-l)u(k-3) 

(4-76) 

(4-77) 

(4-78) 
(4-79) 
(4-80> 
(4-81) 
(4-82) 

(4-83) 
(4-84) 
(4-85) 
(4-86) 


{y(k-l)-y(k-l)+q. [y(k-2)-z(k-2) ]+q, [y(k-3)-z(k-3) 1 } 

v(k-I) ■ = 5 

lh+iijZ^(k-2)+u2z"(k-3)-h5ju'^(k-2)+P2U (k-3) } 

a,(k) “ aj(k-l) + |.j^v(k-l)z(k-2) 

02 (k) ■ a^(k-l) + ii.,v(k-l)z(k-3) 

Bj^(k) « ti^(k-l) + p jV(k-l)u(k-2) 
p 2 (k) ■ ^^(k-l) + p,v(k-l)u(k-3) 

z(k-l) - a^(k)z(k-2) + ,^(k)z(k-3) +hj(k)u(k-2) +p2(k)u(k-3) 
Designet-selected variable restrictions; 

112 0 

pj > 0 

p2 > 0 
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Re < 


1 ^ ^ -2 

1 + qj^z + q^z 

1 - 


> 0 


V|zl * 1 


(/I-87) 


qj could bf chosen as -0,98 .ind q., as zero for the rapidly sampled 
exanqiles 5 and h and (A-87) would be satisfied. 


(A-88) 


h Is nominally chosen as one unless such a value would cause ' 0 or 

(k)/.',(k) - 0. Tlien h is chosen within the 
range of (/t-88) such that the absolute value of both terms • 7 • 0^. 

Each ot the last three algorithsts in (A-50)-lA-88) provides a^(k), 
a.,(k), !' 2 tk) for use In (reler to (3-27)-( 3-301 ) , 


(k-1 ) ^ (kl 


(4-89) 


vj(k-l) 


[r,j(k-na2(k) + (nj(k-l)ajtk)-a 2 (k)+y 2 H<*j(k)-a 2 (k) 8 ^(k)/ 82 (k))l 

lhj(k)a^(k)+ 82 (k)-a 2 ^*^^^l" 


(4-90) 


(4-91) 


n,(k-l) - (k)Hj (k-11 "*■ >2 “ hj(k)vj(k-l) - -i,(k) 

C2(k-n - T^(k-l).i2(k)/p2(k) (4-92) 

for jMrvmeterIzat ion I'f 

u(k-l) « *rlk-l) + 'i_^(k-l)u(k-21 + n2(k-l)u(k-3) 

+ v^(k-l)y(k-2) + 72(k-l)y(k-3) (4-9 1) 

lor application to the difference equation representation oi (3-17) 
y(k) ■ ^ ‘21 ^ Oj^^)y(k-l) + (»22 “ ‘‘ll*’’! Uplylk-Z) 

— a ..^^2 1 ) y (k— 3) — ^22*^22^^^""*^^ 

+ (Bj^^ + “ ^’‘'n'^21 " *^12 ^ ^2l'^l] " 

(Bjin22 ■*■ ‘' 12*21 ^2l‘'l2 ^ 22 ^^ i)“(k-3) 

— + f 2 T'*|,^) u (k— 4 ) 


(4-94) 
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gen«r«Cing y(k). Itot* that all of tha pracadlng «cha«wt in (4-i^)-(A“9 3) 
ua« the current Inforwatlon as efficiently as possible and therefore 
result in causal rather than strictly causal adaptive controllers [17]. 
(See Appendix F for programs of each of these controllers.) 
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V. TEST FORMATS 


HavlnK speclflsd the nnn~adaptlve and adaptive control strateglea 
and objectives the remaining Itame to be selected for the testing of the 
guidelines In the second section are the reference inputs, adaption step- 
size constants, parameter estimate initializations, and performance 
measures. 

Selection of the reference signal is important fur two reasons: 
richnesH and magnitude. The magnitude of the reference signal is important 
diw to the nonlinearity of the adaptive control system and therefore the 
anticipated transient (at least) differences due to different signal 
levels. The richness is important because it is well-recognized {38} 
that a sufficient "input" richness is necessary to perturb all of the 
modes of the system for their "identification." This "input" differs for 
the Indirect .and direct approaches .and the equation and output error 
formulations. For the equation error Indirect approach the plant input 
and output must be suftlciently rich. Since the plant is linear this 
tr.inslates to a sufficient frequency content of the plant input (or control) 
signal, of which the reference is only one component. The difficulty of 
assuming adaptive control input richness based on reference signal rich- 
ness is addressed in [30]. For the output error indirect approach the 
plant input and Identifier output roust he sufficiently rich. In an open- 
loop output error identification task the Identifier output richness can 
be translated [39] solely to input richness necessity under the assumption 
that the identifier output asymptotically converges to the plant output, 
which requires sufficient identifier order absent from ROAC. For direct 
ad.iptive control the control parameter identification problem must be 
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recast as an open-loop identification task as in [19] and (37). In th« 
equation error case the richness rcqulre»ents are converted principally 
to conditions on the a>odel output to be tracked or, if generated by a 
linear, time- invariant model, to the forcing or reference signal [39]. 

Once again this assumes sufficient controller order. The requirements 
for output error based direct adaptive control, though as yet unspecified, 
are assumed to be similar. 

All of the preceding richness requirements are for complete identifi- 
cation and assume non-reduced adaptive model order. As is well-recognized 
[22] [34] [37], for a control objective coaqslete identification may not be 
nect*Hs.iry. Consider the control of any order system with a single constant 
output feedback gain and a reference signal gain. To achieve aHymptotlc 
convergence to a particular I>C level set-point any feedback gain stabiliz- 
ing the system In conjunction with the appropriate reference gain will be 
adequate. This applies readily to examples 1-4 in section three. As 
noted in [40] fur open-loop output error (and by implication equation 
erroi ) tdent If icvit Ion the richness of the input signal lor stable reduced- 
order identification is dependent on the reduced-order model dimension, 
not that of the higher order plant. From a frequency response point of 
view, the viewpoint promoted in [41] for reduced-order modelling, the 
minimum number of sinusoids providing this richness matches the number of 
points at which the controller can specify the steady-state frequency 
response. Therefore, If this reduced-order controller stabilizes the 
system, asymptotically convergent steady-state control seems possible up 
to the richness level of the reduced-order model. Wlien the signal to be 
cracked (or equivalently the reference signal) has a higher frequency content, 
reduced-order tracking will onlv be approximate at best. The reference 
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aignal should be chosen to test this boundary condition, with the 
principal case of interest being reference signal over-richness. 

'fherefore the folltiwlng reference signals (all xero prior to k • 0) 
will be considered; 

INPUT 1 (unit step input): 

rj(k) - 1 (S-1) 

INPUT 2 (single low-frequency sinusoid): 

r^(k) - 2 ••ln(ku) 2 T) (S-2) 

Note that the sample period T is assusied to be one second for examples 1-4. 
Tlu* frequency I** chosen as 20% below the desired cutoff frequency for 
ex.uiiplcM 1-4, I. e. 0.178'il, 0.1A4f>), 0,S'>A“>2, and 0.A0866 radians per second, 
respectively, and 20% bel«iw the dominant natural frequency for examples 
b and 6, i.e. O.AO radians per second, with T * 0,S seconds, 

INPUT 3 (single high-frequency sinusoid); 

r^(k) » 2 sin(kwjl) (5-3) 

Again the sample period Is assumed to be one second for examples 1-A. The 
f requeni'V Is double the di-sired cutoff for examples 1-4, I.e. respective- 
ly 0,4Ab29, 0.8f)157, 1.38629, .ind 1.02170, and double the dominant natural 
frequency for examples 5 and 6, i.e. 1.00, with T ■ 0,5 secondn, 

INPUT 4 (non-zero mean, uniform, white noise): 

r^(k) r l-0.5,1.5J (5-4) 

This is clearly the over-rich input of most interest. 

Since the rate of convergence and parameter estltnate variance near 
convergence are conaaonly acknowledged as being affected by the adaptive 
step-size (SS) coefficients and pj , two cases will be considered. 
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SSj^ « 1; 

- pj - 1 (5-5’. 

SS2 - 0.1: 

Uj ■ Pj * 0.1. (5-6) 

Given the approxlinate nature of the gradient formulas for GOl embodied in the 
required p and p Bmallnoss noted in (4-20)- (4-21) and (4-71 )-(4-74) the p and p 
(5—5) may lead to unstable behavior for (.01 while (5—6) should prove more 

satisfactory. (Note that h in the normalizing terra of all of the algorithms 
is selvvcted as unity). 

For output error type algorithms such as SOI and SOD the error 
smoothing coefficients q in (4-24), (4-4?), and (4-77) are to be selected 
in satisfaction of (4-32), (4-48), and (4-b7), respectively. As noted in 
[42] selection of these q to make (4-32), (4-48), and (4-87) equal unity 
does not offer the most rapid convergence. Furthermore, in a reduced- 
order application these q can influence the mean convergence point. Such 
factors may be used to advantage in (4-48) where f is known a priori, but 
in (4-32) and (4-87) the strictly positive real (SPR) condition is dependent 
on unknown plant parameters. Therefore several values will be tested. They 
are summarized in Table 5-1. 

Since only local convergence is currently anticipated three .settings 
of the initial plant (and therefore controller) parameter estimates (EST) 
will be considered corresponding to the nominal values from control design 
based on neglecting the second mode (i.e. e * 0 for examples 1-4 and * 0 
for examples 5 and 6) and ± 20% of these values. These values are summar- 
ized in Table 5-2. The controller parameter Initializations for the direct 
schemes are based on the corresponding plant parameter initializations and 


Example ( m) 

Algorithm 

Error Smoothing 
Coef f lcient(s) 

LaIt© 1 



-0.9 

sc- 


SOI 

-0.95 

SCO 

1 


-0.97 

SC4- 


-0.7 

sc- 


SOD 

-0.8 

SCO 



-0.9 

SC+ 


SOI 

-0.9 

SCO 

2 

Si*D 

-0.85 

SCO 

-- 





sol 

sou 

-0.9 

-0.5 


— > ■ 

- 

* ‘ 




-0.8 

sc- 


SOI 

-0.9 

SCO 



-0.95 

sc+ 

4 


-0.4 

sc- 


SOD 

-0.6 

SCO 



-0.8 

sc+ 

5&6 

SOI 

qi— 0.9,q2"0 
qj«-1.846,q2*0.9515 




qj«— 0.96,q2*0 



TARl.K 5-1; Error Smoothing Coefficients 
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(3-8) and (3-9) for exaoplea 1-4 and (3-27)-(3-30) for examples 5 and 6. 
Note the instability of the EST ■ + 20% initializations. 

Since the control objective of all of the examples in section 3 is 
model-followinR, the tracking error is the principal performance measure. 
For the expectation of asymptotic tracking, the figures of merit should 
be divisible into short and long term quantities. The mean and variance 
of the tracking errors should be normalized and tabulated for comparison 
among the examples and between the fixed and adaptive contro’ schemes. 
Similarly tlie Input cost and controller (and plant, for indirect schemes) 
parameter estimates should be observed. Therefore the following quantities 
will be tabulated at k^ ■ 1 , 2, 5, 10, 20, 50, 100, 200, 500, .md 1000 
iterat ions. 

Instant.meous desired output: 

IDO(k) * s(k) (5-7) 

Instantaneous controlled output: 

I0(k) - y(k) (5-8) 

Normalized instantaneous tracking error: 


DITEtk) . 


Normalized average segmented tracking error 

k. 

f 1 ) 

NASTE 




I NITE(j) 

j-ki_i+l 


(5-9) 


(5-10) 


Normalized segmented tracking error variance 


NSTEV 


[n(n-l) J 


k.^ 

n I NITE^(j) 


r ^1 

- r ^ 


12 ') 


NITE(j) 


j-ki_i+l 


(5-11) 
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(Note th.it n Ih «inu8 the nuiaber of tine* «(k) < 10~^, which 

are excluded due to (^-9). If NSTEV is negative due to round-off, it 
la printed .ih ziTo.) 

Seginenti-d aver.jge squared Input 


k -V- ): u“(i) 

1 ‘'i-l + l 


Segmented average controller parameter estimate for f 


r 1 I I 
sAi I'K F - 1, *: V "- I 

^*^1 *^1-1 ^ j-kj j+1 


(Similarly for g, a, h. a, t , v, and n.) 


Segmented controller parameter estimate variance for g 


SOPFV C. 


(kf - k,^j)(kj - 1) 






-I I h])]' 




(Similarly loi g, a, b, «, S, v, and n.) 

These quantities are tabulated lor each of the combinations of example, 
adaptive control algorithm, reference Input, step-size weights, error 
smoothing coefficients (If necessary), and par.ameter estinwite Initiali- 
zation. (See Ai'jiendix K lor listings ot the programs compiling these 
figures of merit, 1 Table '>-3 summarizes all the various combinations 
that were simulated. 






T— 

Control Alg. 

Inputs 

nonadapt ive 

all 4 

CEl/SFl (1) 

all 4 

GOI (2) 

all 4 

SOI (3) 

all 4 

(JED/SED (A) 

all 4 

SOD (5) 

all 4 

nonadapt ive 

all 4 

CEI/SKI (1) 

all 4 

COl (2) 

all 4 

SOI (3) 

all 4 

GED/SED «i) 

all 4 

SOU (5) 

all 4 

nonadapt ive 

all 4 

GEI/SEl (1) 

all 4 

GOl (2) 

all 4 

SOI (3) 

all 4 

vJKD/SKD (■'.) 

all 4 

SOD (b) 

^all 4 

nonadapt ive 

falTT 

(d l/SEI (1) 

all 4 

U'l U) 

all 4 

SOI (3) 

all 4 

GED/SED (.4) 

all 4 

SOD (5) 

all 4 

nonadapt ive 

all 4 

GKI/SEI (1) 

all 4 

GOl (2) 

all 4 

SOI (3) 

all 4 

nonadapt ive 

ail 4 

GEI/SEl (1) 

all 4 

GOl (2) 

all 4 

SOI (3) 

all 4 


I’ar. K«t 
Inlt. 


KST - 0 
dll 3 

all 1 

all 1 

all 3 

all 3 


EST - 0 
EST - 0 

EST - 0 

EST - 0 

EST - 0 

EST » 0 



all 3 
all 3 


TABLE *>-3: Tested Combinations. 
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VI . KK SUl.TS 

IXje to the number (1008) of printoutM for the ex.vwpleH «hovm In 
Table 5-i the tabulated reaulta are under separate cover. (See Appendix F 
for pertinent prop,ram listings.) The next section provides an evaluation 
of these numerical results. 
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VII, TEST RESULTS INTERPRETATIONS 

Refer to eectlon 2 for the etateieent of the objertlvee of this 
study. Tlie four questions raised in that section will be addressed 
in this section with reai>«et to the slnulations of Section 6. 
liUE.STlON (i): 

Are heavily d.imjjeil open-loop modes negleot.ible rel.it ive to more 
lightly damped modes? In particular: 

a) Are the accepted criteria for neglecting mt»deH, based on 
relative damping, realistic for the examples cimsldered here 
in a llxed ( non- adapt Ive) reduced-order controller? 

b) For these examples, how do the various adaptive control 

schemes lare under the reduced order modeling "rules" 
concerning relative damping? More specific questions might 
be: Do the reduced-order adaptive control systems remain 

stable? Can the adaptive roeehanism compensate for mls- 
roodeling to yield Improved performance over non-adaptive 
(modal) controllers of the same order? 

In example 1, the open loop and closed loop modes are aepar.ihle 
based on the accepted criteria for relative damping, so the good closed- 
loop (non-adaptive) reduced-order performance was to be expected 
(despite the fact that In this example and all the others the reduced-order 
modal model was based strictly t»n the partial fraction expansion term 
without DC gain correction). F.xnmple 2 represents a marginal Iv separable 
system in closed loop, yet the (non-adaptive) reduced-order controller 
t racking errors are only slightly worse than in example 1. This indicates 
that, for these examples, this boundary between separable and Inseparable 
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modts is reasonable hut fuzzy (as expected). 

Ex.tnpltf 1 poMeesees neither clo*t»d loop nor open l«»op separability, 
but a pole- zero ne.ir-i .incel lat ion should lc.tve one »ode doisinatinK the 
system bahavior. As shtnzn in the slsmlat ion, t^’e reduced-order 
Buidelinjt causes significant errors. This is presuswbly due to the relative 
closeness ol the "cancelled" pole to the unit circle cosip.ared to the 
"dofflin.int" pole. Siaill errors in cancellation in this case may have 
large effects isi system behavior and the perfonsance of reduced-order 
t ont ri>l . 

The approximation of a pair of complex poles by a simple real 
axis pole is studied in example U. Here, no first order model can be found 
to cliiselv represent this second order system, and considerable degrad- 
ation In perliMmaace by reduced order control should be expected. The 
slmulat imis .show this to be true; tin* steady state tracking errors are 
sl gnl f leant . 

ot the five adaptive c.mtrol algorithms, all but the second one, 
the gradient -output -erri'f- indl rect (Uol) scheme, were staMi*. In I'xample one, 
which has the no.st accurate reduced order nKidel, a larger number of test 
combi nal ii'ns t>f inputs, step sizes, and initial parameter estimates 
produce stable responses than the other examples. Even so, the only 

consistent results were for zero Initial parameter error (witli re.spect to the 
extracted modal model) and the step and random Inputs. In these cases, the 

step re.sponse is improved over the non-adaptlve case to essentially zero 
steady state error, while the random input causes a stable yet more widely 
ranging response. The other four adaptive algorithms produced essentially 
equivalent responses for roost Inputs, step sizes and smoothing coefficients. 
Parameter convergence occurred between 50 and 500 iterations for 
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noHt runn, the random input caae being oonaiitcntly faster. Also, tUe 
error between the desired and actual outputs always resulted in Iie- 
provement over the non-adaptive output error. Some minor differences 
in responses were: 

• Example 3 usually responded f....ter, due to the dominant pole 
being more highly damped. 

• Algorithms (stobixity-output -t«. lor-direct, SOD) prod red 
consistently smaller steady sta‘e output tracking errors. 

• Smaller step sizes improved the ateady state tracking error 
variance, but lengthened convergence times. 

• The smootiiing coefficient .SO In algorithm 3 produced more 
accurate steady state tracking than the other smoothing 
Coe f f Ic lent 8. 

• The adaptive parameter estimates have higher variances tor 
example 2 than example 1 indicating more active seart hing due 
to the marginal dominance of the desired reduced-order 
objective in the problem setup. 

Except for the GOI algorithm (^ 2 ), all adaptive schemes provided 
Improved output tracking of desired models over non-adaptive control 
in examples 1-4. This is the result of the adaptive algorithms in- 
clination to select some controller parameterization to "better lit" 
some control objective. In the indirect algorithms, the Identified 
plant often has parameters that do not match either those of the actual 
plant or those of the assumed reduced order (modal) model. These effects 
are explored further In the discussion of the third question, on choosing 
the "best" reduced order model. 
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.1 

1 

QUESTION (ii)j ; 

j 

i 

If actuator dynaalcs are Im ludfd In thf iiyatem aurh that the ^ 

residual nodes have frequencies within the stop band of the actuator I 

i 

frequency response, does isqiroved reduced-order control result? 

examples 5 and 6 were used to test the effects of actuator in- 
clusion. Example b included no actuator dynamics and wu« seen to be 
unstable for a large number of runs, particularly those with initial 
parameter estlnwte errors relative to the dnninant mode. Examplt* h 
included the three actuator configurations shown in Figure 3-12, where 

actuati’r configuration 2 in the siMulations corresponded to a< lu.Jtor . 

(n) In the figure, i to d>), and 4 to (i‘). 

In the non- adaptive simulations, inclusion of an actuator provided 
Improved tracking fi>r those runs of example b that were stable based cm the 
reduced- order moual model. For the fixed (initial) reduced order mticlel 
parameter estlra.ites that resulted in an unstable control system, actuator 
inclusion generally did not provide stability. Actuator conf lgur.it Ions I 

2 and 4 provided up to an order oi magnitude improvement in tracking 
compared to 3 and 1 (no actuator), which were similar in response. 

Under adaption, the moat improvement was obtained, conversely, 
by including an actuator as In configuration 3. This arrangement 

produced stable simulations for all but a few runs, vrfiere the output I 

was apparently unstable but still bounded within the 1000 iterations 

tested. For the stable adaptive usage of actuator configuration 3, 

nearly all cases provided between 1 and 3 orders of magnitude improvement 

in tracking error over the other three actuator conf Igurat Ions. 

Actuator arr.ingements 2 and h behaved almost identically during 
adaption under all condititms. They were seen to improve the response 


I 
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ov«r that without actuator Uynanlcs tor auat runs, but ocfoaional ly their 
Inclusion wuraened the output tracking’ The only conaistent lapruveaent 
over configuration J was for Input J, the high I roquency einueold, where 
an order of aagnltude difference was aoaetlaea seen. 

Among the three algorithms used on examples 5 and 6, the first 
(GED/SED) and third (GKI ) behaved essentially the same, giving the 
responses discussed above. Algorithm 2 (001) caused unstable responses 
for all combinations of inputs, step sixes. Initial parameter estlputes, 
and actuator conf Igurations. 

For these examples. Inclusion of a finite bandwidth actuator 
generally Improved the perforowtu'c of the reduced-order control. However, 
tliere were some cemhlnat (ous of inputs. Initial parameter estisuites, and 
adaptive step sizes that c.iused poorer behavior with an actuator than 
without, live ictu.itor at l angeraent that allowed both tlve plant Identifier 
and controller to sample the actuator output was clearly the best choice 
for the majority i»f cases liere. This structure effectively includes 
the actuator .is part of the c»mtroller. As noted In section i, this ia 
not the most physically teaslble configuration since the actuator outputs 
can be Inaccessible. 

QUESTION (111): 

Does an "optimal" reduced order controller result from retaining 
the dominant modes of the plant for the reduced-order plant model, or 
by choosing a model (or controller) parameterization that does not cor- 
respond to plant modes but is based on matching (in some sense) a reduced 
order control objective? 

As phrased, the latter alternative Intuitively seems "better". 

This question can be answered for these examples by comparing the "zero 
Initial parameter estimate error" parameters with the steady state parameters 
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In the adaptive runa. If the difference In paraiaetera la allKht, then the 
aelectlon of the dominant plant mode lor the control model would seem to be a 
rtaaonnble choice with reapect to the adaptive parameter estimator update 
criterion. If tli»‘ parameter difference la conalderable and the .idaptlve rcnponhea 
are deemed belter than the non-.'idapt 1 ve onea • then a reduced order model 
based on those parametera selei ted by the "behavior approxlmat inn" adaptive 
algorithms would seem the better cliolcc. 

As woulu be expected, the ad.iptlve runs shi*wlng the greatest improvement 
over the fixed controller runs were those with the largest differences In init- 
ial and steady st.Ue par.ipieters. The implication is that choosing dominant 
nu^des out of the set of plant modes for the reduced order model in often not as 
"optim.!!" a.*? cljoosing model "modes" based on some performance criterion. 

QUEST! ON (iv); 

Are indirect .md direct adaptive controllers equivalent and 
liitercluingeable In reduced-order appl Icat 

The absurdity of such an expkectatlon is underscored by a coaparlson 
of the algorithms In section IV. Note that indirect schemes involve 
division by estlnuJted quantities as In, e.g. , (4-7) and (4-fi), but that 
direct schemes do not. Clearlv, if the estimated quantities ever 
approaih zero the resulting <iuotIent Is quite large. A similar mechanism 
does not seem to appear in the direct algorithms. The claims of 
equivalence seem to rest on the (linear?) transformabi 1 ity between the 
indirect and direct approaches in simple cases, such as inverse 
or model-following control. If the algorithms are not Identical then, 
as clearly deim'nst rated by this study, the tnansie tt differences in the 
full order case translate into steady-state differences In the reduced- 
order case. Unfortunately neither indiiect or direct adaptive control seems 
consistently superior to the other in all situations. 
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VII 1. CONCI.USION 

The bulk i>f the data generated fro« the 1008 alwulatlona .itteata 
to the difficulty of drawing definitive cuncluaiona fro« auch aimulation 
Btudiea. This la further coaiplicated for atudiea of adaptive contrullera 
by the number of varlablea, e.g. plant parameters, control objective, 
initial parameter eHtiauttes, atep-alae weights, error smoothing coef- 
ficients and input richness, effecting overall performance. Conclusively 
examining all possibilities is clearly unfeasible. Even n foc„rfcd study, 
such as this one with respect to the four objectives conatented on in 
the preiedlng section, can only be expected to provide heavily qualified 
remarks. At best, such a studv can t>nly be used to test .m>l .lugment 
the "folklore". As particisl.ir effects are Isolated without contradict l<*n 
further simulation studies are In order to test seeming generalizations, 
ihe ones tnat survive closer scrui inv pri'Vlde guldeposts tor restricted 
an.ilvtl.al studies, that c.>uld pr«>vc equal Iv fruitless If not initially 
We 1 1- locused. 

The data generated by this study nuiy provide a dat.ibase for 
judging eventual conclusions of ftiture studies. A different format, 
e.g, graphical, for even the tests nin here might facilitate evaluation 
of different conclusions. The one general (tautological' conchmion that 
seems well documented by this report is that all present adaptive control 
schemes will fail in some reduced-order application. Though most of the 
(stable) adaptive controllers tested seem to better the performance of 
fixed modal controllers this does not present a conclusive argument for 
unsjipervlsed adaptive control. It is trite but true that further study 


Is necesrarv. 
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Appendix A: Root Locuk Program Listings and Tabular Outputs 


Contents: 


. A-1: Program determining p^ and p^ in (3-12) for Examples 1 

A-2; Pj^, p^t h, and f for Example 1 (Pig. 3-1) 

A-3: pj^, P 2 » h, and f for Example 2 (Fig. 3-3) 

A-4: Pj , P 2 » h, and f for Example 3 (Fig. 3-5) 

A-5: Pj^, P 2 » h, and f for Exmnple 4 (Fig. 3-8) 

A-6: Program determining and q 2 (3-15) 

for Examples 1-4 

A-7: qj^, q 2 » h, and f for Example 1 (Fig. 3-2) 

A-8: q^, q 2 » h, and f for Example 2 (Fig. 3-4) 

A-9; q^^, q 2 » h, and f for Example 3 (Fig. 3-6) 

A-10: q^, q 2 » h, and f for Example 4 (Fig. 3-9) 

A-11; Program determining “ a^, q 2 ^f, 

i.e.(0 fj In (3-14), for Example 3 
A-12: q^,, q 2 » h, and f for Example 3 (Fig. 3-7) 

A-13: Sample Plot Program (Fig. 3-1) 
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13 
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100 
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fOU 

300 
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WAIF lV»'^>*-^v,PACE*100, IIMt'lC'O 

RTLl: 

ROOTlOCUS INTcKPWETATIUN OF THu UNADAPTtt PLANT. 

T.F.:V(2 )/M 2)>^( (B^E )»Z-(l3»A2^t*Al ) )/( Z**r-(Al4Ar-»f n'»Z-»/ l*A 

24F*(B*A1+L»A1 ) 

DIMENSION t (A) , XX I (500) ,XX2( bOO) . YY 1 ( b 00 ) »Y YK 1 0( ) , A Z ( boO ) , F2 ( bO ) 
DIMLNSIUN Al(A) ,A2(A),bU) |YM1( bOo ) ,YH1 ( i OO » 
complex &C,FC,ZC(2» 

DATA E/O.Olt-O.Ul.0.O2f-0.3/ 

DATA Al/C.S>bt0.;»0.Vt0.9/ 

DATA AlVO.^.O.l.O.Bf 0.7/ 

DATA B/0. 065, 0.1,0. 01,0.5/ 

DO 300 J*1,A 

WRITE! 6,200 )A1 (J) ,A1 (J) ,B( J) 

FORMAT (IX, »A1»* 13.6, 5X,'A2**, FI 5.6,5 X,*B»* ,F 13.6) 

WKIiC(o,AOO)E( J) 

FORMAT!//, • t-»,Fl3.6) 

DC 100 Ksl,bO 
IF (J. EC. A) bd TU 13 
F?(K)*-0.2*(K-1 ) 

c -0 ro 1 

F?(X)sO.?*(K-l) 

G«(Al( J)4Al(j)4Fl(K)*(b( J)4L (J) ) ) 
F=^(AI(J)*A2(J)4F2(K)*(B(J)*A2Tj)+E(J)«A1(J) ) ) 


OCs^CMPLX(b,0.0) 

FC=CMPLX(F , J.O) 

ZC (1)*TgC4LSLRT (GC*+1-a.G*FC ) )/2 .0 
ZC (?) = (GC-CSLHT (GC’»4?-A.0>»FC ) )/?.0 
2R1=TEAl(ZG( 1) ) 

ZP7 = R6AL(ZC(?) ) 

ZI1=AIM/G(ZC(1)) 

ZK=A1MAG(ZC(4. ) ) 

XXKK) -2R1 
XX? =ZR? 

YYKO =Z 12 
YYKO-Zll 



YM1(0=YMAG1 
YM?(K) *YMAb2 
V=P(J)4E (J) 

W* (B( J )*A? ( J)4l ( J)*A1(J) ) 

AZ CO^W/V 
fONTlNUL 
WR1TE(6, ot> ) 

FLRMAT (/ ,5X, 'REAL* ,1 X, • 1 MAbl N • , 8 X , • MAG 1 * ,FX, 'lUAL • , « X , • I MAG IN • , 3X 
•MAG2* ,t-X, 'ZIK')' ,8X, 'FELDBACK* ) 

DO 500 lsl»50 

WRITE(6,19)XX1( I) ,YY1( 1 ) ,YM1 { 1) ,XX2( I ) ,YY1( 1),YM1'(I),AZ(I),F2(|) 
FORMAT(rF13.6) 

CuNTlNUt 

Cl'NTINUt 

STOP 

END 
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lOU 


lU 


2C 


N«i>0 

00 190 K*ltt>l 
FilKl«-0.2*}K-U 
G»UUA2«F2 )) 

F-«Al*A2«F2U)«(ti«A2«e*Ail 

GC*tMFtXlGtO.O) 

FC«CMPLX(FtU.O> 

ZC 1 1 I«TgL*C :»QRT CG(.»*2 -A. 0*FC 11/2.0 

iCl2)«IGC-CSQKTIGC*»^'-A.0*fCII/2.0 

ZtU-RtALUCUn 

ZR2«Ki:AL(ZCl2l) 

Z12«A1MAGTzU 2) I 
ZIl*AlHAG(ZUin 

xxuk1«zu 

XX24r|>212 

VVl(Ri«ZRi 

VY2(KI«ZK2 

V«BAfe 

W»0*A2«b*Al 

YY3IkT*R/V 

XXXiK >*-0.2*tR-l I 

cx»xxxm 

tY«VV3UI 

CONTINUE 

CALL PLOTXtO.U|O.Ot!>U 
CALL PL0TT5.0t6.0t-3) 

CALL FACT0R10.5) 

01*XXX*51)/10.1* 

CALC AXIS lO.CtO.O. *FbLUbACR PAR. • fl3f lO.OtO.O tO.C 1 01 ) 

CALL GRID (O.Ot-5.0tl0.1.0tlOtl.Ot-30S8A) 

Yl»AbSIYMAXi 

IPtABSlVMlNI.GT.Yl ) YTMBSIYHIN) 

YSTAR|*-VI 

UV*YT/S.O 

CALL PCOTtXXxm/OT tYYSIl )/0Yt3 I 
UO 10 l*2tN 

CALL PLUnxXXm/DTtYY3«II/DYt2l 

CALL AXIS (C.0t-5.Ot*PLANT PULLS AND ZLRU* t20tlO.Ot90.0tYSTART tOY 


30 


♦ ) 

CALL 
CALL 
00 cO 
CALL 
CALL 
DC 30 
CALL 
CALC 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
*£}*it* 
CALC 

♦ ) 

CALL 

CALC 
♦ o5 
CALL 
CALL 
STOP 
LNU 


NLHPLNO I 

PLOT (XXx{l)/UTtYYlU)/0Yf3) 
l«2tN 

PLUTiXXXm/OTtYYUl )/0Yt2) 
PCUT(XXX(U/0TtYY2U)/0Yt3) 


1»2 

PLOT 

PLOT 

PLOT 


N 

XXX(l)/0TtYY2(I l/O^'tZ) 

— 0 a 5 1 —6 a 5 1 3 I 

-0a2t-6a5t2 I 


NbWPCNtl) 

PL0Ti-0aSt-7.0t3 ) 

PCOT(-Oa2t-7aOt2) 

PCUT40aOtO.Ot3) 

SYMBOL 4-Oa 15 t-Oa 5t0. I5t *PUCbS*f OaO»5 ) 

SYMBUU-0al5t-7.0.0ai:- 

SYMb0L(-0a5t-7a5t 


__ I •ZtROSO.Of A) 

atnoui. t -Ua? t-*r a? t Jal5 a * PLANT * Y 4Kl«i AIaA 2 )VIK-l )-| Ai*A2 ) YfR- 
k)»U4R-l )-4Al«L«b<»A2)^lK-2)*tO.Oa7U 

SYMbULf-Oa5|-BaO tUal5t*CUNTR0Lt UIK I *F1*R I K) ♦FZAY IK | • tOaOt 2V 


SYMBUL(-0a5t-b.5t0.l5|*RLSUCTlNG 
)^Al»A2«PA(B*A24LAAl ) • 


#0.0.641 

SYMBOL (-0.5 t -9. 0t0.15t* PAR AMefeRS: Al»0.95 
L«0.01*t0a0tAV) 

PLUT(Lb.0.0.0«-3) 

PL0T^0.Ut0a0#^V99) 


CHAR. tON.: Z**2-Z*(A1«A2«F 


A2S0.2 
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Appendix B: 

Supporting Alebra for (3>31) ■> (3-43) and Sinulated 


Check 

Contents: p. 

Supporting algebra 

B-4-5 ; 

Step response test prograa 

B-6 : 

Step response natch 


1 


from ( 3 - 31 ) 




(v,« -f v^)l(e,,z 6 ,^)(* 


(t^^) ^ (6«,x B^«)(z 


a^z 4* a^z + a^z + a^z + a^x + aj^z 4 a^ 


Kivaluating the numerator 

N(z) - <fiz^ “ ^u“'’l*" ” ®ll' 22 ^ ^ ® 12 *^ ~ 

" ^ 12^*22 ^ *^ 21 *^ ” ^ 2 l‘'ll’^ " ^ 2 l“l 2 * 


^ 22 ^ "* ^ > 2 '^ 11 ^ "* ^ 22 ^ 12 ^ 


1<S( ^ *'*‘^'’^Hi“ 21 ^^12 “ ^2l'*ll '*’ ^22^^’’' 

+ (6(-B^jU 22 - ®i2"21 " ^2l“l2 " ‘'22“ll^^* 


+ l6(~tJj2“22 * ‘^22“l2^* 


yields ( 3 - 32 ) - ( 3 - 36 ). Define the denominator as 


D(z) • Dj(z) - D^Cz) 


where 

2 2 2 
Dj(z) ■ (z - rij^z - n2>(z ’ " ® 21 * 


D-(z) » (v, z + v„) 





I>2 


Expanding D^(z) yields 


Dj(r) . - nji - nj)<«‘ - «2jJ^ - * “u“2l‘^ 

+ *11‘»22* " “l2* ■*■ “l2"‘21* “l2'*22^ 

2 3 2 

* Cz •• t| j z •“ ^^2^ ^^4 "" ^^21 ** ^^22 ^*11^21 ^12^^ 

+(ujj.«22 ^ "‘i2"‘21^’' ■*■ *12“22^ 

“ ^"21 ■" “ll> - <“22 ' “ll“21 ■" ^ <“ll“22 “12“21^’' 

+ + '>1<“21 * + '>1<“22 • '>U“21 * ‘iz’*’ 

^l<“ll“22 * “l2‘‘2p^ ~ '’l“l2‘*22’' " '^2<“21 '*' “ll^^ 

+ H2(‘»22 " “ll“21 ■*■ “ 12 ^^ ~ '^2<“ll“22 “l2“21^* “ '^2“l2“22 


z + (-U2^ - aj^ - n^)z + (-a22 + * '^l“ll “ ^ 2 ^^ 

*<“n“2.> * “12“21 * ’’l“2:"''l‘"ll“21*''l‘>12 * ’’2“21* ''2“u>"^ 


’■‘ ‘l2“22 • "l“u“22 ■ "ri2“21 * "2-22 ’ "2-u“21 * "2"l2>‘'‘ 


+(-n .J2-22 • "2-U-22 • "2-12-21'' - "2-12-22 ' 


(B-6) 


Expanding 1 ) 2 ( 2 ) j^iven (B-2) yields 

1 ) 2 ( 2 ) * (VjZ + V2)((t<jj + ■*■ <"^ll“21 '*’ ®12 " ^2l“ll '*’ ^22^^' 


+ (- - 1'i 2“21 “ ^2l“l2 ” ^22“ll'^ 


+(- 6j 2''‘22 “ ^22“l2^ ^ 


B-3 


- (Vjfiu t Vj 021^*^ ^ ''l®12 " ''l ®2l“ll **■ ''l ®22 

■'■ Vll ■*■ ''2 ®21^* ■*■ ^'^l®ll“22 ” ''l^X2®2l " ''l ®2l“l2 
" ^1 ^22®11 " ''2®ll“21 ■*■ '’2^12 " '*2 ®21®11 ''2 ®22^*^ 

+ (“Vj6j2®22 “ ''1 ^22®12 ” ''2®11®22 " ''2®12“21 

-V 2 621 “i 2 ” ''2 ®22“U^* '*■ ^”''2^12®22 ‘ '^2 ®22“l2j* 

ttsibtracting (B-7) form (B-6) as in (B-3) yields (3-37) - (3-43). 

For example 5 the unit step response from zero initial conditions of 

Y ( z ) 

in (B-1) with the parameter definitions of (3-?.2) - (3-43) was compared 
with that of (3-17) using the control law of (3-25). The program listing 
and outputs for this comparison complete this appendix. Note that the 
outputs are identical, which verifies (3-32) - (3-43). 


B-4 


UOb 


IJ 






) fVCUOO 
/u,r./,GMl 
/#k6/0,6/» 


WATFly»^P*2‘^fHAGf -lUO.TlML»U»0 
DIMiNSiCJN Gl(PliFl<?TtZLll?)»WI?lfAlC?»»G?l2»fFL'<2>,Ati: 

♦) # VC(200) .DVJ200) 

DATA T/O.E/fAL/l*Ot l.C/,E/l .0/t^E T/0 .2 f 0. Oi/ , W/ U. »,t.O/ 

♦/I. 687/ fGM2/-C.7*«l/,K/l. J/,YI/'J.U/,Y:/0.0/,V3/<'.( /|YA/( 

♦/f U2/J.O/fOi/U. 0/»OA/u,O/f Kl/U.O/ |K2/0.0/ t^3/0,0/ ff'A/L 
♦Ro/O.O/ 

00 10 I*lt2 

Gi( n* 2 *tx»M-ztn n*wm*n*cus 4 wm*T*swKii i«zlt n »** 2 ) i 

IM l.U.nAlUGl(l) 

IF( l,EC.2lA<:lsol{j{ 

in*-Expi-2*2trn i^wu I*!! 

IH l.tW.UAl^»(.*2lll 

1 »■ t = »A22=o2l^l 

FU I)s( AL U I/WI n ♦♦2»’»Cl-fcXM-2tT ( 1 l*W( 1 |* U ♦ ICl S C W I U* l♦SO■^ T ( l-Zt 

♦T n )** 2 » ) ♦! ZL II I }/:>;Km-ze II 1 »♦♦?) l♦SIN^ wiimi ♦S i.vii i-ztii n**?) i 

♦ I J 

iPli.Eo.mn-Piin 
IH J .EJ.2I' LJ U) 

F2ll»s|ALin/wm^'*2l*IEXPI“ZETm*WlI»*l»*l£XPI-ZlTII Uwl j)*n-cc 

♦ Sj W| n*T* 5 >v.Kl I l-z- I 1 1 }**2M ♦! 2t I I n/SA« 1 1 1-Zt II II **?n*LlNtWI 1I*T* 

♦s;kiii-zlti I 

IMI.Ew.l )LK=H^m 
IM I.cQ.i I I. I 

CLM l.NUf 

wf n 1 1 o , I & ) 

FORXAll lx , * A ii - , U»X,* AU» , U'X, • 1 1 1 » , 1 * • L- 1 ? • , 1 OX » • A 2 1 • , 1 ') 


*Kf*f22* tlC'X 


, • I : I • , H'X , • 


fFl lt 6 l< f/<-l tkt. 2 tC 


:1 


WRIlE(^»tG)A3 IfAl 
FOi^MATI//,. f 
•;T1=C*^1-Ali 

ANU1=I t U»A l’4|t I I*A1 All-Al ? *011 /a LI I / I ” 1 1 » A 1 1 ♦ ;« I L-A 1 

tK = All»t n4oNL-J>ll*A/.L'l~Ai: 

ANlu' = L T?^AU/:a? 

A'RlTb (ofir ) 

f iamak I't ♦iox,»( II* ,ijx,ni?',lox,*Nii*.i(X,»Mur») 

KPI7‘^lfa,20)Ltf TUlT^,ANU 1.ANU2 
FLr.MATI//,t-Fl '.ol 
L>U 3*5 Jt 1,^(0 

V = IA214AI I )♦Yl♦^A<^^-An♦A^|♦Ali)♦Y2-f Al I * A22 ♦A I 2* A? 1 ) ♦YD” I A1 2*A22 I 
♦ ill^t Ol,?! I ♦Ul~iail*A2l-ai2*t*y?l*A ll-F*822 I*U2-I f 1 1*A22+H 12 *A 
♦2l4£*821*AlL^l ♦t'22AAil)*U3-lbl2*A2^:^LwB22*Al2|*UA 
U = L ♦K+L 11 *U 1 4 t T ^ ni*.* frNUl *Y1 ♦ANL-L* Y2 
YC IJI=Y 
AUIJ|=U 
YA = Y3 
V?" Y2 
Y2= a 
Y1=Y 
UA = U3 
U 3 = U 2 
U2=U1 
01 -U 

con-tinul 

3 ^=^C*IB 11 *L*H 21 ) 

a4 = D»l-0ll»A?l>bl?-t*02l*AU*E*Ll?2) 

F3 = D*l-t ll** ^L»Fl2'^A2 1-fcAF21*Ai2-t*322*All) 

8 L = L'*l-blc*A 4 -t.-t*b.. 2 »Al 2 | 

31-0.0 


40 

60 


55 

50 


ll-6T2-ANUl*Hll-ANUl*E*B2l 
l*A2Uil l«Et|*An«ANU 

-ANUl* £ •B22-ANU2 *H U -ANU2'PE*&2 1 
!♦£ f?*A??-f T2*All»A2l^tt2*AI2*ANU 
•A12 aANUI*L»D2..*AIUANU2*HI*A21- 

!2*A2I«ANU1»B1 2*A2L'+A>lUl*t*B22*Al 
£*B21*A12«ANU2«E*L22*AU 
♦ E*b22^Al2 


Bo«ei 

A6«1.0 

A5*-A2l-AU-tn 

A4*-A22^AU^a;1-A12^E I1*A21«ET1*A 

A3*All*A22*Ali»A2Utn*A22-LTl*Al 
♦ l*Bll*A2l-AN0l*fal2^ANUl*i:*B2l*AU 
A?*Al’*A?2-l r l*All4A??-ttl^Al?*A? 

♦l*Bll*A22*ANUl*B12*A2l^ANUl»E*e2l 
♦AMU2*0124ANU«:*t»32l*All-ANU2*E*b2 
Al*-ETl*Al2*A22“tT2*All*A22-i:T2*A 
♦?^ANU?*bll*A?r^AMU2*ai2*A21»ANU2*l 
A0*-£T2*AU*A..2^ANUi*tl2*A22*ANU2’ 

Vl*U,U 
Y2*0.0 
f 3*0.0 
YAsO.O 
Y5*0.0 
Y6«0.0 

DO 40 K*l ,200 

Y*-A5*Y1-A4*Y2“A l♦Y3-A2•Y4-Al♦y5-A0•Y6♦B^•Rl♦B4•R2♦t3♦R^♦92♦R4♦8l* 

♦R5^bO*R6 

YO(K)*Y 

Y6»Y5 

Y5*Y4 

Y4*Y3 

Y3*V2 

Y2*Yl 

Vl*Y 

R6*K5 

R5*R4 

R4*R3 

R3*R2 

R?*M 

R1*R 

CONTIMUt 

MRITE(6,60> 

FOKMAK 1X,»YC: CUNTRULLED SYS*,5X,*VUJ UVtRALL 1. F • ,12X, * I YC-Vul ♦* 
♦2 • ) 

DO 50 L*1,?00 

OYIU«l YC(U-V;MLn**i 

WKJTE(6,55)YC(L) fYOU)tOYIL) 

FORMAT t lX,Fl3.6,i0X,Fl3,o,l0X,Fl3.6) 

coNmuc 

STOP 

END 


A21 


vet 




8—6 


All 

A 12 

8 1 1 

bl2 

1*895669 

-U.9U4037 

0.120319 

0.116359 

A22 

bTl 

b22 


0.909837 

0.0/0016 

0 .067121 


DU 

6T1 

fcI2 

NUl 


-1.523571 


NUi 


O.SODUOO 
CONtRULLkC 
O.OOUOOO 
O.OV516 5 
o.2;:9<.98 
0.^49P52 
0. 73199 3 
0.972126 
1.269010 
1. 50971 2 
1. 72996 
i. 990312 
2.0V3265 
2.213959 
?. 305132 
2. 356959 
2.915576 
2.9172*>1 
2.921 103 
2.923309 
2.392859 
2.365598 
2.399551 
2.325536 
2.307992 
2.273036 
2.263695 
2.291999 
2.226110 
2.222950 
2.206096 
2.207929 
2.203193 
2.198592 
2.205502 
2. 19V992 
2.205691 
2.206909 
?. 206796 
2.2197Q3 
2.212575 
2. 216139 
2.219666 
2.217291 
2.222535 
2.220832 
2.221516 
2. 229159 
2 221069 
2.223919 
2.22271 5 
2.221739 
2.223776 
2.221090 
2.222990 
2.222060 
2.220692 
2.222396 
2. 22099 0 


- 0.158069 


$VS 


VO: 


-0.O7l5t7 
UVLRALL 1.6^ 
O.UUUUUU 
0.095165 
0.229998 
0.999052 
0. 731993 

O. 972126 
1.269011 
1.509711 
1.729961 
1.990311 
2.073263 
2.213952 
2.305181 
2.356950 
2.9155/2 
2.917256 
i!. 928098 
2.923298 
2.3^2852 
?. 38559? 
2.399595 
2.325586 
2.307938 
2.2/3030 
r .26 1839 
?. 291936 
^•226102 
2.222949 
2.206080 
2.207929 
2.2U3189 
?. 198537 
2.205999 
2. 199988 
2.205638 
2.203902 
?. 206799 
2.^19701 
2.212573 
2.216133 
2.219666 
2.217290 
?.??2535 
2.220881 
2.221515 
2.4.29151 
2.221067 
2.22J912 
?.???713 
2.2.:1736 
2.223779 
2.221037 
2.222978 
2.222057 

P. 22063‘5 

2.22239 

2.220410 


-O, 47/284 0.5ft6t3 

I VC-V0|**2 

0.000000 

o.ouocoy 

0.700000 

o.ococro 

0,000000 

o.oooooo 

O.OOUOOO 

o.oooooo 
o.ocoooo 
0.000000 
0.0 00 con 
0.000000 
o.uooouu 
o.oonooo 
o.oooooo 
o.oooooo 
o.ooocoo 
o.ococoo 

O.UOUOOO 

0.000000 


o.ouoc on 
o.oonono 
0.0{»0(>0o 
o.ocoooo 
0.000000 
o.ooocco 
o.oooooo 
0.000000 
O.OOUOOO 
0.000000 
o.oooooo 
O.OOOOGO 
O.OC OC 00 
O.OCOuOO 


0.000000 

O.OOUOOO 

o.ooocoo 

o.oooooo 

O.OOOOOO 

o.oooooo 

o.oooooo 

o.oooooo 

O.OUOCCO 

0.000000 

0.000000 

O.OOUOUO 

o.oooooo 


o.oooooo 

o.ococoo 

o.oooooo 

0,0( OI'CjO 
O.OuOOOO 

o.ooocoo 


O.OUOCCO 
C. 000 coo 

0.900009 

0,0l‘0('0(i 




2.220990 
2.220992 
2.22099 0 

2.220991 
2.22099 0 
2.22C99 0 
2.220989 


. 22098? 
2,220979 
2 . 2 J09. 3 
2. 220979 
2.220981 
2.220*582 
2.22097V 


O.OCOCOO 

a . 000000 

0.0 cor no 

O.OuuOuO 

o.oococo 

0.000000 

o.ococoo 
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IMSL ROUTINE NAME - ZRPOLV 


C 

c 


COMPUTER 
LATEST REVISION 
PURPOSE 


- 1SM/OOU6LE 


- JANUARY If IVTd 


ZEROS OF A POLYNOMIAL MITN REAL 
COEFFICIENTS ( JENKINS-TRAUBI 


USAGE 

ARGUMEN1S 


•> CALL ZRPOLY tAfNOEGfZtlER) 


A 

NOEG 


INPUT REAL vector OF LENGTH NDEG^l 

CUNTAINING THE COEFFICIENTS IN ORDER OF 

decreasing powers of the variable. 


CUNT A] 


INPyf INTEGER DEGREE OF.POLYNOMI AL 
NDEG MUST BE GREATER THAN 0 AND 


THAN I 01. 

OUtpUl COMPLEX VECTOR OF LENGTH NOEG 


LESS 


lER 


CONTAINING THE COMPUTED ROOTS OF THE 
POLYNOMIAL. 

NOTE - THE ROUTINE TREATS Z AS A REAL VECTOR 
OF LENGTH 2*NDE6. AN APPROPRIATE 
EQUIVALENCE. STATEMENT MAY BE REQUIRED. 

SEE DOCUMENT EXAMPLE. 

ERROR PARAMETER. lOUTPUTI 
TERMINAL ERROR 

IEP«IZ9f INDICATES THAT THE DEGREE OF THE 
POLYNOMIAL IS GREATER THAN 100 OR LESS 
THAN I. 


1ER«130| INDICATES THAT THE LEADING 
. COEFFICIENT IS ZERO. 
IER*13lf^INDICATES THAT ZRPOLY FOUNI 


^ FOUND FEWER 

THAN NDE6 ZEROS. IF ONLY M ZEROS ARE 


FOUNDt Z(J)jJ»M4I|...,NDEG ARE SET TO 
POSITIVE MACHINE INFINITY. 


PRECIS I ON /HARDWARE 
REOD. IMSL ROUTINES 
NOTATION 


SINGLE AND D0UBLE/H32 
SI NGLE/H36 tH46 tH60 


UERTST.UGETlOfZRPQLB.ZRPQLCfZRPQLOfZRPQLE 
2RPQLF fZRPQLGtZRPQLHfZRPQLI 


INFORMATION ON SPECIAL NOTATION AND 

CONVENTIONS IS AVAILABLE IN THE MANUAL 
INTRODUCTION OR THROUGH IMSL ROUTINE UHELP 


COPYRIGHT 

WARRANTY 


- 1978 BY IMSLt INC. ALL RIGHTS RESERVED. 


IMSL WARRANTS ONLY THAT IMSL TESTING HAS BEEN 
APPLIED TO THIS CODE. NO OTHER WARRANTY* 
EXPRESSED OR IMPLIED* IS APPLICABLE. 


SUBROUTINE ZRPOLY 


INTEGER 

DOUBLE PRECISION 


INTEGER 

REAL 


(A*NDEG*Z*IER) 

SPECIFICATIONS FOR ARGUMENTS 

NOEGflER. 

Auitzm 

SPECIFICATIONS FOR LOCAL VARIABLES 


REAL 

DOUBLE PRECISION 


NjNN*Jj|JJ*I*NMI.ICNT*N2.L*NZ*NPI 
EtAtRMRE*RlNFP*REPSP*RAOIX.RLO*XX*YY*SlNR. 
COSRtRMAX*RMlN*X*SC*XM*FF«OX*OF tBNOfXXXfARE 
PTI lOU 


1 

2 


DOUBLE PRECISION 


TEMPI 101 )«P (101 )*QP( lOlItRKdOI )*QK(101)» 
SVKdOll 


LOGICAL 

COMMON /ZRPQLJ/ 


SR»SI*U*V*RA*RB»C*Of A1*A2*A3* 

T|ZA»§l!c^?^A£foR!REp£RI*zlRb*^ 

P»QP*RK*QK*SVK»SR|SI*UtV*RA»RBrC*DtAltA2*A3*A6* 

A7»EtF*G*HtSZR*SZi«RL2R*RLZI»ETA*AREtRMRE*N*NN 


»FN 
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fHE 


iHi.PIIOCMM* THt HEANIMS 


^00| J^pSTANTS AR| - HEPsSl THt 
m SnAUtSt ROSlflVE FLOATING 

£glSfinSS8 r*^ 

RINFF |HL LARGEST FLOAT 1NGH»UINT 
NUNS ^ % 

•'^?a.iviL«!ftvisf*jySirii!SHs 

IXFONENT RANGE DIFFERS ' 


IS 


AM) DOUBLE 

AND HINFF SHOULD INDICATE 


DATA 

DATA 

DATA 

DATA 

DATA 


lEK * D 
IF (NDEG 


^ IN SINGLE 

FRECISION THEN REFSF 
>HUULD INDICATE THE 
SNALLEH RANGE 
RADIX THE^BASE of the FLOATING-POINT 
NUMBER SYSTEM USED 
KlNFP/iTEFFfFFF/ 

REPSP/ZOO 100000/ 

RADI X/l 0.0/ 

REPSRI/Z3AIOGOOOOOOOOOOO/ 

ZERO/O.ODO/tONE/l.CpO/ 

ZRPOLV USES SINGLE PRECISION 

CALCULATIONS FOR SCALING* BOUNDS 


AN^^ ERROR CALCUL 


ETA 

ARE 
KMRE 
RLO e 


>GT. lUO .OR. NDEG 


REPSRl 
ETA 
' ETA 

HEPSP/tIA 


FIRST EXECUTABLE 
•LT. 11 GO TO 16S 


SCALING, 

_ATIUNS. 

STATEMENT 


XX « .707 1068 
VV « -XX 
SINR » .9V756AI 
COSR * -.069756<*7 
N » NDEG 
NN B N«1 


IF U(1 I.NE.ZERO) GO TO 
lER « ISO 
GO TO 9000 


» IF (A(NN).NE.ZERO) GO TO 10 
J « NDEG-N^l 
JJ K J4NDEG 
Z(J) * ZERO 
2(JJ) « ZERO 
NN * NN-l 
N * N-l 

IF INN. EC. U GO TO 900i» 

GO TO S 

10 DO IS 1>1,NN 

P(I) ‘ Ail) 

15 CONTINUE 




CONSTANTS FOR 


ALGORITHM FAILS IF THE LEADING 
CDEFFICUNI IS ZERO. 


REMOVE 

ANY 


The zeros at the origin if 


20 


IF 

IF 


(N.GT.2) 
(N.L1.1 ) 


GO U) 30 
GO TO 9005 


MAM A COPY OF THE COEFFICIENTS 


START THE ALGORITHM FUR ONE ZERO 


CALCULAU THE FINAL ZERO UR PAIR OF 
ZEROS 

IF IN.EQ.2) GO TO 25 
Z(NDEG) « -PI2)/PU) 

Z(NOEG^NOEG) « ZERO 
GO TO 195 

25 CALL ZRPQLl ( P ( 1 ) * P 12 ) * P (3 ) , Z t NDEG-t ) , 2 INDEG^NDEG-U ,Z iNOEG ) , 

I ZINDEG^NDEGI) 

GO TO 195 

FIND LARGEST AND SMALLEST MODULI OF 
COEFFICIENTS. 

30 RMAX « 0. 

RMIN * RINFP 
DO 35 1*1, NN 

X * AbSisNGHPi 1 ) n 
IF (X.GT .RMAX) RMAX « X 
IF (X.NE.0..AND.X.LT.RM1N) RMIN * X 
35 CONTINUE 
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c 

c 


c 

c 


c 


c 

c 


SCAU IF THERE ARE LARGE OR VERY 
OEFFICIENTS OF THE FOLVNOHIAL. 


, _ OF THE FOLVNOHI4 

OvlRlLoii^ ANO^ lo^vS iS %nK{!{EcT EO 
UNDERHOH IHIERFERING WITH THE 

THl°?l!l?S"il THE 


RASE 


SC « RLO/RRIN 


IF ISC.Gl.l.ti) GO.TO.AO 

If Irmax.11. 10, 


IF ISC.EQ.O.) « 

GO TO AS 
IF IRINFF/SC.LT.RMAXI GO TO 

FACTOR « OftLE (RADIX )«*L 
«l,NN 

FAC10R*P( II 


GO 10 t>S 

REPSF*RA01XARA0IX 


40 

45 


50 


55 

<>0 


00 50 1«1,NN 
Fin s 


60 m.NN 

FI (II • AeS(SNGL(P(in I 
PTINNI « -PKNNI 

COMPUTE UPPER estimate 
X » EXPIIALOGI-PTINNII-ALOGIPKIII l/NI 
IF (PTINI.EO.O.I GO TO 65 


55 


COMPUTE LOWER ROUND ON MODULI 
ZEROS. 


OF 


OF ROUND 


C 

C 

C 


65 

70 

75 

ao 

85 

90 


XM 

IF 

XM 

FF 


* -PKNNI/PKNI 
(XM.Ll.XI X XM 


X*.l 

PI (II 


IF NEWTON STEP AT THE ORIGIN IS 
BETTER, USE IT. 

CHOP THE INTERVAL (0,Xl UNTIL FF.LL.O 


DO 70 I«2,NN 
FF « FfAXM^PKII 
IF (FF.LE*O.I GO 
X * XM 
GO TO 65 
DX » X 


TO 15 


IF (AbS(OX/XI.LE..005l 
FF « PI II I 
OF * FF 

DO 85 nZjN , _ , 
ff. m FFAX^PIIII 
DF * DF*X^FF 
CONTINUE 
FF * FF»X4PT(NNI 
OX * FF/OF 
X * X-DX 
GO TO 80 
BND * X 


DO NEWTON ITERATION UNTIL X 

CONVERGES 10 TWU DECIMAL PLACES 

GO TO 90 


NMl 

FN 


s N-l 
ONE/N 


COMPUTE THE DERIVATIVE AS THE INTIAL 
K POLVNf ■ 

NO SHIF1 


K^PQLYN^MIAL and DO 5 STEPS WITH 


95 


JO ^95 1«2.,N 


II * (NN-I |4P(II«FN 
RKIil « Pill 
AA * FINN I 
BB « PINI 

ZEROK s RKINI.EQ.ZERO 
00 115 JJ«i,5 
RK(N) 


C 

C 


CC 

IF 


(ZERUKI GO TO 105 


100 


I * “AA/CC 
DO 100 

J * NN-j 
RKUl e 
CONTINUE 
RKIII > Pill 

ZEROK « 0ABS(RK(NII.LE.DABS(BBI*ETA«10. 
GO TO 115 


USE scaled form of recurrence if 

VALUE OF K AT 0 IS NONZERO 


»NM1 

TPRKlJ-mPI Jl 




•*.' 


rrn nonrrio no r*t^i 


C-6 


I 

C 


105 

no 

115 

120 


OU 110 1»1|NM1 
J • NN~l ^ 

^ HMJI « RKIJ-1> 
CONT INUt 
HKIU ■ 2kHU 

^ Z|R0K « RRlNI.EO.ZfeKO 
CONTiNUt 

00 120 I«ltN 
TEMPin « HKIU 

00 140 ICNlKltZO 


XXX « CU5R4XX-SINR4YV 
VV « JilNK*XX4C05H*VV 
XX « XXX 
SR > BN04XX 
SI « BN0»YV 
U * -SR-SR 
V « bNU«BNU 

CALL ZHRCLb (ZO^lCNIiNZI 
IF <NZ.bO.OI GO TO l30 


L 

C 


J * NDtU-N^l 
JJ « J^NObO 
ZIJ) K SZK 
2<JJ) K SZ! 

NN ® NN-NZ 
N ^ NN-l 
DO 125 l»l.NN 

125 PUJ * 0P<1 > , ^ 

IF (NZ.bO.n GO TO 20 
2< J«n « RL2R 
ZiJJ^U • KLZl 

GO TO 20 


150 DO 135 I«1(N 

U5 RK<n UMP< U 

140 CONI INUb 


lER « 131 

145 00 150 l^^ltNDbG 
NPI* NOfcG^l 

p<n * ziNPu 

ISO CONIlNUb 

N2 * NOtO^NObG 
J K NDIG 

DO 155 l*l,NDtC ^ 

2<N2-lI » ZU» 

Z(N2) P<J) 

N2 * N2-2 
J c J-l 

155 CONTINUE 

IF HER .to. 0) GO TO 9005 

N2 » 2*<NDtG-NN»+3 
DO 160 l»ltN 

Z<N2) « RINFP 
Z(N2tU • RINFP 
N2 « N2^2 
160 CONTINUE 
GO TO 9000 
165 lER » 129 
VOOO CONTINUE 

CALL UERTST < UR t6HZRPULY) 
9005 RETURN 

bND 


use UNSCALED PORN OP RECURRENCE 


SAVE K PUR RESTARTS WITH NEW SHIFTS 

LOOP^TO SELECT THE (lUADRATIC 
CORRESPONDING TO EACH NEW SHIFT 


QUADRATIC CORRESPONDS TO A 
SHIFT TO A NUN-REAL POIN 
COMPLEX CONjUGATE. THE Pi 


.DOUBLE 
I.AND I 
OINT HA 


MUUIILUS BND. AND AMPLITUDE RUTATEO 
BY 94 DEGREES FROM THE PREVIOUS 
SHIFT 


SECOND STAGE CALCULAT lUNt FIXED 
QUADRATIC 


THE. SECOND STAGE JUMPS DIRECTLY TO 
ONE OF THE THIRD STAGE ITERATIONS 
^ AND RETURNS HERE IF SUCCESSFUL. 
DEFLATE THE POLYNOMIAL. STORE THL 
ZERO OR ZEROS AND RE TORN TO THE 
MAIN ALGORITHM. 


U THE ITERATION IS UNSUCCESSFUL 
ANOTHER QUADRATIC IS CHOSEN AFTER 
RESTORING K 


return with FAILURE IF NO 
CONVERGENCE WITH 20 SHIFTS 

CONVERT ZEROS (Zl IN COMPLEX FORM 


SET UNFOUNO ROOTS TO MACHINE INFINITY 


«B>iS'3r fy\ MS nm'sa m imi lUirtss 
w«;r« 2 prM m PI5 \m t)P (•r.c 


p>» £9Uo.Hn^ p9li>t» >>irv» Con t'w llbory :»lt!:lQ9 .»i'>ll9'nl In tV? 

Ii>L ttbr.ny InCoriMt Ion joimnl, 1. R>? I'-VIL routlnsi nro prni^r idtnry* v»l 

tlv3 uTstrTmitnr hn* plnol UbrToUovlni roitcl«tU»ni »>n V\a\t 

* Oi»ly Intornil mn;ja nt ona nlta In nUoiiirs>t. Hf;L pr'vU>:t» my nnt 

IntJf.Hlamily ba .mla ivnil-ibl j t > ot'uc onfxitulotvil oltci. 

^rra(\ 3 anants with in'iit ba if ip.JUc.itlan pro|r«ia3 sontnlntn] anlyrtlal 

aodn Crotu I41L library rmitlnss .tro ta ba Jistr lbvita«l ovitsiilo VlritnU Taah. tniL# 
tns., has th> Intant af balaj ajnjaablo In this casa IC th? appUsation pfojraafi nro t> 
ba «iia a/allnblo it na chir|o» b\it tha (ilstrlbutor Jaslron Hnnwl'xlia af »ir:b 
'Slstrlbatloni. IC t!w projruu ii‘2 t.» bv< mb ny/allnbla Cor a chirji, aontrvst 

• torn IF-N00r») covers thj rai'ilral financial arranjsnont. 

* When rint, library sobrovitlnat are unxl In rosairch »vjck ml rasaltlnj ptibUtatloni 
roiolrof In th» opinion oC t)u aottwr, tivi reproilurtlon oC listings of psrtlnmt 
library source cole In thi rosearch ixibllcatlon, tMlL rc'|>*3sts rofjronne# a»vl jrmtn 
:>er!ni salon. It only the julifoiitlntj <k) 2 anant>tlon la nijusiry or a.iCClclnt, l^*L 
on'Tourajes thit mob uC usi. Vun a.^alluatlon proirans djvolopKl n a by-?rolu:t >C 
rosearsh v<rark contain inbsJlel library routines » theaj a^>pllcstlon pro] ram my be 
taken, is Is, by ths autlwrs on thslr «bpirtnro fron the rosoarch Institution Cor uaap 
lt» their now locations. H'it’a only restriction la that any anbeittlai library rootlnen 
bo use‘1 only as orl jlnelly InrjrpjratO'l In the applUatlon projraa a.vl thit they not be 
re;iV)/e-l for usije In any other seCtwire de/elopnont. 

* l*r>t soCtwira .ftsy not bo tmvjCnrcJ to other conpiAtera wltlvjut H'JL’s ipeelflo (by 
c>’b nmol pernlaslon. 

In mmy cases ^raJuite utolent * uie library olenonta to altl In satlsCylnj rasenrch 
ro.jolronents Cor their cbjroea. IC, In tim opinion of the author, Ustinia oC Hit 
library source cole aist bo roprv>le:ol in the resultwiet thesis or dissertation In order 
that the rosults bo dopllCvatablo by otivera, H'.»L irants peruHtslon Cor such reproduction. 
In coses of sueli 'isajo, HJL rev>?sts that the author nako the follov/lnj points In the 
research docunent In a .mnner noticeable tv> sooeone li^tareste.) In the llstinjs: 

* r.»o listtyl projrvin Is part of a proierlotary product bolon]ln] to Intirnatlo.nal 
»^athe.natlc3l and Statistical Literal las, Inc., In ihuston, rvj<a*. 

* The llstlnjs aro raproJueoJ wltt) the pernlinlon of I'Ut. 

* Vive llstlnjs may not tea oxtracuel Cor oUver pirjeoios or uaal ai the ?>asl3 for any 
U'ftvaro de^alopnent. 


I ha/u reed the ebo/o restrlctlonJt on I ajrae to abl b by then, Ploa.ie 


furnlsa e copy the IdjL routine 


Ma«e 


j \tj b^jpertiwent 



( ^A(t<LA 


Con.'xiter Center .accoiwat I to be tjiilel Cor cost oC run 7 . ^6 uffi^ 

'rtec charge for this sendee is the cost of the noi plus $2.00. 
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otiiGH^AL Page E 3 
OF POOR QUALITY 




C -9 


UNI1 ClKtU: 

irtis j*hu( 5 Kam 

SVSUM iN n 


A UNIT tlKCtt: AND LOCAlt^ »>UUb C Zi-HUi. OF A 

>X 3 (t)DQl» fX^(D 0 l)l)fVilbD 0 U»V^U« 00 ntV 3 ^ 
\AV 2 i 3 » 



DATA 

DATA 

DATA 


XZ/-D.lAAU;>,-t'.lAAl.: 3 #~D. 9 t- 9968 / 
• ~ ■^^ 401 * 0 . 0 / 


V 2 /t. 9 AL' 4 Olf- 0 , 9 A» 

N«U0C0 

ymax«i.u 

VM 1 N« 5-1 .1 

DU Z!> i>-<.«buDl 

X^ll >sX^'U~l |~U.U00^> 

V^n»--Ali.lStKT U 6 .C'-X^« 1 » 


J^OuOU”* I 

vi(jU“T;n~ii 
xut jU-x^i 1-11 
Yj(J>*¥. Il-l ) 

X9( 1 »S-X.M 1 1 
YA< 1 1 = 1 ) 

kb tUNTlNUL 

V.AU L P t UT b 1 0 .0 I C .U , bl »D 1 I 

CALL PLlJMb.L.t .U»-4 ) 

CALL pACTDKie.CM 

ASTAKls-1,1 

L'1-l/'*.C 

CALL AXlb I — *i»c»0*U»*KtAL*f"^«»P»C|C»C|Xi>lAKltt'lJ 
VisAbbI YDAXl 

YHINJ .U1 .Yf » YUAbblYMlNl 
VbTAKI a-Yl 
UY=Y t/H.C 

CALL AX lb IC.C i-A.g, • 1 f^ADlNAKY * , 9t U . C# ‘K. .0 tY bl AK I ♦ L Y > 

CALL NLWPhNl-U 

CALL PtiiTIXl (1 )«Y1 t Ui4) 

t'U 3D 1=4»N 

4U CALL PlliHXUntYll 1 
CALL PLi.Tl X2ii)»Y3l 1 
lU 3 d l=?,N 

JL CALL PLl 1(X3intY3il ii<’) 

CALL PLUTlX3mtYil llf'il 
LU HD i=«.»N 

HD CALL PLDK X3(l ) tV3l 1 ) »3l 
CALL PHiT(XAiI)»YH( l),3) 

Di.i be Is^fN 

bO CALL PLuTlKHlUfYAn )»3> 

DX— C»«,b 
UY» D.I:t) 

UU ICO lal,A 

Call LYhHuLi xpm/uXfVpm/:uY,f'.3» H,o.Of*i > 

Uv‘ CUNTINUt 

Du LCu J»l*3 

AXi < J)«XZt JJ /DX-U.Ult)lb*A,C 
AYZ ( J) = YZI J 1 /DY-0, O: 7 . t 

CALL LYMBULi AXZLJI |AY2(J) |U.^tll2f U.Ct-ll 
3 Cd CUkllNUt 

CALL PLUTl Ifc.U.D.O,-^) 

CALL PLut I D. D» c, Df W9 ) 

SlUH 

tND 




C-IO 


c 

c 

t 


2 *> 


i-'O 

30 

<♦0 

b\i 


100 


200 


UHAW!k A UN11 ClKCLt AND LOCATE PULES C ZEKOS OF A 


UNIT C1HU.&: 

THIS PRUCKAM 
SYSTEM IN lie 

DIMLNSIUN Al|DiiOU(A2(OOUniX3Hi)Ul).X<i|bOontVi(OUUntV2il»UUntV3 

♦ IbOOll tV<*(bOOU,XiMo)tViMo) tX2(Slf V2lLi)tAXnt>) tAV^I^I 
UAIA Xl( l)tX2U|f X3l U tX4( l)/O.OtA»OtU,Ut-<«.0/ 
data YltntY 2 ( litY 3 ( l »» VA ( U /^. Uf C . O »-'*. OtO . O / 

AP/0,lVVt03t-0..'bV02^,-0. 72 V^^St*"C. 72 V7*»!> fC ^0 fO .8^23^0/ 

YP/O.D*w.U»U.ODO'772f-U.6U07 72 tO. I 71 lb I «-U* I 71 1 01/ 
A 7 ./ C # O » 0 . Ut - O . lH ^ l 2 5 »- w * IhA IZit-U.^tjVVCV/ 
Y 2 / U . 0 | 0 . UtU . V « ciDlt *‘ 0 . V < t 2 :»Ul «( . 1 / 


data 

U7l A 
DATA 

data 

DATA 

N=but0 
YMAAftl.O 
VMlNs-l.U 
DO Zt> l*2fbOi l 
Xkll laXAli-i |-0.»^DDt 

V2III=-AbiCSUkmo.0-X2iil**2 II 

Js{jOC3-I 

AUJ)«X2C i-1) 

YHJ)*-T2U~I) 

X3(J|s-X2(l-l) 

V3U) = Y2U-1 I 
XAU >*-X2in 
Y4U I--Y21U 
CUNT INuL 

CALL PLUTSlU.UtO.Of JOJII 
CALL PtUUS.0»<»,0,-3l 
CALL FACTDKtO.bSI 
XSlAH|s-1.0 
Dl =1/A.U 

call axis 1— <,UtU, 0, »«t AL» t-HtO. Cf 0,0 *AST AKT ,U1 » 
VT=Abl»ntKAX| 

IKAbSlYMl-^1 .oT.Yl » YT sAliS|YMI.‘7) 

VSlAKT^s-YT 

UV=YI/A,0 

call axis iU.Uf— t.O, •1MAD1NAKY»» V,«.OtVO.UtVSlAKl»DVl 
CALL .4twPr>Hji» 

CALL PLUKXltUtYimtS) 

DO 20 l-4» i 

CALL PLUliXiU )tYUl)»2l 
CALL PLUIIX2U Ky2U )»3) 

DU 30 i=2tN 

CALL PLt.r IA2tUfY24i If 21 
CALL PLuTlX^ti )fY3( lit 3) 

DU '♦O l = 2fN 

CALL PLLi tXj(l |fY3( 1 It2l 
CALL PLUl 1AA(1 ItYAll )t3) 

DU SO l»2fN 

CALL PLuUxAU)fYH(i Ifl I 

D X • 0 • 2 S 


DY=C.2S, 

DO loO 1 = I ft' ^ . 

CALL SYMbOL 4 XP 1 i I /DX t Y P U I /uV tO, 2 t A 1 0 .0 1 “I I 
CUNT I >»Ot 

UU 200 J-ltl) . ^ 

aX 2 4 J)=XZ4 JI/DX-D.0lslS*A.C( 

AY24 JI-Y24 JI/DY-0.02 72 7*A.O „ 

CALL SYMbJL4AX2IJ)f AY24Jlf 0.2f ll2fC*0i-l I 


CONTliOt 

CALL PLUT4 Ib.D fO. C f-3) 
CALL PlUT 40. 0,0. Of ♦‘7*79 I 
ST UP 
LNO 


Appendix D: Overall Transfer Punceions for Exanple 6 


Contents: p. D-1-2: 

D-3-4: 
D-5-6: 


Supporting algebra 

Simulated step response check for (3-46) 
Simulated step response check for (3-61) 


Note : 


a*8 and b's in printout correspond to p*s and q's In (3-48) 
and (3-61) 


D-1 


Consider the control cooflguretlon In Flg» n (or c) with a convergent 

identifier, l.e. the time- Invar lent, aeyieptotlc case. For the plant In (3-17) 
with the controller In (3-46) and the actuator of (3-45) with 

then, since 


and 


Uc<z) 

6 K(z) + 

(Vj z*^ + V2 2"^)Y(z) 

(1 - 

“1 ,~2. 
z - H2 * ) 

liil- 

« UilL. 

ItsL 


U^(z) 

C 

U(z) ♦ 

Y(z) 

(1 - Hj^z’ 

- H 2 z"^) Y(z) 

U^.(z) 

6 R(z) 4- 

(Vj^z“^ + v^z”^) Y(z) 


^12 , ^21 *^22 

(1-8 

2 ■*‘2 

z - g 

Lz - a^^z - a ^2 ^ - “ 21 * ‘ “ 22 - 

1 


(D-2) 


(D-3) 


(D-4) 


Cross-multiplying in (D-'4> and rearranging yields 


Y(z) 

R(z) 




2 2 2 
- g)(z - ” “ 21 * ” ” 22 ^^* ” ’’l*' ’ 


^(z — B/v''- ” “ll*" ~ “12^'^ ~ '^21** ~ *22^'*" ~ ''l*" ~ ' 2 ‘ 

2 2 

-(Vj^ + V,)(l - g) UBj^jZ + * “ 21 ’' ” ““ll’'”"!^^ 

(1 - g) N (z) z^ + q^z^ + q^z^ + q 2 Z^ + qj^z + qQ 

(z - g) Uj^(z) - (1 - g) D 2 <z) Ip^z^ + p^z^ + + p^z^ -f PjZ^ + p^z^ 

+ PjZ ■»■ PqI. 


wiiere N(z), Dj^(z), and DjCz) are from (B-2) , (B-4) , and (B-5), respectively, 


Note that the of (D-5) ore equal to the (i-«) of (3-32) - (3-36) oi 
Bhown la (3-A9) - (3-53). From (b-6) and (B-?) the of (3-54) - (3-CiO) 
are readily fonned. 

1 

Similarly for tlu‘ control con fi Herat ion of Fig, 3-l2b after identifier 
convergence, from moving the pickoff point acroan the actuator 


U^.(r.) • 5 Hz) + (n^( ) «““1 U^(r.) 

■F [v^ + z"**! Y(z) (D-6) 

and 


Y(y.) , / l.“il X 
U^,(2) ^ z-g ^ U(z) . 

Therefore 


|z^- gz^ - n, (1-g) ^ - n,, (l-g)l Y(z) 

- - (1-g) 

Hz) + [V.R 4- V.J Y(z) 

X. ^ 

(D-7) 

Again croHS-multiplying and rearranging yields 


"ll" * "l2 


'^21“’ »22 
2 

z - a , , z - 


I(s?) 

Hz) 




(5 (l-g)UBjj ^2 + “ ‘'21“' “ ^’‘22^ ^ fiooHz* - vijjZ - Oj,)] 


r 3 2 2 ^ 

i(z — gz — ( 1 — gl^ — — ^^ 11 ^ ** ” ^22^ 


tv^z -F vJ{l-g)l«ilH -F - “ 21 " “ ‘^22> *• <^21" "22^'"-‘' ’ 


. (klilLNliL). 

/ 3 2 '* ^ 

'(z - gz - Hjid-g) z-n 2 <l-g)](z - rtjjZ - “ ‘*21*’ '* ‘’*22^ -(l-g)D 2 (z) 

, i 

(I)-8) 

where N(z) and D.,(z) are. from (B-2) and (B-7) , respectively. Kxpansion 
*• 

yields (3-62) - (3-68). 


0-3 


lU 




70 


^0 


10 , 


data r/U.t>/fAL/i.O#i.l'/ f -£b1/U.ti( .0? /fH/l*.‘>|6»*(./tU/t.»»/.Z£lA/l 
UATA Vl *VZ#Y^tV^t7»»UifU*JfU3| U** iM^/0#vf l)*0|v*0t O* U#0»l**0«0f0*if#0*0 

UATA yt f YEif Y(,7/0.u«a.uta*0t0.0tiuut0*0t0.0/«KifH 

3 f Kh |Kb^Kf» *k 7 »H/ i>« 0 f V • 1 * 1 1 < • C* f C » 0 1 0 •(' * 0 • I* # 0 • t* t 3 • E«/ 

H?uV*l*cxf I-/K 1 U I 1*1 »*ou.‘jwn u 

iFU.! «.l' >A^l*C U2» 
ui:m«-ixp|- 2 »zHii i«^wu >*! I 
iFu.cO.UAi^ 

! Ft I *c:g.z )A<.4.n><.t<. ) 

♦It 

*> t 

iM i.tw.ni tn 

n«i Au* uK>iaU{*.')»uxiM-ziin )fwtn;i tyiixFi-zkYu jfwt n*i >-CL 
— - ,Tt u/i.aKm- 2 iM n*Y.!M*siNiwn »*T* 


( 1 i»f AL U »/»^l i >♦*►4 > '»U-LXF1 -7tl t I » •Wl n*l tctt: twt 1 t l-Zt 

U»»ZM4l7tUl t/iOKU l“7tn I ♦t»lNlWtl )*Mi.wKll l-Ztf tl 


4i.twm*iyt>oKi u-7£i iu*»4iM*t zt 
♦:»«Km-7i.Tt u'»»m > > 


lFn.tU.ni;l<‘>F«tll 
If tl.tU.ZH ;^«.-F4^tii» 

CUNY INUt 

0Ml*^*bXPt-Z6lA*Wll»*1 J*OOStWUI*Y*!>OKl 1 1 .0-ZL I AW*^' ) I 
0 MZ*-hXPt-L*zLiA*wt n*i I 

EU»uMl-All 

ANUl = iLlltAU4 it I l4Ai |Al AU-Al2*t)ll/Ul2> >/ Ibl l*All4bl^-Al 

Ll2»All*LlltOK<:-bi4*ANOl-Al.; 

ANU 2 =t;n’*AlL/ll«: 

0»1>PI-M0»1 » 

Ht> = U*( 1 .0-0 )*■ it* H ♦b^ I ) 

B«*sD*( i .C -0 »♦ i-l‘l l*Ai. l4b,U-b2i*AU4b<2) 

»j>*0*il .0- 1 |4 i-t ll*A..fc-tJl2’*A0l-»2l4Al2“B22*An ) 

B2»L»( I .u-t l‘*i-L l2>*-/ tL-y.fi ^/U ) 

L I s U . 0 
.Uttt I 
A?*l .0 

A 0 -- 0 -At l-A U-t I I 

A t> s i t » t A 1 1 ♦ A I 1 ♦ 1. ? I > - A A 1 1 v^A 1 - A I 4 ♦ t 1 1 4 A J 1 ♦ C 1 1 *A 1 1 -I 12 | 

AAs i-0* i-A.,l4 AU»*AJ l.l4 t- 1 14 A,. 14 1 1 i»Ai l-L U I ■* A H V7 22 ♦ / I i’4« A J l+ll 1 

♦ 4A24-U I'fAil*#.*: j 41: n4Al2*l 1 44 A2i4U 2^X1 l-ANUl<‘U 1 1 *\ 4 I »»! 1 .0-0) ) 

A3* l~t.* i A U*A. ^ 4A vA 2 l4L ! IVA 4 . -I- I l*All fA.. 14L 1 1 fA 12 1*-*A2 !♦! I 2# A| 

♦ I ) ♦ A12* / . . -L I 1 ♦»* I I 4 A 2 I.-C 1 i *» \ A,. !♦? 124 A 2 2-fc \ d* f I 14 A2 !♦). I 2*A 12-ANu 

♦ 14 i-‘.ai < a 4 i*A 1 1 ) + i 1. w-o)-A;vt.2 4tb I ui r i ) <*( i .o-o ) ) 

A2» i-V.4i Alt I 14 a 1 14 A 22 -b 1 14A 12*A2l ♦ 1:1 2» A22 -t 1 2 I l*A2l4l I 2*Al 

♦ 2 )-t U*Aic>tA2t -fc I24AI 14A2t“t I L 4 A I 24 A t !■» ANU !♦ ( F 1 1 *> At * •» b 1 2 ♦ A 2 1 41 1 1 4A 

♦ I2*t 2t4An )4i 1 ,u-0 )-ANI*4 4 i-l*U 4A2l4!‘224Ll2-b2l4Al 1 » >1 1 .0-0 M 

A l» lw.4| t I l4Alt4At2-*l.‘i .,4Al l*A 2 2*t 12'4A 124 Al l I-L12 4A 12»/2l’4ANu2*ibll» 
♦A22 4ttl2 4Atl4Uit l4„lt*bt2*An )4 n ,U“l ) ♦ANUl»^t* l44A22■^^'.•2♦Al2»4^1.0- 
♦0) I 

AU=0*LT24A K>A2. ♦Al4024(31.i4A2t4tt224Al2 ) ♦! l.C-t ) 

wRll U(6 »e.u )i b ,i*H»uife; ,b i ,i»o 


i-UKMAl i iXt * 1 1 >* • t i V.b, tXt • t HV, t,f 2 X|* I A = SFV.t I tXf 'bZe* tFV. 6 » 2 Xf 

♦ *bls • .f y.O.tX , *t 0-* .FV.fif 2X ) 


WRilttbf70)AytAo.Ah»AAfA^tAt|Al.A0 
f UKMAH lx, •/,/=• tFV.C,2Xi 'Ats* Ab- 


♦ *Ai = *fFV,6f2X|*7*; 


• .f-V.5.2X Al 


•ffV.t , 2 Xt*A^s*,Fy.bfiX, 


»^KlTL(bt30)Lr ltl.l2tAtMiJl»ANU2»OMl,OM2 


• fF9.6,2X,*A0 = » ffV.b»2X) 


MU 4 : f un i t 

• ’ti i~’.rv.o»i;A.'t;12=*.r‘#.t»f2X t* NUl = • t .fc f t X # * NU2= • . f 9 .6 
SF 9 . 6 » 2 Xf •bM 2 =» tfV. 6 » 2 X I 


fURMAfi lx, ‘ti Is • jHV.to ,2X 
♦»2X.»GM1**,F 
»<Rl 1 e( 6 ,lU 0 ) 


SYS.. • ,3X» *V2:avtKALL 1 .F . • ».^X , M Y 1-Y2 )4»2 


FUKHATi lx, *Yl rCwNiKLLLLU 

4* ) 

DU 3b J=*l 1 1 00 

Y»(04A2l4AU)*Yl-|b*i A2l4All)-A224All*A2l-Al2l*V2-i AU4A224A124A21 
♦-b*i-A>224All4A.2l-Al2 ) )♦¥*-! Al 24A22-&4IA 1 1 *A224A 124A2U l♦V9♦G*Al2♦ 
4A22»Yb4 ll.O-O )*ibll43tl )*Ut*l I .0-0 1 ♦ l-B 1 1*A2 l4L*22 4b 12-02 l*Al I j*U3* 
♦i I .0-0)4l-fcll»A22-al24A2l-b2lAAl2-b22yA ll l*U94| 1 . 1 - 0 ) ♦ I -B 1 24A22-B2 
♦ 2 *Al 2 )*Ub 

UsD*K4En*0i4tl2*U24ANUi4Yi«AN02*Y2 

Yb = Y9 

Y4 = Y3 

Y2*V2 

Y2=Yl 

VlsV 

Ub=OH 

0*.=-w3 
U3 = U2 

U 4 SUI 

U I S‘0' 




D-4 




Y ^ * -Au* V 1 1 “/ V Y -A * yt > -A J <*r t, **-A c’ * V (- ^- A I * Vt t>-A D* VC 7 ♦ t t)* K i ♦HA *R 3 ♦H 

VC7»fC6 

VC(>«Ytt> 

VLt«VCA 
YCasYCa 
Y C3*YCi» 

YCilsYCl 

YU*YC 

K7*Ko 

K*>»Rh 

KA*Ki 

K3«H^ 

k2»HI « 

Kl*R 

uV-iV-Yw)**^ 


WKlU(6t&t>)Yt YtfUY 

PUUMaTJ lA»MA.s,tVX,hl.A.J,bX|F i3,d#‘jX) 


♦ u«n#» 111. 
L^JlUMuU 


ftiR 
kNU 


at»« 

0.03 V A*tA 

kA- 

C.UAt /3t 

b3» u.(j*«a3{}0 

l*i = 

fc iK 

0.U0( (iOU 


( .1 O0u(jC 



A7* 

l.OCuuOJ 

A6 = 

-«J. VVUi'Oi 

At *•"'1.' • ? t>4 

#.<♦ » 

A^s 

t'.VAOLl'* 


Al»“l.’.U4,.bbd‘J> 

Al-s-< 

IM 1 i 

t l.fcVlOl 


.--l .7^u*.ilb 



Ell 

e-U. iiX't»6 / 

i. t 

L- --U.U /IAI (• 

NUl «-U.M 

»' i t N 




Yl :tUNTKLLLLw bil . 

Y.:LVLK7UL T.b, 

1YI-Y2JA*2 

C.UUOOUUOi/ 

(..Ot'CCOUOO 

0 , OuO 0 0 0 0‘J 

0.( (jtOl CUv. 

( ,t( tl wuLO 

0 ,000( 00 ( (j 

O.OiVifA^A'* 

1 ,C J /aaaa*^ 

0 .OoOUOUOO 

u . n AO I t^t 

». , 1 A >' 0 1 |j U 

O.OUOOgUUO 

0. 4,Hb6tOUO 

0 . 2 •• 0 0 1 0 1. 1 > 

0 .OliUOOOOO 

0. AAb^iblOO 

L .AAO Jb loU 

0 .OUOOOOgO 

o.oovtjV- ic 

f • 0 6 Vo 7 a (J (1 

3 • V JO 00 0 ^ >• 

C. va il tCbt 

t .V diC 1 c 0 1 0 

O.tOOl.OOt ( 

1. i'mtAO j 

I.IVU3300 

o.oyooeouu 

l.Ai'rgtOOU 

1 »**3vy**VUD 

U.OOO'OOOOd 

1 . t) V UV dVO 0 

1 .b7CVb6Uu 

0 .OUOOOOOO 

l.dtt IbAOO 

l.dt3l(32C’U 

O.UOOOUOOO 

2 . tl ^'jodAOC' 

2 .0 106320C 

0 ,( C OOOOOu 

2.U (’i i'cOO 

2 . i I 7373UU 

o.couoooou 

2,13tjk:l>0u 

2.lbb32iOU 

0 ,00006013' 

2.;2blHbJ0 

. 4 2 1 1 3 0 0 0 

0. 000 C 60 00 

2,22b7ttOO 

2. .,07 3900 

O.OUO jOOUU 

2.222^‘#70(l 

2 2230 100 

0 .OOdOOOOu 

2.l9c9t30b 

2 .1 9LV2 700 

0 • t-OO'OOOOO 

2. l7l2tV0C 

;:.i /i.^6A00 

o.oocuoooo 

2. lAVIddOU 

2.1AVib3Q0 

0.0000 JO OKI 

2. l2A'^C'4.0t 

2. l.AVObUU 

O.OOOOGoCO 

2 • 1 1 to f to U 

2,1 l3 0d7U0 

0 • 005)000 UO 

^.IOVhVVOI 

2 ,U'9‘, 7voo 

0 ,ooo;c juoo 

2.112l2bOi.> 

1213200 

0.600000(0 

2. l2t^l50U 

2.123..17UU 

O.UOOO'JOOO 

2.13777700 

2. 137/3300 

o.ouuoooou 

2. 13072900 

2,1 3373300 

0 ,000 OUO ( 0 

2. 177a370o 

2 ,17 /AAbUO 

O.OOOUUOUO 

2.193HC30U 

,. , 1 93a 1 600 

O.tOOl OOlOl 

2.212VA30C 

2,2 13VAOOU 

0. coo COO 00 

2.22A97«0U 

2,22AVV100 

O.OOUCOUOU 

2.2303t bOO 

A .2363 ?20O 

0. OUO too ('6 

2.2A2UA00U 

2 , 2A tUbAuU 

0. coo 000 00 

2.2AA2 2HH 

2. 2 ‘♦‘♦2 A 000 

C.GOVOUoUSO 

2..,a302oOO 

2 ,2A3t'AU00 

0.0 (,* 0 ( 0 ( ,* G 

2.2A.;0lt0v 

2. 2 A2 03100 

0. coo GOO 00 

Z.ZAUilUO 

• 

2,2Au l2tOU 

• 

u.(;uo( JO to 

• 

t 

2. c32fc 3t0o 

• 

2,2 0. bAAO 0 

, 

O.OOuGOOv.'U 

2.23^db70U 

2.23gbVbOU 

0 ,0'UOOUUI.'U 

2. 2 3cbo ?0C 

2 ,23^6/300 

0 .OUUDgt f’O 

2.td2b 

,2 Aiit'AiOC 

(1.0’CO l.O'U •, 








D-5 


u#.1A 1/U.!«/tAL/UO.i.i>/,iifc1/u.^»U.(>«V.U/a.!i.S.0/tU/0.5/.Z{:TA/0.fr/ 
UAlA VI tV/tVStV^t y&»UAti:<t»U*itUAtUf»/U»0t0.4)«Ci.0t0»0t0.0t0»(»t0*0t&*0 

^ I mit/ 

UAIA Vi.ltVc2tVi.itVi^fVl.dtVt6«VC7yu.UtO.UtO*0«0*0»p.OtO.OtO.O/»RlfR 
♦ 4, * AA (K 1> t R** • R 7* i/*C ttf*OfO*U*Ci«0#v»0#0«0#v»0# Z* C/ 

l^lfau.l jA^i»bm> 






(»/ms--kx»»i-2*ztTi I 11*1 1 


}{: 


• kb.l IAUeu<;4 II 
iW.;’!>A^2>o.:i2) 


17 
171 

« 14 11>»».^))^U7I Ul/SUKI I Wil II )«SINlM|i)»1*SQR14l*2Cl41)**2)) 

♦)) . . 

72l!)«4Aiil »/Wut*i2»»4iX^I-27lll I*WI1I*1 )• lfcX7*l-lliT 1 1 )*NI I )*1 IhCO 
«:>lki4il*UbWKl4l-/.LI4U«*21)*4ZLlll)/SQRm-^ZCT|l)**2)|*SII»lWllUT* 
«bukMJ-XLl411»*<;l))l 
IFll.iQ.l )t>i/-724 1 1 
l741.kg.^ll.22s^F<:44.) 

K iUNllNUi 

b7a»< *fc A7 4-27TAAM4 1 1«1 l«tUb4M4 i 4 1 .O-Ze i A**2) I 

bh2«->tX74-<.»2t1A‘»w4 n« 1 1 
I ll^bMl-AU 

mNU 1^4L ll*Ai.i« li.{ l*All->Al2«i^H2 t»4Al l»Al2*bl 1/B12 ) 1/ 4 bl l*AI 1«B12-A 1 

♦ <.<4:ll**i/bl21 

L 1 c ~A 1 m 1 1 -41 1 1 *AMU l-A I 

ANUb«Ll£*Al«./b)4. 

t,s|X7 4-i.lC*11 

kH *l>* 4 1 *U~u I * 4 “hi I *ki'i ♦! 1 2-»ji*l* Al 1 ♦b22» 

b-i^l»*4l .0 -ol*4 “bil*A<t4C~bi/*A2l-b2l*AI2-»2 2*Alll 

I *2 4 1 . t - - 1* i -b I />» A*'4' -b2 2* A1 2 I 

l>t>O.U 

bu»b I 

A i K 1 . u 

At.«-A2i-A II “i» 

Al> *“A^2 *A 1 1 *A 1 -A 1 i ♦b'*" 4A2l4Alll~fcli*4 1»0~GI 

A-«>AU»H24.4Al<;«A<.l>b«4-A«l24All«A«.l-A12)^k11*ll.G-G)*4A2i4All)>fe12* 
« i I .C'-b>-AN0l»4 l.u~GI'»4ull«l»2i 1 

AisAl2«A^ . ~b*4 AAl%Ak**«Ai24A2i 1 “tl I* 4 1 .0-G I* l-A ?2*A1 1*A21-A12 » -€ 12* 

♦ U .0~b»«'4-A/l-AU 1-ANUIAI I.U-GI*l-GII*A2l4B2A«B12-ll^^l*Aii)-AMU2*II 

♦ • L -i» I ♦ ( bi 1 ♦ bn , i I 

AA*-b4AU-«Ai^-t 1i*4 l.u-4.|*|Ali*A22*Al2*A2ll-6T2*ll.O-GI*l-A22AAU* 
«ALl-Al^MMNUl«11.0-b)*44ill«A22«bl2*A2l«b21*A12«^G22«All )-ANU2*II.Q- 
•4-bl l*AAl+b*.'-f*t I A*bkl*All 1 

il«41.b>b)4AlA»A22-l.)^>»4l.U>GI*IAll*A224AI2*A2n4^ANU2*4l.O-Gl 
**ttiiUAiA4blA»‘AAl *l)<l<»Al2*B2*lAAll )*ANU1*4 ♦! li 12*A22^B22*A12 I 

A0=-11a^4 i .b-o*^Al2«A22*AMU2*i l.U-Gl*lbl2*A22*B22»Al2l 
wKlIk tk itvib::>«bAtB39b2»bl»tiU 

tv. 7UKMAmA,»Bbs».7y.ttAX,*BA^*,7V.6,2X#»l43«»t7 9.6,2X,»B2«*tF9.tot2At 

♦ •bl=-*t7V,fc,.^Xt •7t» = ‘ t7V,6#2Al 
WKIH 4tt7l)lAlt At>« A b t AA t A.7 tA«.'tAlfAO 

/b buKkAl 4 IX * 

«*A3-»,7V 


X>»A/s»,kV.6t2X»»A6=*,7V.fc,2X*«Ai»««fF9.6,2X,*AA»**7».6f2A, 
,Z,2X. «AZs:*.7y.tt2X.*Als*|7y.6t2Xt*A0»*t7 9.6»2Xl 

A . f c ir I . * 1 A.tk i A 1 k 1 Ai. k ''k 


3L 


i bb 


WRl It 46 tioUl I »tr2tANUI»ANU2tGMItGH2 


.6t2Xf«NUl»S79.6»2Xf *NU2>^**79.6 


lUKhAU lX,*4.1i^*,l V.tt4X|*L12f •|7Y 
1 2 X t * 4»H 1— * 1 79 »6 f X t *(»M2 *■ fr9«6t<tX I 

7uKjAV4Ui*{l:bUNTRUlL4.D SVS** t 3Xt •y2tOVfcRAlL 1.7 . • t iX # * I VI-V2 1 ♦♦2 

♦ • 1 

bU Jt J=1*I0U 

V^46*A2HAlll*Vl-464lX*lkAli l-A22«All*A 21 -A 12 l*Y2-l All* A22« A1 2*A21 
♦-6*4-A22^AU=«'A2l“Al2»)*V3-|Al2*A22-G*4All*A224Al2*A21i 1*YA^GAA12* 
♦A^:2»Yb*4l.O-b»‘*4Bil^B2llkU2^II.O-&l*4-4Jll*A2l*B22*Bl2^2l*Ain*U34 

♦ 4 l,b“G)*4 -Gil *A4:*.-bl<^*A2l-BAl*Al2“622kAll l*G44| 1.4J-G1*! -B I2*A22-«2 

* U^b*A-^**G»»Hl*ANUl*Y 4- ANUI*G*V2*ANU2AV2-AN4J2*li*Vj»*&*UI*fe 11*1 1.0-61* 
iU-cH ')2*^4 1 

Yb-VA 

VA-VJ 

Y3=Vt 

VZ-Y4 

Yi=Y 

Ub-UA 

0A=U3 

0 “U *C 

U1 s-k> 




D-6 



Jlt 


Vtl«VCt» 

VC6»¥t6 

Vt5«-VC4 

VC4»Vti 

Vt3»VC4. 

\rC^»Vtl 

VCi»VC 

K|«A4» 

Rc»»k& 

H4S.A3 

K3*R4i 

RlaK 

UV»IV-¥C>**i 

kKltbU>.(.*>)YtVC«UV 

HlHkAl(lAtFi3.U,9X,f^l3«D,!»Xfl IB.bt&XI 
CUNI InUfc 
&10P 
I NO 


bl»* d3* O.U<i4ifil» bi* U*0324»62 

0*000000 00 * 0*000000 

l.tuiouu A6KoO*V2bbdca At>>>0.744*»70 A4* 0.360U2 

»3. 0.M6T0T «2>KI.»»;7IT *|. O.C24S«3 «0>^.0t3201 

OMl* l»O0?101 (»N2«->0* 740810 

t tl4'*-C.Ci#l46fa Mll»>U.4»71>2li3 NU^» 0.!»0&740 


vislunihuluo bVS. 

i'.uOUUUUUO 
u.UUOCCuOU 
u. 037^444^ 

I 

g.26t>9lUA0 
u.^tbAi/Sbu 
0*?32U<»loO 
i.uli»t7i<uo 
I.2927bdU0 
i*!>!>0V3bC0 
l*Vb2A!>800 
i«scO:.'Ati>o 
< .lC'239b00 
2* 1 V3^VC'wO 

«..2bV 17900 
2.26iVt>i>UU 
<^»iA7761?00 
3..2221b^UU 
2* l9at)*>7oO 
2*171^0100 
2* A*>OVb2O0 
Z.lb<*lb200 
• il>3 33200 
2.1 bo 4 loot) 
^•17140600 


V2SUVLRALL T.b. 
0.00000000 
O.OUOoOOOO 
0.03744444 
o.iloybzdo 
0.26691040 
0.4841I480 
0.73204660 
1*01607700 
l*2921bb00 
1*66093200 
1*78246000 
1.96063200 
2.10238100 
2.19346800 
^•24463400 
2.26916200 
2.26192700 
2.24773600 
2.22216100 
2.19663200 
2.17726600 
2.16891 300 


:.16411300 


2.16329200 

2.16836800 

2.17136700 


IYI-Y2I462 

U.OOOOOOOO 
0.00000000 
0.00000000 
0.00000000 
0.00000000 
0.00000000 
0 .00000000 
0.00000000 
0.00000000 

8 : 81888 Sg 8 

0.00000000 

6.00000000 

0*00000000 

0*00600000 

0.00000000 

0*00000000 

o.oouooooo 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.06000000 

0.00000000 

0.00000000 

0.00000000 


« 


# 

« 


2.23287vu6 
2.23287506 
2.232619U0 
2.232 91600 
<.'•232 6220 c 
4.3:3267600 
2.23267300 
2.23262666 
4. 432&9800 
2.23282400 
4.23206700 
2 .23266600 
Z...3262800 


2.23284000 
2.23283600 
2.23277600 
2.23487200 
2.2 3217 700 
2.23283206 
2.23282900 
2.2 3278200 
2.23286600 
2.23278600 
2.434 64900 
2.23283000 
4.23279300 


0.60000000 
0.00060000 
0.00000000 
0.0^00000 
o.oooooooo 
0 .00006000 


0.00000000 

1.00000000 


o.oooooooo 

o.oooooooo 

o.oooooooo 

o.oooooooo 

o.oooooooo 




Appendix E 


Contents: 


Factorization of (3-4i8) and (3-61) for various A 2 iw 
Table 3-1 

p.E-1-2; Sample program for (3-48) for Xj ■ I 

E-3-4; Coefficients of (3-49)-(3-60) and singularities 
of (3-48) for ^2 • 20, and 50 

E-5-6: Sample program for (3-61) for ^2 “ ^ 

E-7-8: Coefficients of (3-62)-(3-68) and singularities 

of (3-61) for ^2 " 10. 20, and 50 


Note: a’s ami b's correspond to p's and q’s In (3-48) and (3-61). 
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Appendix F: Sample Progranui for the Five Adaptive Controllere 

In Section 4 and the Run Statletlca In Section 5. 
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bLOCK DATA 

IMPLICIT RtAL*t» (A-fl,0-SJ) 

HLAL*M 100(10) ,10 (1U) ,NITM (10) ,SASI ( 10) , SAI'I KA ( 10) , 

• SAPPFB(IO) ,SPPKVA(10) ^SPPLVH (10) ,SACrEP ( 10) ,SACPLO(10) , 

• SCPBVF(IO) ,SCPBVG(10) ,MAFTE(10) ,KSTE¥(10) 

INTEGER K0(12) 

COM BOM /STATS/IDO, 10 , K ITf , iJAfel .JiAPPEA , EA H'Eb , 

• SrPBVA,SI*PtVH,RACPBI%SACI>Ea,SCPEVF,aCPEVG,NASTE,NSTKV, 

• S TE.SSPI A, SSPI B, .SSCP P, SSCI‘6, M,KD 

DATA KD/0, 1,2,b, 10,20,50,100,200,500, 1000,2000/ 

END 


IMPLICIT RtAL*H (A-H,0-?,) 

BXTfcUNAL OVPL 

Hfc‘AL*8 A1(5) ,A2(5) ,0 (5) ,K(5) , AMP (4) , N (4, 5) 

RKAL^U C(5) ,D(5) ,AHAT(4) ,DI!AT(4) ,PilAT(4) ,«HAT(4) 
hLAL yniN(4) ,YMAX (4) 

hBAL*0 IDO (10) ,10(10) ,NITfc (10) ,SASI ( 10) , RAPPEA ( 10) , 

♦ f>APPkI»(10) ,SPl’BVA(10) ,SPIi;V» (10) ,SACPKP(10) ,SACPEG(10) 

• SCPKVF(IO) ,SCPRVG(10) ,NA.STK(10) ,NSTFV(10) 
iihhl*a BSTHUL(O) 

iNn;aBr< pctl,ptkm ,1'Xchg,bbtstk(3) 

INTIJGBR KD(12) 

COMMON / ST A T S/ i DO , I U , N 1 T E , S A S 1 , S A I' I' B A , S A I ' P E I) , 

• bPPKV A ,SPPKW li,SACrRK,3ACl 1;G, SCPRV F, SCI'KVG, NASTK, WSTEV , 

* STE,SSPPA,SSPPD,SSCJ P,SSCPG,M,KD 
COMMON /HIN/LP 

COMMON /CO NT E L/F , G , A H A T , D H A T , F HA T , G H A T 
DATA K/4/ 

DATA A 1/0 . 95 DO , U . 5 1'O , U .91)0 , U . 9 DO , 0 . H i DO/ 

DATA A2/U.2D0,0. IDO, O.HDO, 0.71)0, 0.3D0/ 

DATA U/0 .065D0 ,0 . 1DU ,0 .0 BDU ,0 . 5D0 ,-0 . 193 DO/ 

D A TA B/0 . 0 1 DO , -0 . 0 ID 0 , 0 . 0 1 DO , - 0 . 3D0 , 0 . 20 01)0/ 

DATA C/0.4D0,0 .7D0 , 1 .0 DO, 0 .ODO ,0.000/ 

DATA D/0. 0U0,0 .b5D0,0.5D0,U.6D0, 0.700/ 

DATA YMlN/-1.0,-4.0,-4.0,-1.0/ 

DrtTA YMAX/2. 0,4. 0,4. 0,2.0/ 

DATA MAX, NSRBD/1000, 13579/ 

DAI A AMP/ 1.000,2.0 DU, 2.0 DO, 2.0 DO/ 

DATA W(1,1),«(1,2) ,W(1,3) ,W(1,4) ,W(1,5)/ 

• 0. 17851DO,O.344b3DO,U.55452DO,O.4OUfa6DO,O.20534DO/ 

DATA H(2,1) ,H(2,2) ,W (2,3) ,N(2,4) ,M(2,5)/ 

♦ 0.a4629D0,0.0b157D0, 1.38fa29D0, 1.0217D0,0.71335D0/ 

DATA H (3,1) ,W(J,2) ,H (3,3) ,N(3,4) ,H(3,5)/ 

* 0. 1705100,0.3444300,0,5545200,0.4004600,0.2853400/ 

DATA W(4, 1) ,W(4,2) ,W (4,3) ,W(4,4) ,W(4,5)/ 

• 0.44629DO,O.H6157DO, 1 .3862900, 1.0217D0 ,0.71 J35D0/ 

DAI A ESTMUL/O.ODO, l.ODO, 1.200/ 

DATA PCTL,PTEMP/* 1», 'O*/ 

DATA BSTSTR/»-20X» ,* 0 »,»420*»/ 

LP = fa 

CALL ERHSET(207, 1000,-1, 1,UVFL,0) 

WE1”^K(6, 10) 

10 roi<MAT(» PRINTOUT iNTthVAl- ^ •) 
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, WKITR(6»20) 

I 20 FOHnATf* STARTING AND ENDING EXAHRLBS « •) 

1 hEAD(5,*)NEXST,NEXEND 

HR1TB(6,30) 

I 30 rORHAT(* SVARTING AND ENDING ALGORITHMS » •) 

I READ(&,*)NALGST,NALGED 

NSC 2 

. DO 1000 NBX»NEXST«NnXEND 

DO 100 NALU«NAL6ST,,NAtGEl> 

^ DO 100 NBST-1»3 

DO 100 1NP»1,4 

I AHAT(R) A1(NEX)«1:.STRUL(NRST) 

I AHA'i’(K*1) * AHAT(K) 

r 1 AT(K) » B(NEX)*ESTnUL(NBST) 

, UHAT(K-I) » DHAT(K) 

i GIIAT(K-I) * C(NEX)/D (NEX) 

t5HAT(K-1; * GHAT(K~1) 

, «IIAT(K-2) » (illAT(K-l) 

6 e GHAT(K-1J 

‘ PHAT(K-I) > (D(NRX)-AI (N£X))/G (NEX) 

FHAT(K-I) » FHATJX-1) 

1 FllAT(K-2) » fHAT(K«1) 

I r * PHAT(K-I) 

WRITE (10,60) PCTL, NEX, 1NP,ESTSTR (NEST) 

50 K0K«AT(A1,15X, •EXAMPLE*, 12,* INPUT*, 12, • NON-ADAPTIVE * , 

* • BST *,A4/) 

WRITE (6,50) PCTL,NBX,INP,ESTSTB (NEST) 

PXCfIG * PCTL 
\ PCTL * PTEttP 
^ PTEMP = PXCHG 

CALL SinOL(A1(NEX) ,A2(NEX) ,D(NEX) ,E(NEX) ,AMP(1NP) ,»(INP,NKX) , 
I • NSEE0,1NP,NALG,NEX,NSC,C (NEX) ,D(NEX) ,YniN(lNP) ,TRAX(1NP) , 

I * INTV,HAX) 

100 CONTINUE 
» 000 CONTINUE 
I KNDFILE 10 
* STOP 

END 

I SUBROUTINE OVFL 
LOGICAL OVFLOM 
COHKON /BOBB/OVFLOW 
WHITE (6, 10) 

! 10 FORnAT(/» OVERFLOW *♦♦•/) 

OVFLOH .TRUE. 

RETURN 

END 



i 
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BLOCK BATA 

I IMPLICIT BBAL*e (A-H,0-Z) 

K£AL*8 IDO(IO) *10(10) ,N1TL|10) *SAS1 ( 10) , SAPPEA ( 10) * 

I • SArPEB(IO) *sa>EVA(10) *SPPEVIi (10) *SACPEF(10) *SACP£G(10) , 

I * SCPBFP(IO) *SCPEVG(10) *NASTE( 10) *BSrEV( 10) 

INTEGER KD(12) 

. COMMON /STATS/IDO, 10, N1TE*SAS1,SAPI'EA, SAPPED, 

( ♦ SI PEVA,SPPEVB,SACPEF,SACPEG,£C1‘EVF,SCPEVG,NASTE,NSTEV, 

' • STE,SSPPA,SSPPD,SSCPP,SSCP«,M,KD 

DATA KD/0, 1,2,5, 10,20,60, 100,200,600, 1000,2000/ 

END 

IMPLICIT HfcAL*b (A-H,0-Z) 

1 EXTERNAL OVPL 

I LEAL*e A1 (6) ,A2 (6) ,D (6) ,E(6) ,AHP(4) ,U (4, 6) 

REAL*8 C(6) ,D(6) ,AHAT(4) ,BHAT(4) , THAT (4) ,GHAT(4) 

, REAL yniN(4) ,ynAX (4) 

REAL*8 IDO(IO) ,10(10) ,N1TE(10) ,SASI(10),SAPPEA (10) , 

' * SAPPED(IO) ,SPPEVA(10) ,SPPEVD (10) ,SACPEF(10) ,SACPBG(10) , 

« SCPEVF(IO) ,SCPEVG(10) ,NASTE( 10) ,NSTEV( 10) 

NEAL*B ESTHUL(3) 

INTEGER PCTL,PTEHP,PXCHG,ESTSTR(3) ,SCSTR (3) 

INTEGER Kt(12) 

( COMMON /STATS/IDO, 10, NITE,SASI,SAPPBA, SAPPED, 

* SPPEVA,SPPEVD,SACPEF,SACPEG,SCPEVP,SCPEV6,NASTE,NSTEV, 

* STE,SSPPA,SSPPB,SSCPF,SSCPG,H,KD 
COMMON /LUN/LP 

1 COMMON /C0NTBL/P,G,AMAT,BHAT,PI1AT,GHAT 

^ DATA K/4/ 

DATA A1/0.96DO,0.9DO,0.9DO;;0.9DO,0.83DO/ 

I DATA A2/0.2D0,0. lDO,0.bDO,0.7DO,0.3DO/ 

1 data B/0.06bD0,O,lD0,0.O8D0,0.6D0,-O.193D0/ 

DATA E/0.0 1 DO , -0 . 0 IDO , 0 . 02 DO , “0 . 3D0 , 0 . 20 8D0/ 

DATA C/0. 4DO,0.7dO,1 .0DO,0.BD0,0.6D0/ 

I DATA D/0 ,0D0,0 ,66D0,0.6D0,O.6D0, 0.700/ 

* DATA yHlN/-1.0,-4.0,-4.0,-1.0/ 

DATA YNAX/2. 0,4. 0,4. 0,2.0/ 

DATA MAX, NSEt.D/1000, 13679/ 

' DATA AMP/1 .0DO,2.ODO,2.0D0,2.0DO/ 

DATA W (1,1) ,N(1,2) ,W (1,3) ,W(1,4) ,H(1,6)/ 

* 0. 17861DO,0. 344U3D0, 0.5646200,0. 40«GbD0,0.2O534DO/ 

DATA W (2,1) ,H (2,2) ,H (2, 3) , H (2, 4) ,W(2,6)/ 

* O.44629DO,O.0bl67DO, 1.38b29D0, 1.021700,0.7133600/ 

DATA H (3,1) ,W (3,2) ,W (3,3) ,W (3,4) ,H (3,3)/ 

* 0. 1786100,0. 344b3D0, 0.6646200,0. 40i3bbD0, 0.2853400/ 

DATA W(4,1),W(4,2),M(4,3),H(4,4),H(4,5)/ 

* 0 ..4462900,0. 8b167D0, 1. 3862900, 1.021700,0.7133600/ 

I DATA ESTMUL/0.8D0, 1.0D0, 1.200/ 

DATA PCTL,PTEMP/» 1*, 'O*/ 

DATA ESTSTK/*-20%* ,• 0 •,'^20*»/ 

DATA SCSTR/»SC- »,*SCO •/ 

LP = 6 

, CALL ERHSET(207, 1000,-1, 1, OVPL, 0) 

1 HR1TE(6,10) 
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1 10 POHHAT(* PRINTOUT INTERVAL * •) 

* READ fS,*) INTV 

I NRIT£(6,20) 

1. 20 FONHATC STARTING 4N0 ENDING EXAHPLES » M 
READ (5,*) NEXST,NEXiP.ND 
{ NRITE(6,30) 

I 30 PORHAT(* STARTING AND ENDIMO ALGORITHHS •) 

READ (5,«) NALGST,NALGED 
. NR1T£(6,40) 

40 FORBAT(* SflOOTHIHG COEFFICIENT * •) 

‘ READ (ft,*) NSC 

DO 1000 NEX«N£XST,NBXEND 
I DO 100 NALG«NALGST,NALGED 

I DO 100 NEST* 1,3 

DO 100 INP>1,4 

I AHAT(K) » A1(NEX)*BSTMUL(NEST) 

AHAT(K-I) - AHAT(K) 

BHAT(K) * B(NEX) •ESTHUL(NEST) 

DHAT(K-I) s BHAT(K) 

GHAT(K-l) - C<NEX) /B (NEX) 

GHAT(K-I) ~ GHAT(K-l) 

GHAT(K-2) “ GHAT(K-I) 

6 = GHAT(K-I) 

FHAT(K-I) (D(NEX)-AI (NEX))/B (NEX) 

FHAT(K-I) = PHAT{K“1) 

I FHAT(K-2) a PHAT(K-I) 

I F - FUAT (K-1) 

WRITE (10, 50) PCTL,NEX ,1NP,NALG,ESTSTR (NEST) ,SCSTR (NSC) 

50 rORHAT(Al, 15X, 'EXAMPLE*, 12»* INPUT', 12, • ALGORITHM' ,12, 

• • EST » *,A4,* SS2 ',A4/) 

WRITE (&, 50) PCT L,NEX,INP,NALG,E STSTR ( NEST ) , SCST R ( HSCI 
PXCHG = PCTL 
PCTL = PTEMP 
PTEMP = PXCHG 

CALL SinUL (A1 (NEX) ,A2(NEX) ,B (N EX) ,B (NEX) ,AMP(1NP) ,W(INP,NEX) , 

• NSBED,INP,NALG,NEX,NSC,C (NEX) ,D(NBX) ,yBlN(INPJ ,YMAX(INP) , 

• INTV, BAX) 

100 CONTINUE 

iOOO CONTINUE 

ENDPILE 10 

STOP 

END 

SUBROUTINE OVPL 
LOGICAL OVFLOW 
COBMOW /BOMB/OVFLOW 
WRITE (6, 10) 

10 FORMAT (/* ♦♦♦ OVERFLOW ♦•♦*/) 

OVFLOW = .TRUE. 

RETURN 

END 
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SUBROUTINR SIHUL (A 1, A2,B ,E ,ARI> 

« HODB, M ALG , NBX , NSC, C , D, YR IN , Y RAX , INT V ,N AX) 

IRPLICIT BBAL*U (A-H,0-Z) 

BBAL*6 RU,U(4) ,Y(U),H(4) ,5(4) 

NEAL*8 AHAT(4) ,BHAT(4) ,PRAT(4) ,GHAT(4) 

REAL S16(4),TIRB 
INTBGBN SYNDOL(4) 

hBAL*8 IDO (10) ,10(10) ,NITt (10) ,SASI ( 10) , SAPI BA (10) , 

♦ SAPPED (10) ,SPPEVA(10) ,SPPKVB (10) ,SACPBP(10) ,SACPEG(10) , 

♦ SCPEVF(IO) ,SCPEVG(10),NASTB(10),NSTBY(10) 

REAL*8 T1(4) ,T2 (4) ,T3 (4) ,T4 (4) ,TS(4) ,Yb(4) ,T7(4) , 

♦ T8(4),T9(4) 

REAL*8 SC(4,3,2) 

INTEGER KD(12) 

LOGICAL OVPLOH 

CONWON /STATS/IDO, 10, NITE,SASI,SAPPEA, SAPPED, 

♦ SPPEVA,SPPEVP,SACPEF,SACPEG,SCPEVr,SCPEV6,NASTE,NSTEV, 

♦ STK, SSPPA , SSPPD , SS CPF, SSCPG , H, KD 
CUnnON /LUN/LP 

COnnON /BORB/OVFLOH 

COMMON /SIGNAL/K,U,Y,S 

COMMON /CONTRL/P,(?,AHAT,HHAT,FHAT,GHAT 

COMMON /MISC 1/T1»T2,T3 

COMMON /MISC2/T4,T5,T6 

COMMON /M1SCJ/T7,TH,T9 

DATA K/4/ 

DATA SC(1,1,1) ,SC(1,2, 1) , SC ( 1 , J, 1) /-O . 9, -0 .9 b, -0 .9 7/ 

DATA SC (1, 1,2) , SC (1,2, 2) , SC ( 1 , 3, 2) /-O .7, -0 .8 ,-0 .9/ 

DATA SC(2,1,1) , SC (2,2,1) , SC (2, 3, 1) /-() .9, -0 .9 ,«0 . 9/ 

DATA SC (2, 1,2) ,SC(2,2,2) ,SC (2 , 3, 2) /-O . 65 ,-0 . 6b ,-0 . 65/ 

DATA SC (3, 1,1) ,SC(3,2,1) , SC (3, 3, 1) /-O .9, *0 .9,-0 . 9/ 

DATA SC(3, 1,2) ,SC(J,2,2) ,SC(3, 3, 2) /-O. 5, -0.5, -0.5/ 

DATA SC(4, 1, 1) ,SC(4,2, 1) ,SC(4, 3, 1)/-0.8,-0.9,-0.95/ 

DATA SC (4, 1,2) , SC (4,2,2) , SC (4 , 3, 2) /-O .4, -0 .6 ,-0 . 8/ 

DATA SYMDOL/'U •, *K», »S», *V V 

OVPLOW = .FALSE. 


DO 100 1=1,4 

U (1) 

= O.ODO 

Y (I) 

= O.ODO 

H (I) 

= O.ODO 

S (I) 

= O.ODO 

T1(1) 

= O.ODO 

T2(I) 

= 0.000 

T3(I) 

= O.ODO 

T4 (I) 

= O.ODO 

T6{1) 

= O.ODO 

T6 (I) 

= O.ODO 

T7(l) 

= O.ODO 

TU(I) 

= O.ODO 

T9(l) 

= O.ODO 

100 CONTINUE 


NSD * NSEED 
M = 1 
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DO 200 X*^1,10 
i IDO(l) » O.ODO 

10(1) » O.ODO 

I SASI(I) « O.ODO 

i NASTB(I) « O.ODO 

NSTKV(I) » O.ODO 
, N1T£(1) O.ODO 

SAPPBACD « O.ODO 

‘ SAPPED (1) * O.ODO 

SACPEP(I) « O.ODO 

I SACPE6(1) « O.ODO 

I SPPEVA(I) « O.ODO 

SPPEVD(l) » O.ODO 

SCPEVP(l) » O.ODO 

SCPEV6(1) » O.ODO 

200 CONTINUE 
KOUNT » 0 
STE » O.ODO 
' SSPPA « O.ODO 

SSPPB • O.ODO 
I SSCPP « O.ODO 

SSCPG » O.ODO 
CY1 A1^A2 
j CY2 * -A1*A2 

COl » B4E 

CU2 « -P*A2-E*A1 

, IPBT = 0 

INIT = 1 

nU = 1.0D0 

I lUlO * 1.0D0 

I ’ 

H « 1.0D0 

j 1F(NALG.E0.3) 0 « SC (NBX ,N?C, 1) 

I 1P(NALG.EQ.5) 0 *- SC (N EX , NSC, 2) 

j DO bOOO N=1,HAX 

IF(OVFLaW) GOTO 6000 

Y (K) * CY1*Y (K-1)^CY2*Y(K-2)4C01*0(K-1)^CU2*U(K-2) 
S (K) « C*H (K-1) ♦D*S(K-1) 

: COHHENT OUT THIS COfiPUTED GOTO AND ALSO THE 

FOLLOWING CALL'S FOR THE HON-ADAPTIVE CASE 

GOTO (210,220,230,240,2b0) ,NALG 
i : 210 CALL ADAPT1 (C,D,MU,RHO,n) 

’ J GOTO 300 

C 220 CALL ADAPT2(C,D,HU,RH0,H) 
j ; GOTO 300 

I : 230 CALL ADAPT3(C,D,nU,RH0,H,Q) 

C GOTO 300 

j: 240 CALL ADAPT4(C,D,MU,HH0,H) 
i : GOTO 300 



2S0 CALL ADAPT5(C,D,RU,RH0,H«C>) 
iOO CONTINUE 

CALL lNP(nODE,N,Anr,H,NSD,K(K) ) 

0(K) » G*R(K)*P*Y(K) 

TIRE PLOAT(N) 

1 

I CALL STAT(N,KOUNT) 

1P((N/1NTV)*1NTV.NE.N) GOTO 10U0 
SIG(1> « SNGL(U(K)) 

S1G(2) » SNGL(R(K) ) 

S1G(3) • SN6L(S(K) ) 

! S1G(4) » SNGL(Y(K)) 

CALL TTPLOT(TIRE,S16,4,SyHBOL,yHIN,YRAX,lPRT,lNIT) 
j 000 CONTINUE 

I R(K-2) » R(R-1) 

B(K-1) e R(K) 

, U (K-2) * U (K-1) 

I U(K-1) » U(K) 

' Y(K-2) » Y(K-1) 

Y(K-1) » Y(K) 

I S(K-2) = S(K-1) 

f S(K“1) = S(K) 

bOOO CONTINUE 
) GOTO 6500 

I >000 CONTINUE 

IDO(R) = O.ODO 
I 10(N) « O.ODO 

SASI(M) * O.ODO 
' N^STE(M) = O.ODO 

NSTRV (H) = O.ODO 
I NITE(R) = O.ODO 

I SAPPEA (H) = O.ODO 

SAPPEB(H) = O.ODO 
I SACPEF (H) * O.ODO 

I SACPEG(R) - O.ODO 

SPPKVA(H) « O.ODO 
j SPPEVB(H) - O.ODO 

SCPBVF(H) » O.ODO 
' SCPEVG(H) = O.ODO 

6500 CONTINUE 

I HR1TE(10,7000) (KD(1) ,X«2,6) 

1 7000 F0RHAT(11X,» K • , 6X ,5 (2X , I5,4X) /) 

MlvITE(10,7010) (IDO (1) ,1*1,5) 

I 7010 P0BnAT(11X,» IDO*,4X,5(1PD11.3)) 

NBITK (10,7020) (10(1) ,1*1,5) 

7020 F0RttAT(1lX,» 10* , 5X, 5 ( IPDl 1 .3) ) 

HR1TE(10,7030) (NIT£(1) ,1*1,5) 

7030 P0RnAT(1lX,» NITE • , JX ,5 ( IPD1 1 . 3) ) 

' WRITE (10,7040) (WASTE (I) , I* 1, 5) 

7040 F0RnAT(11X,» N ASTE • , 2X ,5 ( 1PD1 1 .3) ) 

} ■ HR1TE(10,7050) (NSTEV (I) ,1*1,5) 

i 7050 F0HnAT(11X,» NSTEV • , 2X , 6 ( 1PD 1 1 .3) ) 
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WKITB(10,7060) (SASl(X) ,1 = 1,5) 

:)60 r0UMAT(11X,* SASI*,3X,5(1PD11.3) ) 

COnnEHT OUT TIIBSK WHITE'S FOR THE NON-ADAPTIVE 
CASE BECAUSE THE PARAHETERS DU NOT CHANGE 

WH1TE(10,707U) (SACPEF(I) ,1 = 1,5) 

)7U FORMAT (1 IX,* SACPE F » , 5 ( 1PD1 1 . 3) ) 

HR1TE(10,7080) (SACJ'EG(I) ,1=1,5) 

)80 FORMAT (11X,* SACPE G » , 5 < 1PD11 . 3) ) 

HhITE(10,70y0) (SCPtVF(I) ,1 = 1,5) 

'90 FOHMAT(11X,* SCPEV F • ,5 ( 1PD1 1 , 3) ) 

HR1TB(10,7100) (SCPEV6(1) ,1 = 1,5) 

100 FORHAT(1lX,* SCPEV G • , 5 ( 1PD11 . 3) ) 

1F(NALG.GT.3) GOTO 75U0 
WRITE (10,7 110) (SAPPEA (I) ,1 = 1,5) 

110 FORMAT (1 IX,* SAPPE A * , 5 ( IPOl 1 . 3) ) 

HH1TE(10,7120) (SAPPEB(I) ,1 = 1,5) 

120 FORMAT (I IX,* SAPPE fl * ,5 ( 1PD11 . 3) ) 

WH1TE(10,7130) (SPPEVA(l) ,1 = 1,5) 

IJO P0HMAT(11X,* SPPEV A*,5(1PU11,3) ) 

WRITE(10,7140) (SPPEVH(I) ,1 = 1,5) 

140 FORMAT (1 IX,* .SPPEV B • , 5 ( 1PD 1 1 . 3) /) 

GOTO 7800 

tOO CONTINUE 

WRITE (10,7750) 

50 FORMAT (////) 

00 CONTINUE 

HRITE(10, 7000) (KD(I) ,1=7,11) 

WRITE (10,7010) (IIH) (1) ,1 = 6, 10) 

WR1TE(10,7020) (10(1) ,1=6, 10) 

WRITE (10, 7030) (NlTE(l) ,1=6, 10) 

WR1TE(10,7040) (NASTE (1) ,1 = 6, 10) 

KR1TK(10,7050) (NSTEV (1) ,1 = 6, 10) 

WKITE(10,7060) (SASl (1) ,1 =6 , 10) 

COMMENT OUT THESE WRITE STATEMENTS FOR THE 
NON-ADAPTIVE CASE, SINCE F,G,A,b DO NOT VABV 

WRlTb(10,7070) (SACPEP (1) ,1 = 6, 10) 

WRITE (10,7080) (SACPEG(I) ,1=6,10) 

WR1TE(10,7090) (SCPEVF(l) ,1=6,10) 

WRITE (10,7100) (SCPEVG (I) ,1=0,10) 

1F(NALG.GT.3) GOTO 8000 
WKITE(10,7110) (SAPPEA (I) ,1=6,10) 

WRITE(10,7120) (SAPlKB(i) ,1 = 6, 1U) 

WR1TE(10,7130) (SPPEVA (I) ,1=6, 10) 

WRITE (10, 7 140) (SPPEV B(l) ,1=6,10) 

RETURN 
)0U CONTINUE 

WFITE(IU,7750) 

RETURN 

END 



I SUBhOUTlME INPCnOD£,K«AnP,W,NSKCD,U) 

I inPLIClT KEAL*» (A'H.O-Z) 

REAL RANMO 

, IP(K.LT.O) GOTO 100U 

GOTO ( 100 , 200 , 300 , 400 ), RODE 
* 100 U « ARP 

RETURN 

I >00 CONTINUE 

{ iOO U « AHP*DS1H (H*PLOAT (K) ) 

RETURN 

j 400 CALL RANOU (NKBED,H,RANNO) 

NSBED » H 

U *• ANP*(RANNO-O.2b) 

RETURN 

j 000 U s O.ODO 
RETURN 
END 
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SUBBOUTXNG STATlK^KOUNT) 

IMPLICIT RKAL^B (A-H,0-56) 

R£AL*e R|4),U(4),y <4),S(4) 

HI!AL*8 XDO(IO) * 10 ( 10) , NITE ( 10) ,SAS1 ( 10) , SAPPEA (10) , 

4 SAPPBB(IO) ,SPPKVA(10) ,SPPEV» (10) tSACPEP(IO) ,SACPKG(10) , 

* SCPt!VF(10) ,SCPEVG(10) »NASTE(10) ,liSTBV(10) 

RBAL*8 AHAT(4) ,BHAT(4) ,FliAT(4) ,6HAT(4) 

1NTKGB8 KD(12) 

COnnOM /ST ATS/IDU , 10 , N ITE , SA SI ,S APPE A » SA PPEB , 

4 S PI'BVA , SPPEVB, SACPEE , SAC PEG , SCPEVF , SCPEVG , BASTE, NSTEV , 

♦ STB,SSPPA,SSI'PB,SSCPF,SSCP6,H,KD 
COHHOB /C0BTRL/P,G,AHAT,BUAT,PHAT,6HAT 
COnnON /SlGBAL/R,U,y ,S 

DATA R/4/ 

1F(DABS(S(K) ) .LT.1.0n*03) GOTO 100 
TE * (S (K) -y (K) )/S (K) 

GOTO 200 
100 TE * 0.0 

KOUNT KOUNT*1 
200 BASTE (M) » BASTE (M) ♦TE**2 
STE = STE+TE 

SASl(n) > SASl (N) 4U (K) «*2 
SAPPEA(B) » SAPPEA (H) ♦AHAT (K-1) 

SAPl*EB(n) * SAPPKB («) ♦BHAT (K-1) 

SACP£P(H) » SACPEP(M)*P 
SACPEG(H) « SACPEG («)♦(! 

SSPPA = SSPPA*AHAT(K“1) ♦♦2 
SSPPB » SSPl >P*M!AT (K-1)**2 
SSCPP = SSCPP4P**2 
SbCPG " sscpg-m;**2 

300 TIHE = FLOAT (N) 

1P(N.LT.KD(B*1)) BETUUB 
IDO(H) = S(K) 

10(M) » y (K) 

RANGE FLOAT (KD(M+1) “KU(H)) 

DIV = (RANGE-FLOAT (KOUNT) ) ♦(RANGE“FLOAT(ROUNT) -1) 

KN6 > HANGE'FLOAT (KOUNT) 

NSTEV (H) » 0.0 

IP (DIV. GT. 0.0) NSTEV (M) = (ENG ♦N ASTB (H) -STE**2) /DI V 

IP (NSTEV («) .LT.O.ODO) NSTEV (M) »0.0D0 

NITB(n) * TE 

DIVSOR * RANGE* (RANGE-1) 

IP (DIVSOR.LE.0.0) GOTO 1000 

SPPEVA(N) * (RANGE*SSPPA-SAPPEA(n) **2) /DIVSOR 
IP(SPP£VA(H) .LT.O.ODO) SPPEVA(n) * 0.0D0 
SPi>BVB(B) * (RANGE*SSPPB-SAPPEP(H)**2) /DIVSOR 
IP (SPPEVB(H) .LT.O.ODO) SPPEVB(H) »0.0D0 
SCPKVF(H) = (RANGE*SSCPF-SACPEF(n) **2) /DIVSOR 
lF(SCPEVF(n) .LT.O.ODO) SCPEVF(fl) » O.ODO 
SCPEVG (H) - (RANGE*SSCP6-SACPE6(H) ♦♦2) /DIVSOR 
IP(SCPEVG (M) .LT.O.ODO) SCPEVG(M) «0.0D0 
GOTO 2000 

1000 SPPEVA(B) = O.ODO 
SPPKVB(M) » O.ODO 


SCPKVr(H) » 0.000 
SCI>L'VG(R) » 0.000 
000 COMTINOB 

SA$1(H) » SASI (H)/HAN(;e 

SArOKA(H) » SAOOEA <n)/RANGi; 

SAPPBB(H) » SAPPED (H) /RANGE 

SACPEP(H) « SACPEP(H) /RANGE 

SACPEG(H) * SACPEG (H) /RANGE 

NASTE(n) » NASTE (R) /RANGE 

KOUNT s 0 

STE « O.ODO 

SSPPA « O.ODO 

SSPPR < O.ODO 

SSCPP « O.ODO 

SSCPG » O.ODO 

H » 

RETURN 

END 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

6EX/SEI ALGORITHd 

EOOATIOMS («-4) - (4-0) 

SUBROUTINE ADAPT1 (C,D,KU,R||0,ii) 

IMPLICIT REAL*8 (A-M,0-Z) 

REALMS HU,B(4) ,U(4),Y(4) (4) ,AHAT(4) ,BHAT(4) ,E(4) 

RKAL*e PHAT(4) ,GHAT(4) 

COnnON /SIGNAL/R,U,Y,S 

COMHOM /CONTRL/P,c;,AHAT,BHAT,PHAT,GHAT 

DATA K/4/ 

AHAT(K-2) * AHAT(a-I) 

AHAT(K-I) * AHAf(K) 

MHAT(K-2) BHAT(K-I) 

UHAT(K-I) ® BHAT(K) 

E(K-2) * E(R“1) 

E (K-1) « E (K) 

FHAT(K-2) « PHAT(K-I) 

FHAT(K-I) a PHAT(K) 

UHAT(K-2) » GHAT(K-I) 

GflAT(K-l) = GHAT(K) 

E(K-1) - Y(K-1)-AHAT(K-1)*Y(K»2)-BHAT(K-1)*U(K-2) 

DENOM s H »MU*V (K-2) **24RH0*U (K-2) ♦•2 
AHAT(K) « AHAT(K-1)4MU*y (K-2) ♦E(K-1) /DENOM 
DHAT(K) » DHAT(K-1)^RH0*U(K-2) *E(K-1) /DENOM 
GHAT(K-I) e C/DHAT(K) 

FHAT(K-I) = (D-AHAT(K) )/BHAT(K) 

¥ * PHAT(K-I) 

G » GHAT(K-I) 


RETURN 

END 
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GUI ALGORITHH 


EUUj^TlONS (4-12) - (4-19) 


»UBi^OUTlNE ADAFT2 (C, D,HU ,RHO,H) 
inPLiciT rba:.*8 (a-h,o-z) 

liEAL«8 nU,R(4) ,U(4),Y(4) ,S(4) ,YHAT(4) , AH AT (4) , DHAT (4) 
hEAL*B £(4) ,FHAT(4) ,GHAT(4) 

REAL*8 LAHDDA(4) fGANHA (4) 

COHHON /SlGNAL/KfUfXfS 

COMHON /CONTRL/P,G,AHAT,BHAT,FHAT,GHAT 
COnilON /niSC1/GARnA,LAnB0A,YHAT 
DATA K/4/ 


LAHHDA (K-2) 
LAHDDA (K-1) 
GA«rtA(K-2) 
GAMMA (K-1) 
AHAT(K-2) a 
AHAT (K-1) 


-2» =. 


BHAT (K-2) 
DU AT (K-1) 
YHAT (K-2) 
YHAT (K-1) 
E(K-2) ^ 
E(K-1) = 
FUAT (K-2) 
FUAT (K-1) 
GHAT (K-2) 
GHAT (K-1) 


» LAHBDA(K-I) 
s LAMBDA (K) 

» GAHHA(K-l) 

» GAnnA(K) 
AHAT(K-l) 
AHAT(K) 

DHAT (K-1) 

= DHAT(K) 

YHAT (K-1) 

« YHAT(K) 

E (K-1) 

K(K) 

= FHAT(K-I) 

*- FHAT(K) 

- GHAT (K-1) 

« GHAT(K) 


YHAT (K-1) » AHAT (K-1) •YHAT (K-2) ♦DHAT (K-1) *U (K-2) 
LAMBDA (K-1) * YHAT (K -2) ♦AHAT (K -1) •LAMBDA (R-2) 
GAMMA (K-1) = U (K-2) ♦AHAT(K-I) •GAMMA (K-2) 

E(K-l) - Y (K-1) -YHAT (K-1) 

DENOM * H^MU*LAMBDA (K-1) ••2^RHO^GAMMA (K-1) ••2 
AHAT(K) = AHAT(K-1)^MU^LAMDDA(K-1)^E(K-1)/DENOM 
DHAT(K) - BHAT (K-1) ♦KHO^GAHHA (K-1) •£ (K-1) /DENOM 
GHAT (K-1) C/DHAT(K) 

PHAT(K-l) * (D-AHAT(K) )/BHAT(K) 

F = FHAT(K-l) 

G - GHAT (K-1) 

RETURN 

END 


r>i4 


CCCCCCCCCCCCCCCCCCfCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

SOI ALGORITHH 

fc'OUATIOliS (4-23) - (4-29) 

SUSROUTIME ADAPT3(C,D,nU,KHO,H,0) 

IHPLICIT REAL*H (A-H,0-Z) 

ReAL*P RU,H(4) »U(4),Y(4) ,S(4) «YHAT(4) , AH AT (4) , SHAT (4) 
RFAI.«B V(4) «2(4) , FUAT (4) ,0HAT(4) 

connoN /siGNAL/h,u,y,.s 

COOHON /CO AT ii L/P , G , A H AT , DIt AT , P H AT » GH AT 

COHHON /niSC2/yiiAT,Z,V 

DATA K/4/ 

Z(K-2) * Z(K-1) 

Z (K-1) * Z (K) 

AllAT(K-2) » AHAT(K-l) 

AHAT(K-I) B AHAT(K) 

DHAT(K-2) » BHAT(K-I) 

FHAT(K-I) = FHAT(K) 

YHAT(K-2) » YHAT(K-I) 

YHAT(K-I) » YHAT(K) 

V (R-2) « V (K-1) 

V (K-1) V (K) 

FHAT(K-2) » FHAT(K-I) 

FUAT (K-1) » PHAT(K) 

UHAT(K-2) “ GHAT (K-1) 

GHAT (K-1) = GHAT(K) 

YHAT(K-I) s AHAT(K-I) •Z(K-2) ♦BHAT(K-I) ♦U (K-2) 

DENOn * H4MU*Z(K-2) •*2^ltHO*U (K-2) •♦2 

V(K-1) = (Y (K-I)-YHAT(K-I) ♦0*(Y(K-2)-Z (K-2)))/DEHOM 

AHAT(K) » AHAT(K-1)tHU*Z(K-2) *V(K-1) 

DHAT(K) » DHAT(K-1)4RHO*U(K-2) *V (K-1) 

Z(K-1) * AHAT(K)*Z (K-2)^BHAT(K)*U(K-2) 

GHAT (K-1) - C/BHAT(K) 

FHAT(K-I) * (D-AHAT(K) )/miAT(K) 

P = PHAT(K-I) 

G » GHAT (K-1) 

RETURN 

END 


GKD/SED ALGORITHH 

KOUATIOHS (4-i4) - (4-36) 

SUQROUTiNE ADUT4 (C,D,nU,RU(i,H) 
inpLIClT KL';*t*8 (A-H, 0 -Z) 

R£AL*6 nU,R(4) ,IJ (4),Y(4) ,S (4) , V(4) ,FHAT(4) ,6MAT(4) 
REAL*e AHAT(4) ,6IIAT(4) 

COnnON /SIGItAL/RfUfYtS 

CORMOH /CORThL/r,l?,AHAT,lsllAT,PHAl,KHAT 

DATA K/4/ 


GHAT(K-2) * aiiAT(K>1) 

GHAT(K-I) * GHAT(K) 

FHAT(K-2) » FHAT(K-I) 
fllAT(K-l) * FHAT(K) 

V (R-2) * V (K-1) 

V (K“l) » V (K) 

V(K-1» * C*h (K*2) ♦U^Y (K»2) -y (K-1) 

Dfe,NO« « H*(1.04RU*y (K-2) **2*RHO*U (K-2)**2) 
GHAT(K-I) * GHAT(K-2) ■»RHO*H (K-2) *V (K-l)/DENO« 
FHAT(K-l) » FHAT(K“2) ♦nU^y (K-2)*V (K-1) /DEROH 
G * GHAT(K-l) 

F » FUAT (K-1) 

RETURN 

END 


SOD AtGORITHR 

EQUATIONS (4-40) - (4-44) 

SUDhOUTINE ADAPTS (C,U,nU,RHU,N,Q) 
inPLXClT RBAL*H (A-H,0-2) 

REALMS nU,R(4) ,U(4),y (4) (4) ,V(4) ,PHAT(4) ,GHAT(4) 

RkAL*0 GAnnA(4) •BETA (4) , AH AT (4 ) , DHAT (4 ) 

COnnON /S16NAL/R,U,y «S 
COHflON /CONTRL/F,G,AHAT,DHAT,PHAT,GHAT 
COnnON /niSC3/DBTA»OA»HA«V 
DATA K/4/ 

GHAT(K-2) » ttllAT(K-l) 

GHAT(K-I) « 6HAT(K) 
rHAT(K-2) « FHAT(K-I) 
fHAT(K-l) * PHAT(K) 

BBTA(K-2) » BBTA(K-I) 
bBTA(K-l) ^ DETA(K) 

6AnnA(K-2) > GAnnA(K-i) 

GAHNA(K-I) - 6AnnA(K) 

V (K-2) » V (K-1) 

V (K-1) • V (K) 

BETA (K-1) » KHO*H (K-2) •♦24MU*Y (K-2) **2 

GAMMA (K-1) * (DtQ) ♦H*BBTA (K-2) *V (K-2) ♦D*GAHMA(K-2) 

DENOn » 1.0 + H*DETA (K-1) 

V(K-1) * (S(K-1)-y(K-1)40*(S(K-2)-Y(K-2))-6AHHA(K-1))/DENOn 
GHAT (K-1) * GHAT (K-2) ♦HHO*R (K-2) *V (K-1) 

FHAT(K-I) - FUAT (K-2) ♦NU*y (K-2) *V (K-1) 

G » GHAT (K-1) 

F = FHAT(K-I) 


RETURN 

END 
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RtALY-b Ai 4P4 f At* 4 b) , bl 4 b ) ,l 4. 4 U» ,t lAl 4 b» ft |A^;(b 4 ,(VU1 4 b 4, NU^, 4b 4 
iNUuLk KU4ii-4 
LUtlLAL LiViLLW 
C 

UUMHUN /I) 4 A I b/il.-Lf iOfNi 1 I » bM bi tbAP4 A J f b AP 1. At t b Ap t r- 4 » bAPL It t 

♦ bPtVAl,bf’LVAt#i:PLk/i.-l ,bPt V!.t,bALtll ibACHUtbAU NliS/ulNt, 

* bCkVtltbCtVUAf but YNl » bUtVUt* 

Y NAblb fNbUVtb'ltt 

♦ bbPAl ibbPAPf bbPul »bbPbt ibbCt If bSLEt tbbCNl, SbLlvit,|M|KU 
C 

CUMMUN /LUN/LP 
LOMMUN /uuMu/LV4LUW 
tUMMUN /blOtkAL/KfUf Yfbfi. t 

LUMMON /CJNIKL/AifAtfbi ,1 . i 1 A 1 , II A t* f NU 1 i NU*: t K AC I 

tUMMUN /Mice 1/1 1 fit fi^ f is fib 
CCJMMUN /MiSCL/lOfl /f lb 
C 

DAlA K/b/ 

UAlA bYML UL/»U» t *K* f *b • f • Y*/ 

UAI A bC 4 1 f 1 4 I bu4 1 f / 4/~C .V f ('.O/ 

UAIA be 4<i;il4 »bC4tf«. 4/-1 .tSfjfU.VbJb/ 

UAlA SC ibfli fSCibft 4/“U»Vbfl)*0/ 

c 

LiVI-LUW = .kALbL. 

C 

UO K^C i«i ,b 
0414 ® o.oUu 
0041 4 - u • UL* U 
Y414 - C.UOU 
k 4 i 4 = O.uLc 

sm = u.oui 
11414 = O.Olu 
1/414 “ O . oL u 

1 3 4 1 4 ~ U • V./ u U 

14414 =• O.ClL 
T b 4 A 4 - 0 . L . 0 
164 14 = O.Oi-O 
47414 = O.CtO 
Ibll4 ~ O.uOO 
ioo CUM INUh 

USU^Ubtc.0 
M = 1 

U 0 lie 1 - i k' X o 

100 4 1 4 - o.uL'O 
lull 4 = o.ouu 

SAbl 4 1 4 - 0 .Ok- 0 

t\Ab1 E 4 1 4 - u . uDo 
NbUVll 4 = 4 .OUo 
N1U41) ~ O.ul'u 


r>o o r ron on or o o ooo 
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liO 


!>AI>eAUi) 

SAHfcAil 1) 
t»Aft bl ( 1 > 
l»AH.b4f j i) 
SALLtl ( 1 1 

SACfcNim 
i»ACLN;^< I ) 
!>KLVA1 ( U 
Sf'bVA;: ( I ) 
^.MbVbl ( 1 ) 
iPtVb*.i i ) 
i>ck^vi.i(n 
hCt\/L*H i) 
:>CtvNi j i ) 
SCbVN.: ( i ) 
LUi\li IwL'b 
KUUNJ 
ilt 

SSHA*> 


u .ODU 
U.tilHl 
c • u k*' u 
l.Ol'C 
0 • l< Du 
V •(.*UU 
uJO 
CUuUU 
U.UDU 
C .OLU 
U«l>UU 
o .ulK< 
u*uuu 
u.ul'o 
u • cuu 
O.Ul.t 


o 

U I 
0< 


SSUl 
bbU ii 
bbCM 

i*^K1 - 

it\n > 


.ut o 
.Ul'U 
O.UL 0 

(j * 0 i. b 

U* UL O 
U • Ubu 
U • UiJU 
V • U u u 
0» utl 


U 

i 


MU( n 

MU (<. ) 
KHUJU - 
KHJik!) t 


*» * 


b I > i. 1 /, 

;> n 1 1 4. 


H - i.UU'U 


UU) s 
uU) = 


bU ( iV>v- 1 1 

b 1. 1 f'l w C f 4 . 


LNUR MAIN iiMULAUUW LUUF 


iJU boi/U Ns- It MAX 


ibiuVbLUKi UUlU (OtiU 
V<K* s U.UL'u 


U ALCULA 1 1 (\»l X r bUAN I Do I f U I 
i.'U I bu i s 1 

I'ji) Y(K) - V4KWAU)\>y(K-n»t lU»U(K“n 
LALLULAlt Ulblkt-D lUU'Ul 

1>IK> * UCLlAviht |A( l,i )>!.K(K.--n*bt rA( 1,4 
* ♦UAMM4, t I ) vi> (K-l MUAHMA )vS,( K“< I 




lUMPOU. blAUUJcl. 4»N0 fhlNl VALl'tb 

11 Mb s bLu*0(M 

UALL b I A'U (NtNUlNU 

U- t (N/lNl tf )v INI V.Wt .ivl l.-U hJ K l U 
bib( u s iN^ u (uLik-i n 
blu( A 1 s bNUMK (K~l I ) 
b i o ( 3 » s b Nv. L ( b » K n 
biUltl = bNk-HV(K)) 

CALL I 1 PLt U U Ml , b 1 o ,H , b Y Ml OL t YM IN , V M aX , 1 » K 1 , INI I 1 

bhlbl VAMALLLb lU ALLOA CALLULAllUN Uf NfcW VALUfcb 
iuuu K(K~A» = K(K~j) 
klK~3) = K(N-?) 

k(K~4i) s !v(N-l) 

U(K-<t» s UlK-bl 
U1K~3I = o(N~.‘ ) 

U(K-i) « U(K~1 ) 

UC(K-A) ~ UHK-i» 

UC(k-3l !s uC(K~4 ) 

Ucik-ki) = Utl\-ll 


oo o no o oo ooo o 



Y(K-j) s V(K -.. ) 
YiK-t| s Y(K-l) 
Y(K“1 » - YiN) 

* Mn-.*!) 


S t K*"i; ) * 

- UK-1 ) 

S(K-U a t,(K) 

CPDAIt PLANl HAhAMclU^ 
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Ll>llMAlLi> 


WC’ I Li ( 1 W t i U t t t N A L I 

^10 CALI. bL i ( f i-t^ fKhiifl ) 

bUlU bUO 

0 C ALL GU i t u a 1 1m f I'lk.i ) K nL r I » t 

GUiU jOO 

2.JO GALL bU i t o A^.(•lM f .Mi,, I ,s I i.i» I i t L I 
GflU JVi‘ 

300 CUNllNUL 



HNU Nt h INI Ul K t.L-n 
CALL iN»' iMCuu »N» AMc ♦bbG ,mN~i M 
lUU iMt W LulGlsJCLtK i'i.»liG| 1,‘bJK — l) ANL ACIUATC^< UUJ^'Ui bCK — l ) 
U «NAC1 .Li. . - » GUIJ tlG 

CALCCLA I Un rtClL«liJiv LUt'<Mv.lf t AND A 

CC(K-H - Ui LI A*K(,.-U ♦u l AH K-ll ♦Uv. ( K-2 >♦11 A2tK-i H'UCtK-- J 
♦ ♦WUMN-nvy tK-..)4MU4.lK-n*»'YlK-3l 

1 h I NAw I » L 1 . i I uulvt Atb 
CLIL Ali 



CALCULA I lew huK AeluAlUA LUiNhlb. J» 

Ale UC(K-l) - ul LlA>i‘MK-lHt lAUK-l )*LMlL-*.Ht lAUK-i I ♦GIK-.*) 
» tNUUK-1 )»YlK-i.HNU..(K-l ^♦Y^K-:>} 

Ait» UIK-l) =• t blCl*G<K~;’HLbli.>.>UC(K-id I 
CC10 bOOU 

440 U(K-n =■ LUK-1 ) 
t)000 CGWIlWUL 
GUlU 6b GU 
6000 oUNIKMUL 

IDUM) - o.uUo 
lUlMl = u.uGU 

SAbilKJ - u.GlJG 
NAbItIM) - C.ULU 
NbUViM) - l.OL/G 
WlH(M) = U,ulH< 
bAPLAl I M I - O .Ul'O 
bA^t;A<i(M) - u.uUU 
bAL>tbi(M) = U.UUO 
bAHtbcIM) = O.UUli 
SACclUMJ - L.UL.G 
bACLL-^lM) - U.UI'U 


bACLNilMl * U.L'LKi 
bACtNe(M) - u.uuu 
bPkVAllMl - U.UUU 
bPfcVAtlH) = U.OUU 
SPfcVbllM) = U.OUU 
SPtYbiMM) - u.UUO 
SCfcVLUM) = U.UUU 
bCLVfcirlM) U.UUU 
b C t V f\l 1 i H I ~ i' . U u U 
bCLVtM^IM) =• U.UUU 
cbOO CCM INUL 


WK n t ( I 0 » V UUU M KU U ) , 1 =*' 1 6 » 

7oou POKMATUiX,* K* |7X. bleX# ib,4X)/l 
WRUfc (lUtVuluM luol lit 1-lf b» 

7010 FORMATIllXf* 100* » bXtb UPUil .j I) 


Wkl ltU0»7UM0H ium,is^ltbl 
/OeO PQRMAlUlX,* IL* ,uX,t>< IPul 1. j) I 
«K 1 1 1: U 0 . VOuU M 0 1 I fc ( 1 » t 1 =1 f b ) 

700U PORHAUnx,* NIU* ,4X.bl iPOl i.jn 
WK 11 1 ( 1 Oi 704 0 M iN/ b I U 1 ) , 1 = 1 , b ) 

7040 PDRMAIUIX,' NAblu* .uX.bllPUll .3)1 
WRlUllO.VObuHNbU VU ) ,i = ltb) 

/ObO PUKMATUiXi* Nb U V • , bX t 11 i*U 1 1 . J )) 
WKl I L a 0 , V vi. u M b^ b 1 1 1 ) f 1 - 1 » b 1 


/ooo PUKMAKliX,* bAbi* f-^Xtb (IPUi l.uM 


wk n t 1 1 G » I U V U M b AC L L 1 ( i 1 » i - 1 » b ) 



V04)0 

70 Vw 

7100 

7iwb 

7 UU 

71 ii> 

7120 

mi) 
V l:>o 

71hO 
i Ih'j 
/t»00 

vvt>t 

?bOC' 


oooo 


;»Ai rt 

i 


i'UKhA Hll>- 
rtKlUtlUiVOVt-Mt Au 
► UKMAJUlAf* It 

rtKm ilO,/ubU) tiiMUtWi 


auhma 


up 

i t 


tuxp ^.AUKl" Wi 


I * If Jt 

IPDIi.-P 

f bUHjIl .c » 

^ Ml !?'i 
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ti iU'Oll ..• 1 
I) « 

f ll‘L'l 1 .J » 
U t P 1 1 1> 1 
i:.( lt*un..u 
lyti-itip 

iNUHUll 

IPDI 

I ) 


f'^i 

1 f'0 1 1 . » 
1 ) t Pi »:*1 
t!:«( li'UU. 
liti-i.i.) 
ft>\in>li.j) 
i ) • PI *t>) 
»;>i ii'uii.3) 
i 

Uull.Ji 
*!•( iH'ii .J) 

h,P 

I 


if Pi ft>i 

t li'U 1 1 . 


) 

f 

1 ) 


f P i t L. i * 
( ll'l^ll . J) 
i P i t L i 
( ibuli .j 


WKi It ^Itf (IiAUtiN*. 

i-DRMAPUXi* i»AU‘U 
WK i U I Of >OVO ) I i 01 VL i I 
f-OKMAUlIXf* SCI tv l.i» 
WKlUUUfVuVS) tSCuVu.; 
HIKKAHliXf* SU'tV i 
WKi ItUOf /lOOH jutVivi 
RUKMAl ( 1 IXf • SCPbV 
WRiU lut aot>MSotVN2 
l-UKMAl I llXf* SU-tV 
WKiTUiUf VilJi (SAPuAl 
l-UkMAMUX,* SAH't Ai 
WKI ItIU'f lllSi 
l-UKMAi mXf • SAPKi. A*. 
WKi1t(iUfil<.ui(SAFl:>ti 
FORMAT (llXf» SAHHb Hi 
hknUiOf Ui^SHSAPtuf. 
FUfvHAMllXf* SAPP= fU 
WKi 1 1 1 P<f /i -U) I SPt VAl 
FUKMA I I 1 iXf • SPPL V Al 
WKifldUf/ijhilSPtVAc 
FORMA MiiXP SFPuV A^' 

W R i tell Of M H U I ( b P r. V e I 
FukMAiniXf* SP» tv bl 
WKntlU-f VPSi U» uVl 
FLkMAT 1 1 1 X f • SPPt V 
u>U I U /bvvu 
OUNI llMUt 
WKiltllOfMSOi 
FUKMA I I ////) 
uUM INL L 

wk 1 1 u 1 1 o I < ooo i I i i ) . 1 - < ♦ 1 1 1 
wKilt iiwf /uiu) lil'C. I i If POfiul 
wHlit ilCf /Vm., ul i it II ifi*l flo) 
WK 1 I L I I t* t / u.‘ U I I >1 P t i i I f i = b f 1 
WR 1 U I I t f ,/t ‘*0 ) ( \A5 I I I i ) f i-f I 
rt K i 1 1 I i I f / 0 1> U I I iM b 1 L V I i I f i - 6 f 
wRilLllUi / uuol ll-Abi (1). i-tfi 
wRlItlli^l /uiu) ISAv-tLil 1 )fl-0 
WKlTtliu, fOlbi ISAcut. /ll ) 
wKiTtllOf IUOOMb#<>-l dli) 
wRi 1 1 ( 1 Of /t bs ) I ..Mut Ui I 1 ) 
wk i 1 1 I i I'f i OVO i I 1 ou V M I 1 ) 

WR 1 I L I I Of /ObS M Sot. Vl <■. I 1 I 
WRllLllCf POO ) I SotVWii 1 ) 

W R i 1 1 1 1 0 f Hub i ( S 0 c V i4<. I 1 i 
WR n t li Of / Uo I I SAPl Ail i ) 

WKi It II Of /Pt>) I SAPlAf.1 I i 
wKiU 1 1 Of I sakl Pin 

w R I I L I i U f / i fc !» i I S A H t !' 1 1 t I 
WRP t 1 1 Of P _>UM SPuVAl I 1 ) 
wklltllOf / PSi ISf'UV'Ai;! 1 ) 

WRi I L I i Of / 1<«U ) I bf'LVt'i I A I 
w R i 1 L 1 1 o t P M b i i b F L V i) 1 1 1 i 
RtUiKN 
CLN I iiMO t 
WR A ll I 1 Of PSO I 

K n 0 tt iM 

lNO 


i/I 


A -t. 
i =o 

A -6 
i~t> 

A-0 
1-0 
' -o 

* ~L 
A 

A < 

A-0 

A-O 


i 

I I 

(>i 

i 

ioi 
lui 
10 i 
U i 
loi 
lOi 
lOl 
Alii 
AOi 
loi 
AOi 

iuj 
10 1 
lui 
loi 
loi 


SUbKUOl iffU INPIMUDl f kf AMP ,0SU Of Oi 
AMPLlon koALvi' U-hfJ-*/i 
RcAL KAki><Oll) 

M - i 

AFlK.Ll.O) t.UlU loot 
0 U I 0> I 1 0 0 f «. OO f _ O i’ f K O O ) f M t I,' i 
loo 0 = AMP 
K t T 0 K N 

fCOO U - AMP ‘f'L S i 'i 1 1 ! Li. I O ,4.0’S'f lUA I I N i I ) 
Kt lUkfy 

300 U = AMe »L S Ji'flUi Lt I 0 .S»FLUAI I k i ) i 
KITURN 

HUO CALL oOOl S lObU fMf KA I 
U - Ai.o * (k AoNl i 1 I -V.' . ) 

KLluKN 

lOOU U - O.ULv» 

KfclUKN 

LNb 


iXCCCCCCCCCCCCCtCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCjCCCCCCCCCCC 

GBl/Ui AL&ORIIHM 

feOUAllONS lA-ftdl “ 

SUtiROUliNfc Uil&ANHAfMUfRHQfH) 

IMPLICII KiALA4» U*HtQ-Il 

RLAi«8 tVt&)t&l8).Al(*»l,A2f»)t8H&»ftt2l&|«UU8l 

RfcA4.*i MU(2ltKHO(2)tGANMA(2> 

R2ALA8 M:>I|V(!»I 

HeAL*6 IEIAib>»tlA2lSltNlillS)»NU2IS| 

COMMON /&ACNAL/K»UtVt8tUC 

COMMON /CCNIHL/Al»A2fU»tB2t£lAl»UA2tNUI fNU^tNACI 


C 

c 


c 

c 

c 


c 

c 


L 

c 


c 

c 


OA1A K/t>/ 

&H1M VAKlAbC£S U) ACCOM CACCUCAliON Of- N£N VACUCS 
Al(A>2f AltK'A) 

AI4K-1I « AKKI 
A2fK«2l « A2IK-1 > 

A 2 (K>n > A 4 i(Kf 
Oi iR -21 a BKK-l ) 
bUK-l> a toUK) 
b 2 (K>^f a b 2 IK-i» 

B 2 IR>U a o 2 iKI 


El 

K**C ) 

£ 

EIK^I) 


R-lf 

an 

LiR) 

Vi 

K*2J 

s 

YlK-n 

Vi 

K-U 


WIKI 


iFINACI uulO !»60 

CACCUCAIIUNS FOR ACTUAfOR CONFIC* 1 AND 
felK-l» = VIR«U-Al(R~ll*VIK- 2 l-A^(K* 2 |AV 4 K-il 

♦ >bl IK-I )«UC|K-2)HS2IK>2I«UC{K*3I 
DCNOM a H«MOa>*V(K-2)**2«MUI2l*YiK-81 ♦•2 

♦ ♦Rhul 1 > *00^-2 »♦♦2♦RHUI2 MUt IK-if l♦»2 
^GIC 57U 

CACCUCAIAUNS FOR AClUAlOR CONFIG. 0 AND 4 
J»GO tIK-n a VlK~l)->AllK-UAVIK-2}-A2(K-2)«Y(K-i) 

♦ -bUR-1 *♦0*^-2) -tt2lK-2»AUIK-i» 

DENOM a H«H0m*Y(K-2l**2^MU(2f*Y(K-3f AA2 
, ♦ . «KHUI1 l*U(K-2|A>»^«KH0l2>»U(K-i)«A2 
i>70 VlK~n a c(K>l)/DC'N0M 

CACCUCAfL NtM FCAN1 PAKAMEUR tSIlMAIkS 
A|.IKf s AltK->I>«MUtl)«V(K~l]«Y4K-2l 
A2lKf a A.:tK-U«MUt^)«V(K-ll«Y(N*il 
bllK) a bl(K-l)«KHOrif*YfK>l )«U(K>2) 
b2IK) a b2tK-l )♦RH0i2f «VlK-l M0IK-3> 

OPDAU CUNIKULLbR PARAME1 lR ES11MA1ES 
CALC CNISc f (uAHHAf 
RE I CRN 
END ■ 


CCCCCCCCCCCCCccCl-LCcLLtLCcCs-CCCU-CCcCCCLCCCCCCCCCCCCCCLLLLCL 

c 

C GUI ALGUkllKrt 

C 

C tCUAllUNL - (h-7c) 

C 

i>UbKLJUI ll\L L l i (w>AMMA»Nil. f KhL i h" ) 

IMPLICIT 

kEAL*b Rll>)tHt>)|Yli>TtUc(t>).i^(i.)»yhA|(i») 

KEAL«b A1 tA^ iMtblU )»bc(PT 
H£AL*J> MU21 »kHi)(2) ,umMMA(2) 

kLAL*b L(b),Vlb)fbTAl(;^.),bTA*.(b) ,NU1 ( t-) , NU2( t) ) 

KEAL*b LAMlti.) ,LAM2 «t>l ,bAMl (M ,C.AM2< i>) 

LOMHUiJ /SiGNAL/K tUt YtS »UU 

COMMUivJ /CbNlkL/Al»A<.»Hl yt)«:tblA] iLfAict I »iNlU<.tlYAL I 
CuMMUIM /Ml bCl/GAMi , GAn2 ,LAM1 ,L*^W2^ YHAl 
DATA K/b/ 

C 

C i>HlM VAKiAL.LLb TJ ALLuv. VaLCcS Tu bt CALtULAitu 

LAMl (k“3) a lamJ (<-/ > 

LAMl (K-t.) = LAMl (K~l ) 


r >r r-mrr'O nr on on nn 
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&AHUK-U > OAKUK) 

&AM2U-3I » OAMi'iK-a) 

0AM2(K-U - LAfU iK.» 

AllK-2) s AllK-i) 

A1 IK-1 1 s Al (K 1 
AilK-^l - Ai^lK-l 1 
A2IK-1) » ALIK) 
bllK-2) » bilK-l ) 
bAlK-l) LilK) 
d2Ik- 2) a t)< iK-i ) 
bt(K-l) a ttIM 
VhAllK-^) = VhAIlK-..) 

VHAMK-LI - VhA I IK-1 ) 

VhAl IK-1 ) s Yr A* IK) 
fclK-L) » fclK-i) 
t U-1 ) s t IK) 

iK IWALT .01 ..: ) ov»liJ /oO 

CALCULA I lUis.L t uK AtlUAloK Cu(<» !(.. 1 ANL> : 

YHAHK-l) = A1 IK-l )*VMA I IK-J )^A.ilK-l )-»YhAl IK-/ ) 

* 4bl IK-i ) vuC IK-/ )♦.»/ IK-I ) VOLIK-J) 

OAMllK-1) = tclK-/)*AllK-l )AoAMl lK-/MA/(K-i )*GAM1 lK-:>) 
l»AKi/|K-l ) - Oci K-J) tAl I K-1 )-»OAf',/. lK-2 ) ♦ A2 I K -I ) « (.Ah/ I K-J ) 

GO It VVt 

LALcOLA I lufd AcKmU.K tuN) 1(.. J A.ijU a: 

/OU YHAllK-1) - Al iK-i I^YHA I IK-1 )^A/ IK-1 )»YHAI IK-«; ) 

* ♦ul IK-1 )vU|K-4.)-»b. IK-l)*tlK-j) 

GAMllK-1) = tlK-/)-*MUK-l)VoAMllK-i)-»A<.IK-l)*tAMl(K-j) 
GAMilK-i) = tlK-J)^Ai(K-i l-t'OAMLlK-D^A^ IK-1 )<-tAM/IK-:>) 

VVt LAMlIK-l) - YhAI IK-/ )-»A 1 IK-1 l-^LAMl IK-/)^A2IK-1 lOLAMl IK-3 ) 
LAM2IK-1 ) - YliAl Ik- 3) + A1 lK-1 )*aAM2 IK-..,) ♦ A/ |K-1 )*LAi'u IK-J) 
LIK-I) = V IK-1 )-YhA i lK-1 ) 

OtNJh - l.'H'.UI 1 ) ♦LAHl I K-l ) #♦. •*Mt 12 ) « L .MU I K-1 ) « V/ 

» ♦Khali )«tAhllK-l)v*/4K»TL)l/ )*bAM2|K-l )»♦/ 

lUtlK-l ).Ll .l.uU-lu) tlK-1) - O.ODt 
IFILLNUM.tl.l.OtUi*) ctNah = l.tliU!> 

ViK-i) - t IK-1 )/Uc.NfjM 

CALtOtAlk. t\lt W KLAUI KAK AMu '( c k tlillMAlLi 
AUK) A1 IK-1 )-»MUl 1 )»V lK-1 ) >»LAhl(K-l ) 

A/lK) 5- A..lK-l)4Mt(/)*VlK-i)ALAh/|K-l) 

BIIK) s bl lK-1 UkHU I il’fVlK-l l-K-AMl IK-1) 
eZiK) = ttlK-1 )4kHUl/MVlK-l )»tAM2lK-l) 

CALtUUAlt Ni W CO'hlKULLLk l-’AKAMlUHb 
CALL CNli,/l I OmMMA) 
ht ItKU 
tNt 


) » 


t ^ 


LAMl IK) 
LAM/ IK-/ ) 
LAM2IK-1 ) 
LAM/ |K) 
GAhUK-/') 
GaMI IK-U 


LAMllK-1 

LAMUK-i 

LAM2IK-2 

LAM2IK-1 

OAMUK-A 

GAMlIK-i: 


SOI ALttKllhM 

fctUAlllJNS Ia-76) - Ia-BI) 

bULKUUi INL ^Ui IbAMMAtMUtKHJyHtO ) 

IMKLltll KL L*ii (A-H,t-2) 

kfcAL*b Kii) »U(^) tY I ;>)#:>( ^ItVHAl 15) #ULI5) 
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